Index: /trunk/psLib/src/astro/psCoord.c
===================================================================
--- /trunk/psLib/src/astro/psCoord.c	(revision 4285)
+++ /trunk/psLib/src/astro/psCoord.c	(revision 4286)
@@ -10,6 +10,6 @@
 *  @author GLG, MHPCC
 *
-*  @version $Revision: 1.71 $ $Name: not supported by cvs2svn $
-*  @date $Date: 2005-06-13 20:16:22 $
+*  @version $Revision: 1.72 $ $Name: not supported by cvs2svn $
+*  @date $Date: 2005-06-16 22:00:03 $
 *
 *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -242,8 +242,23 @@
     PS_ASSERT_PTR_NON_NULL(transform->y, NULL);
     PS_ASSERT_PTR_NON_NULL(coords, NULL);
+    printf("psPlaneTransformApply(): (%f, %f) x:(%d, %d)\n",
+           coords->x, coords->y, transform->x->nX, transform->x->nY);
 
     if (out == NULL) {
+        printf("psPlaneTransformApply(): out is NULL\n");
         out = (psPlane* ) psAlloc(sizeof(psPlane));
     }
+    /*
+        out->x = psDPolynomial2DEval(
+                     transform->x,
+                     coords->x,
+                     coords->y
+                 );
+        out->y = psDPolynomial2DEval(
+                     transform->y,
+                     coords->x,
+                     coords->y
+                 );
+    */
     out->x = psDPolynomial2DEval(
                  transform->x,
@@ -251,4 +266,5 @@
                  coords->y
              );
+
     out->y = psDPolynomial2DEval(
                  transform->y,
@@ -257,4 +273,5 @@
              );
 
+    printf("psPlaneTransformApply(END)\n");
     return (out);
 }
Index: /trunk/psLib/src/astronomy/psCoord.c
===================================================================
--- /trunk/psLib/src/astronomy/psCoord.c	(revision 4285)
+++ /trunk/psLib/src/astronomy/psCoord.c	(revision 4286)
@@ -10,6 +10,6 @@
 *  @author GLG, MHPCC
 *
-*  @version $Revision: 1.71 $ $Name: not supported by cvs2svn $
-*  @date $Date: 2005-06-13 20:16:22 $
+*  @version $Revision: 1.72 $ $Name: not supported by cvs2svn $
+*  @date $Date: 2005-06-16 22:00:03 $
 *
 *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -242,8 +242,23 @@
     PS_ASSERT_PTR_NON_NULL(transform->y, NULL);
     PS_ASSERT_PTR_NON_NULL(coords, NULL);
+    printf("psPlaneTransformApply(): (%f, %f) x:(%d, %d)\n",
+           coords->x, coords->y, transform->x->nX, transform->x->nY);
 
     if (out == NULL) {
+        printf("psPlaneTransformApply(): out is NULL\n");
         out = (psPlane* ) psAlloc(sizeof(psPlane));
     }
+    /*
+        out->x = psDPolynomial2DEval(
+                     transform->x,
+                     coords->x,
+                     coords->y
+                 );
+        out->y = psDPolynomial2DEval(
+                     transform->y,
+                     coords->x,
+                     coords->y
+                 );
+    */
     out->x = psDPolynomial2DEval(
                  transform->x,
@@ -251,4 +266,5 @@
                  coords->y
              );
+
     out->y = psDPolynomial2DEval(
                  transform->y,
@@ -257,4 +273,5 @@
              );
 
+    printf("psPlaneTransformApply(END)\n");
     return (out);
 }
Index: /trunk/psLib/src/dataManip/psFunctions.c
===================================================================
--- /trunk/psLib/src/dataManip/psFunctions.c	(revision 4285)
+++ /trunk/psLib/src/dataManip/psFunctions.c	(revision 4286)
@@ -7,6 +7,6 @@
  *  polynomials.  It also contains a Gaussian functions.
  *
