Index: /trunk/psLib/src/astro/psCoord.c
===================================================================
--- /trunk/psLib/src/astro/psCoord.c	(revision 4224)
+++ /trunk/psLib/src/astro/psCoord.c	(revision 4225)
@@ -10,6 +10,6 @@
 *  @author GLG, MHPCC
 *
-*  @version $Revision: 1.70 $ $Name: not supported by cvs2svn $
-*  @date $Date: 2005-05-25 20:26:55 $
+*  @version $Revision: 1.71 $ $Name: not supported by cvs2svn $
+*  @date $Date: 2005-06-13 20:16:22 $
 *
 *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -256,4 +256,5 @@
                  coords->y
              );
+
     return (out);
 }
@@ -852,6 +853,6 @@
  *****************************************************************************/
 
-psDPolynomial2D *multiplyDPoly2D(psDPolynomial2D *trans1,
-                                 psDPolynomial2D *trans2)
+static psDPolynomial2D *multiplyDPoly2D(psDPolynomial2D *trans1,
+                                        psDPolynomial2D *trans2)
 {
     //TRACE: printf("multiplyDPoly2D(%d %d: %d %d)\n", trans1->nX, trans1->nY, trans2->nX, trans2->nY);
@@ -1040,4 +1041,5 @@
 XXX: What about nRejIter?  Iterations?
 XXX: Use static vectors for internal data.
+XXX: This code has problems with data that corresponds to a non-linear fit.
  *****************************************************************************/
 bool psPlaneTransformFit(psPlaneTransform *trans,
Index: /trunk/psLib/src/astronomy/psCoord.c
===================================================================
--- /trunk/psLib/src/astronomy/psCoord.c	(revision 4224)
+++ /trunk/psLib/src/astronomy/psCoord.c	(revision 4225)
@@ -10,6 +10,6 @@
 *  @author GLG, MHPCC
 *
-*  @version $Revision: 1.70 $ $Name: not supported by cvs2svn $
-*  @date $Date: 2005-05-25 20:26:55 $
+*  @version $Revision: 1.71 $ $Name: not supported by cvs2svn $
+*  @date $Date: 2005-06-13 20:16:22 $
 *
 *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -256,4 +256,5 @@
                  coords->y
              );
+
     return (out);
 }
@@ -852,6 +853,6 @@
  *****************************************************************************/
 
-psDPolynomial2D *multiplyDPoly2D(psDPolynomial2D *trans1,
-                                 psDPolynomial2D *trans2)
+static psDPolynomial2D *multiplyDPoly2D(psDPolynomial2D *trans1,
+                                        psDPolynomial2D *trans2)
 {
     //TRACE: printf("multiplyDPoly2D(%d %d: %d %d)\n", trans1->nX, trans1->nY, trans2->nX, trans2->nY);
@@ -1040,4 +1041,5 @@
 XXX: What about nRejIter?  Iterations?
 XXX: Use static vectors for internal data.
+XXX: This code has problems with data that corresponds to a non-linear fit.
  *****************************************************************************/
 bool psPlaneTransformFit(psPlaneTransform *trans,
Index: /trunk/psLib/src/collections/psPixels.c
===================================================================
--- /trunk/psLib/src/collections/psPixels.c	(revision 4224)
+++ /trunk/psLib/src/collections/psPixels.c	(revision 4225)
@@ -7,6 +7,6 @@
  *  @author Robert DeSonia, MHPCC
  *
- *  @version $Revision: 1.7 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2005-06-11 02:19:05 $
+ *  @version $Revision: 1.8 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2005-06-13 20:18:18 $
  *
  *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -86,4 +86,6 @@
 
     pixels->nalloc = size;
+    pixels->n = size;
+
     if (pixels->n > pixels->nalloc) {
         pixels->n = pixels->nalloc;
Index: /trunk/psLib/src/dataManip/psMinimize.c
===================================================================
--- /trunk/psLib/src/dataManip/psMinimize.c	(revision 4224)
+++ /trunk/psLib/src/dataManip/psMinimize.c	(revision 4225)
@@ -9,6 +9,6 @@
  *  @author GLG, MHPCC
  *
- *  @version $Revision: 1.121 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2005-06-07 02:29:44 $
+ *  @version $Revision: 1.122 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2005-06-13 20:18:18 $
  *
  *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -688,7 +688,7 @@
     for (int j = 0; j < params->n; j++) {
         for (int k = 0; k < params->n; k++) {
-            alpha->data.F64[j][k] = 0;
-        }
-        beta->data.F64[j] = 0;
+            alpha->data.F64[j][k] = 0.0;
+        }
+        beta->data.F64[j] = 0.0;
     }
     chisq = 0.0;
@@ -780,5 +780,5 @@
 }
 
-// XXX" put this in constants.h
+// XXX: put this in constants.h
 # define PS_SWAP(X,Y) {double tmp=(X); (X) = (Y); (Y) = tmp;}
 
Index: /trunk/psLib/src/dataManip/psStats.c
===================================================================
--- /trunk/psLib/src/dataManip/psStats.c	(revision 4224)
+++ /trunk/psLib/src/dataManip/psStats.c	(revision 4225)
@@ -14,6 +14,6 @@
  *      stats->binsize
  *
- *  @version $Revision: 1.132 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2005-06-08 23:40:45 $
+ *  @version $Revision: 1.133 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2005-06-13 20:18:18 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -77,5 +77,4 @@
 psBool p_psGetStatValue(const psStats* stats, psF64 *value)
 {
-
     switch (stats->options & ~(PS_STAT_USE_RANGE | PS_STAT_USE_BINSIZE | PS_STAT_ROBUST_FOR_SAMPLE)) {
     case PS_STAT_SAMPLE_MEAN:
Index: /trunk/psLib/src/math/psMinimize.c
===================================================================
--- /trunk/psLib/src/math/psMinimize.c	(revision 4224)
+++ /trunk/psLib/src/math/psMinimize.c	(revision 4225)
@@ -9,6 +9,6 @@
  *  @author GLG, MHPCC
  *
- *  @version $Revision: 1.121 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2005-06-07 02:29:44 $
+ *  @version $Revision: 1.122 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2005-06-13 20:18:18 $
  *
  *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -688,7 +688,7 @@
     for (int j = 0; j < params->n; j++) {
         for (int k = 0; k < params->n; k++) {
-            alpha->data.F64[j][k] = 0;
-        }
-        beta->data.F64[j] = 0;
+            alpha->data.F64[j][k] = 0.0;
+        }
+        beta->data.F64[j] = 0.0;
     }
     chisq = 0.0;
@@ -780,5 +780,5 @@
 }
 
-// XXX" put this in constants.h
+// XXX: put this in constants.h
 # define PS_SWAP(X,Y) {double tmp=(X); (X) = (Y); (Y) = tmp;}
 
Index: /trunk/psLib/src/math/psStats.c
===================================================================
--- /trunk/psLib/src/math/psStats.c	(revision 4224)
+++ /trunk/psLib/src/math/psStats.c	(revision 4225)
@@ -14,6 +14,6 @@
  *      stats->binsize
  *
- *  @version $Revision: 1.132 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2005-06-08 23:40:45 $
+ *  @version $Revision: 1.133 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2005-06-13 20:18:18 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -77,5 +77,4 @@
 psBool p_psGetStatValue(const psStats* stats, psF64 *value)
 {
-
     switch (stats->options & ~(PS_STAT_USE_RANGE | PS_STAT_USE_BINSIZE | PS_STAT_ROBUST_FOR_SAMPLE)) {
     case PS_STAT_SAMPLE_MEAN:
Index: /trunk/psLib/src/types/psPixels.c
===================================================================
--- /trunk/psLib/src/types/psPixels.c	(revision 4224)
+++ /trunk/psLib/src/types/psPixels.c	(revision 4225)
@@ -7,6 +7,6 @@
  *  @author Robert DeSonia, MHPCC
  *
- *  @version $Revision: 1.7 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2005-06-11 02:19:05 $
+ *  @version $Revision: 1.8 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2005-06-13 20:18:18 $
  *
  *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -86,4 +86,6 @@
 
     pixels->nalloc = size;
+    pixels->n = size;
+
     if (pixels->n > pixels->nalloc) {
         pixels->n = pixels->nalloc;