- *  @version $Revision: 1.106 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2005-06-13 19:39:31 $
+ *  @version $Revision: 1.107 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2005-06-16 22:00:03 $
  *
  *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -671,5 +671,7 @@
 }
 
-static psF64 dOrdPolynomial2DEval(psF64 x, psF64 y, const psDPolynomial2D* myPoly)
+static psF64 dOrdPolynomial2DEval(psF64 x,
+                                  psF64 y,
+                                  const psDPolynomial2D* myPoly)
 {
     psS32 loop_x = 0;
@@ -678,4 +680,5 @@
     psF64 xSum = 1.0;
     psF64 ySum = 1.0;
+    printf("dOrdPolynomial2DEval(%f, %f) (%d, %d)\n", x, y, myPoly->nX, myPoly->nY);
 
     for (loop_x = 0; loop_x < myPoly->nX; loop_x++) {
@@ -689,4 +692,5 @@
         xSum *= x;
     }
+    printf("dOrdPolynomial2DEval(END)\n");
 
     return(polySum);
Index: /trunk/psLib/src/image/psImageGeomManip.c
===================================================================
--- /trunk/psLib/src/image/psImageGeomManip.c	(revision 4285)
+++ /trunk/psLib/src/image/psImageGeomManip.c	(revision 4286)
@@ -10,6 +10,6 @@
  *  @author Ross Harman, MHPCC
  *
- *  @version $Revision: 1.5 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2005-06-11 02:19:05 $
+ *  @version $Revision: 1.6 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2005-06-16 22:00:04 $
  *
  *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -697,6 +697,10 @@
                           const psPixels* pixels,
                           psImageInterpolateMode mode,
-                          int exposedValue)
+                          int exposedValue,
+                          psPlaneTransform *junk)
 {
+    printf("junk(A): (%d, %d) (%d %d)\n", junk->x->nX, junk->x->nY, junk->y->nX, junk->y->nY);
+
+
     if (input == NULL) {
         psError(PS_ERR_BAD_PARAMETER_NULL, true,
@@ -707,4 +711,5 @@
     psElemType type = input->type.type;
 
+    printf("junk(B): (%d, %d) (%d %d)\n", junk->x->nX, junk->x->nY, junk->y->nX, junk->y->nY);
     if (inputMask != NULL) {
         if (input->numRows != inputMask->numRows || input->numCols != inputMask->numCols) {
@@ -733,4 +738,5 @@
     }
 
+    printf("junk(C): (%d, %d) (%d %d)\n", junk->x->nX, junk->x->nY, junk->y->nX, junk->y->nY);
     int row0;
     int row1;
@@ -774,4 +780,5 @@
     }
 
+    printf("junk(D): (%d, %d) (%d %d)\n", junk->x->nX, junk->x->nY, junk->y->nX, junk->y->nY);
     // loop through the output image using the domain above and transform
     // each output pixel to input coordinates and use psImagePixelInterpolate
@@ -780,8 +787,14 @@
     psPlane* inPosition = NULL;
 
+
+
+
+
     #define PSIMAGE_TRANSFORM_DOTRANSFORM(TYPE,MODE) \
     /* apply the transform to get the position in the input image */ \
+    printf("junk(Ea): (%d, %d) (%d %d)\n", junk->x->nX, junk->x->nY, junk->y->nX, junk->y->nY);\
     inPosition = psPlaneTransformApply(inPosition, outToIn, &outPosition); \
     \
+    printf("junk(Eaa): (%d, %d) (%d %d)\n", junk->x->nX, junk->x->nY, junk->y->nX, junk->y->nY);\
     if (inPosition == NULL) { \
         psError(PS_ERR_UNKNOWN, false, \
@@ -791,7 +804,10 @@
     } \
     /* interpolate the cooresponding input pixel to get the output pixel value. */ \
+    printf("junk(Eaaa): (%d, %d) (%d %d)\n", junk->x->nX, junk->x->nY, junk->y->nX, junk->y->nY);\
     ps##TYPE value = p_psImagePixelInterpolate##MODE##_##TYPE(input, \
                      inPosition->x, inPosition->y, \
                      inputMask, inputMaskVal, NAN); \
+    /*    psFree(inPosition); */\
+    printf("junk(Eaaaa): (%d, %d) (%d %d)\n", junk->x->nX, junk->x->nY, junk->y->nX, junk->y->nY);\
     if (isnan(value)) { \
         if (blankPixels != NULL) { \
@@ -800,4 +816,7 @@
         value = exposedValue; \
     } \
+    printf("junk(Eb): (%d, %d) (%d %d)\n", junk->x->nX, junk->x->nY, junk->y->nX, junk->y->nY);
+
+
 
     #define PSIMAGE_TRANSFORM_LOOP(TYPE, MODE) { \
@@ -806,10 +825,19 @@
             ps##TYPE* outputData=output->data.TYPE[row]; \
             for (int col = 0; col < numCols; col++) { \
+                printf("junk(Eb): (%d, %d) (%d %d)\n", junk->x->nX, junk->x->nY, junk->y->nX, junk->y->nY); \
                 outPosition.x = col+col0; \
+                printf("junk(Ebb): (%d, %d) (%d %d)\n", junk->x->nX, junk->x->nY, junk->y->nX, junk->y->nY); \
                 PSIMAGE_TRANSFORM_DOTRANSFORM(TYPE,MODE) \
+                printf("junk(Ebbb): (%d, %d) (%d %d)\n", junk->x->nX, junk->x->nY, junk->y->nX, junk->y->nY); \
                 outputData[col] = value; \
-            } \
-        } \
-    }
+                printf("junk(Ebbbb): (%d, %d) (%d %d)\n", junk->x->nX, junk->x->nY, junk->y->nX, junk->y->nY); \
+                printf("HERE B: Writing (%d, %d) Image size is (%d, %d)\n", \
+                       row, col, output->numRows, output->numCols); \
+                printf("junk(Ebbbbb): (%d, %d) (%d %d)\n", junk->x->nX, junk->x->nY, junk->y->nX, junk->y->nY); \
+            } \
+        } \
+    }
+
+    printf("junk(F): (%d, %d) (%d %d)\n", junk->x->nX, junk->x->nY, junk->y->nX, junk->y->nY);
 
     #define PSIMAGE_TRANSFORM_FROMLIST(TYPE, MODE) { \
@@ -823,7 +851,10 @@
                 PSIMAGE_TRANSFORM_DOTRANSFORM(TYPE,MODE) \
                 output->data.TYPE[y][x] = value; \
-            } \
-        } \
-    }
+                printf("HERE C: Writing (%d, %d) Image size is (%d, %d)\n", y, x, output->numRows, output->numCols); \
+            } \
+        } \
+    }
+
+    printf("junk(G): (%d, %d) (%d %d)\n", junk->x->nX, junk->x->nY, junk->y->nX, junk->y->nY);
 
     #define PSIMAGE_TRANSFORM_CASE(MODE) \
@@ -843,4 +874,5 @@
                     typeStr); \
             psFree(output); \
+            printf("junk(H): (%d, %d) (%d %d)\n", junk->x->nX, junk->x->nY, junk->y->nX, junk->y->nY); \
             return NULL; \
         } \
@@ -848,4 +880,5 @@
     break;
 
+    printf("junk(J): (%d, %d) (%d %d)\n", junk->x->nX, junk->x->nY, junk->y->nX, junk->y->nY);
     switch (mode) {
         PSIMAGE_TRANSFORM_CASE(FLAT);
@@ -857,6 +890,8 @@
                 mode);
         psFree(output);
-        return NULL;
-    }
+        printf("junk(K): (%d, %d) (%d %d)\n", junk->x->nX, junk->x->nY, junk->y->nX, junk->y->nY);
+        return NULL;
+    }
+    printf("junk(L): (%d, %d) (%d %d)\n", junk->x->nX, junk->x->nY, junk->y->nX, junk->y->nY);
 
     return output;
Index: /trunk/psLib/src/image/psImageGeomManip.h
===================================================================
--- /trunk/psLib/src/image/psImageGeomManip.h	(revision 4285)
+++ /trunk/psLib/src/image/psImageGeomManip.h	(revision 4286)
@@ -8,6 +8,6 @@
  *  @author Robert DeSonia, MHPCC
  *
- *  @version $Revision: 1.4 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2005-06-14 02:54:15 $
+ *  @version $Revision: 1.5 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2005-06-16 22:00:04 $
  *
  *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -153,5 +153,6 @@
                                              *  output image shall be transformed; otherwise, entire image generated*/
     psImageInterpolateMode mode,       ///< the interpolation scheme to be used
-    int exposedValue                   ///< Exposed value to which non-corresponding pixels are set
+    int exposedValue,                   ///< Exposed value to which non-corresponding pixels are set
+    psPlaneTransform *junk
 );
 
Index: /trunk/psLib/src/imageops/psImageGeomManip.c
===================================================================
--- /trunk/psLib/src/imageops/psImageGeomManip.c	(revision 4285)
+++ /trunk/psLib/src/imageops/psImageGeomManip.c	(revision 4286)
@@ -10,6 +10,6 @@
  *  @author Ross Harman, MHPCC
  *
- *  @version $Revision: 1.5 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2005-06-11 02:19:05 $
+ *  @version $Revision: 1.6 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2005-06-16 22:00:04 $
  *
  *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -697,6 +697,10 @@
                           const psPixels* pixels,
                           psImageInterpolateMode mode,
-                          int exposedValue)
+                          int exposedValue,
+                          psPlaneTransform *junk)
 {
+    printf("junk(A): (%d, %d) (%d %d)\n", junk->x->nX, junk->x->nY, junk->y->nX, junk->y->nY);
+
+
     if (input == NULL) {
         psError(PS_ERR_BAD_PARAMETER_NULL, true,
@@ -707,4 +711,5 @@
     psElemType type = input->type.type;
 
+    printf("junk(B): (%d, %d) (%d %d)\n", junk->x->nX, junk->x->nY, junk->y->nX, junk->y->nY);
     if (inputMask != NULL) {
         if (input->numRows != inputMask->numRows || input->numCols != inputMask->numCols) {
@@ -733,4 +738,5 @@
     }
 
+    printf("junk(C): (%d, %d) (%d %d)\n", junk->x->nX, junk->x->nY, junk->y->nX, junk->y->nY);
     int row0;
     int row1;
@@ -774,4 +780,5 @@
     }
 
+    printf("junk(D): (%d, %d) (%d %d)\n", junk->x->nX, junk->x->nY, junk->y->nX, junk->y->nY);
     // loop through the output image using the domain above and transform
     // each output pixel to input coordinates and use psImagePixelInterpolate
@@ -780,8 +787,14 @@
     psPlane* inPosition = NULL;
 
+
+
+
+
     #define PSIMAGE_TRANSFORM_DOTRANSFORM(TYPE,MODE) \
     /* apply the transform to get the position in the input image */ \
+    printf("junk(Ea): (%d, %d) (%d %d)\n", junk->x->nX, junk->x->nY, junk->y->nX, junk->y->nY);\
     inPosition = psPlaneTransformApply(inPosition, outToIn, &outPosition); \
     \
+    printf("junk(Eaa): (%d, %d) (%d %d)\n", junk->x->nX, junk->x->nY, junk->y->nX, junk->y->nY);\
     if (inPosition == NULL) { \
         psError(PS_ERR_UNKNOWN, false, \
@@ -791,7 +804,10 @@
     } \
     /* interpolate the cooresponding input pixel to get the output pixel value. */ \
+    printf("junk(Eaaa): (%d, %d) (%d %d)\n", junk->x->nX, junk->x->nY, junk->y->nX, junk->y->nY);\
     ps##TYPE value = p_psImagePixelInterpolate##MODE##_##TYPE(input, \
                      inPosition->x, inPosition->y, \
                      inputMask, inputMaskVal, NAN); \
+    /*    psFree(inPosition); */\
+    printf("junk(Eaaaa): (%d, %d) (%d %d)\n", junk->x->nX, junk->x->nY, junk->y->nX, junk->y->nY);\
     if (isnan(value)) { \
         if (blankPixels != NULL) { \
@@ -800,4 +816,7 @@
         value = exposedValue; \
     } \
+    printf("junk(Eb): (%d, %d) (%d %d)\n", junk->x->nX, junk->x->nY, junk->y->nX, junk->y->nY);
+
+
 
     #define PSIMAGE_TRANSFORM_LOOP(TYPE, MODE) { \
@@ -806,10 +825,19 @@
             ps##TYPE* outputData=output->data.TYPE[row]; \
             for (int col = 0; col < numCols; col++) { \
+                printf("junk(Eb): (%d, %d) (%d %d)\n", junk->x->nX, junk->x->nY, junk->y->nX, junk->y->nY); \
                 outPosition.x = col+col0; \
+                printf("junk(Ebb): (%d, %d) (%d %d)\n", junk->x->nX, junk->x->nY, junk->y->nX, junk->y->nY); \
                 PSIMAGE_TRANSFORM_DOTRANSFORM(TYPE,MODE) \
+                printf("junk(Ebbb): (%d, %d) (%d %d)\n", junk->x->nX, junk->x->nY, junk->y->nX, junk->y->nY); \
                 outputData[col] = value; \
-            } \
-        } \
-    }
+                printf("junk(Ebbbb): (%d, %d) (%d %d)\n", junk->x->nX, junk->x->nY, junk->y->nX, junk->y->nY); \
+                printf("HERE B: Writing (%d, %d) Image size is (%d, %d)\n", \
+                       row, col, output->numRows, output->numCols); \
+                printf("junk(Ebbbbb): (%d, %d) (%d %d)\n", junk->x->nX, junk->x->nY, junk->y->nX, junk->y->nY); \
+            } \
+        } \
+    }
+
+    printf("junk(F): (%d, %d) (%d %d)\n", junk->x->nX, junk->x->nY, junk->y->nX, junk->y->nY);
 
     #define PSIMAGE_TRANSFORM_FROMLIST(TYPE, MODE) { \
@@ -823,7 +851,10 @@
                 PSIMAGE_TRANSFORM_DOTRANSFORM(TYPE,MODE) \
                 output->data.TYPE[y][x] = value; \
-            } \
-        } \
-    }
+                printf("HERE C: Writing (%d, %d) Image size is (%d, %d)\n", y, x, output->numRows, output->numCols); \
+            } \
+        } \
+    }
+
+    printf("junk(G): (%d, %d) (%d %d)\n", junk->x->nX, junk->x->nY, junk->y->nX, junk->y->nY);
 
     #define PSIMAGE_TRANSFORM_CASE(MODE) \
@@ -843,4 +874,5 @@
                     typeStr); \
             psFree(output); \
+            printf("junk(H): (%d, %d) (%d %d)\n", junk->x->nX, junk->x->nY, junk->y->nX, junk->y->nY); \
             return NULL; \
         } \
@@ -848,4 +880,5 @@
     break;
 
+    printf("junk(J): (%d, %d) (%d %d)\n", junk->x->nX, junk->x->nY, junk->y->nX, junk->y->nY);
     switch (mode) {
         PSIMAGE_TRANSFORM_CASE(FLAT);
@@ -857,6 +890,8 @@
                 mode);
         psFree(output);
-        return NULL;
-    }
+        printf("junk(K): (%d, %d) (%d %d)\n", junk->x->nX, junk->x->nY, junk->y->nX, junk->y->nY);
+        return NULL;
+    }
+    printf("junk(L): (%d, %d) (%d %d)\n", junk->x->nX, junk->x->nY, junk->y->nX, junk->y->nY);
 
     return output;
Index: /trunk/psLib/src/imageops/psImageGeomManip.h
===================================================================
--- /trunk/psLib/src/imageops/psImageGeomManip.h	(revision 4285)
+++ /trunk/psLib/src/imageops/psImageGeomManip.h	(revision 4286)
@@ -8,6 +8,6 @@
  *  @author Robert DeSonia, MHPCC
  *
- *  @version $Revision: 1.4 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2005-06-14 02:54:15 $
+ *  @version $Revision: 1.5 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2005-06-16 22:00:04 $
  *
  *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -153,5 +153,6 @@
                                              *  output image shall be transformed; otherwise, entire image generated*/
     psImageInterpolateMode mode,       ///< the interpolation scheme to be used
-    int exposedValue                   ///< Exposed value to which non-corresponding pixels are set
+    int exposedValue,                   ///< Exposed value to which non-corresponding pixels are set
+    psPlaneTransform *junk
 );
 
Index: /trunk/psLib/src/math/psPolynomial.c
===================================================================
--- /trunk/psLib/src/math/psPolynomial.c	(revision 4285)
+++ /trunk/psLib/src/math/psPolynomial.c	(revision 4286)
@@ -7,6 +7,6 @@
  *  polynomials.  It also contains a Gaussian functions.
  *
- *  @version $Revision: 1.106 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2005-06-13 19:39:31 $
+ *  @version $Revision: 1.107 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2005-06-16 22:00:03 $
  *
  *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -671,5 +671,7 @@
 }
 
-static psF64 dOrdPolynomial2DEval(psF64 x, psF64 y, const psDPolynomial2D* myPoly)
+static psF64 dOrdPolynomial2DEval(psF64 x,
+                                  psF64 y,
+                                  const psDPolynomial2D* myPoly)
 {
     psS32 loop_x = 0;
@@ -678,4 +680,5 @@
     psF64 xSum = 1.0;
     psF64 ySum = 1.0;
+    printf("dOrdPolynomial2DEval(%f, %f) (%d, %d)\n", x, y, myPoly->nX, myPoly->nY);
 
     for (loop_x = 0; loop_x < myPoly->nX; loop_x++) {
@@ -689,4 +692,5 @@
         xSum *= x;
     }
+    printf("dOrdPolynomial2DEval(END)\n");
 
     return(polySum);
Index: /trunk/psLib/src/math/psSpline.c
===================================================================
--- /trunk/psLib/src/math/psSpline.c	(revision 4285)
+++ /trunk/psLib/src/math/psSpline.c	(revision 4286)
@@ -7,6 +7,6 @@
  *  polynomials.  It also contains a Gaussian functions.
  *
- *  @version $Revision: 1.106 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2005-06-13 19:39:31 $
+ *  @version $Revision: 1.107 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2005-06-16 22:00:03 $
  *
  *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -671,5 +671,7 @@
 }
 
-static psF64 dOrdPolynomial2DEval(psF64 x, psF64 y, const psDPolynomial2D* myPoly)
+static psF64 dOrdPolynomial2DEval(psF64 x,
+                                  psF64 y,
+                                  const psDPolynomial2D* myPoly)
 {
     psS32 loop_x = 0;
@@ -678,4 +680,5 @@
     psF64 xSum = 1.0;
     psF64 ySum = 1.0;
+    printf("dOrdPolynomial2DEval(%f, %f) (%d, %d)\n", x, y, myPoly->nX, myPoly->nY);
 
     for (loop_x = 0; loop_x < myPoly->nX; loop_x++) {
@@ -689,4 +692,5 @@
         xSum *= x;
     }
+    printf("dOrdPolynomial2DEval(END)\n");
 
     return(polySum);