Index: trunk/psLib/src/imageops/psImageStats.c
===================================================================
--- trunk/psLib/src/imageops/psImageStats.c	(revision 1448)
+++ trunk/psLib/src/imageops/psImageStats.c	(revision 1719)
@@ -10,6 +10,6 @@
 *  @author George Gusciora, MHPCC
 *
-*  @version $Revision: 1.33 $ $Name: not supported by cvs2svn $
-*  @date $Date: 2004-08-10 01:55:34 $
+*  @version $Revision: 1.34 $ $Name: not supported by cvs2svn $
+*  @date $Date: 2004-09-08 06:02:47 $
 *
 *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -185,11 +185,11 @@
 psPolynomial1D* *p_psCreateChebyshevPolys(int maxChebyPoly)
 {
-    psPolynomial1D* *chebPolys = NULL;
+    psPolynomial1D **chebPolys = NULL;
     int i = 0;
     int j = 0;
 
-    chebPolys = (psPolynomial1D* *) psAlloc(maxChebyPoly * sizeof(psPolynomial1D* ));
+    chebPolys = (psPolynomial1D **) psAlloc(maxChebyPoly * sizeof(psPolynomial1D *));
     for (i = 0; i < maxChebyPoly; i++) {
-        chebPolys[i] = psPolynomial1DAlloc(i + 1);
+        chebPolys[i] = psPolynomial1DAlloc(i + 1, PS_POLYNOMIAL_ORD);
     }
 
@@ -333,5 +333,5 @@
     int i = 0;
     int j = 0;
-    float **sums = NULL;
+    //    float **sums = NULL;
     psPolynomial1D* *chebPolys = NULL;
     int maxChebyPoly = 0;
@@ -340,18 +340,17 @@
     float polySum = 0.0;
 
-    // Create the sums[][] data structure.  This
-    // will hold the LHS of
-    // equation
-    // 29 in the ADD: sums[k][l] = SUM {
-    // image(x,y) * Tk(x) * Tl(y) }
-    sums = (float **)psAlloc(coeffs->nX * sizeof(float *));
-    for (i = 0; i < coeffs->nX; i++) {
-        sums[i] = (float *)psAlloc(coeffs->nY * sizeof(float));
-    }
-    for (i = 0; i < coeffs->nX; i++) {
-        for (j = 0; j < coeffs->nY; j++) {
-            sums[i][j] = 0.0;
-        }
-    }
+    // Create the sums[][] data structure.  This will hold the LHS of
+    // equation 29 in the ADD:
+    //         sums[k][l] = SUM { image(x,y) * Tk(x) * Tl(y) }
+    //
+    //    sums = (float **)psAlloc(coeffs->nX * sizeof(float *));
+    //    for (i = 0; i < coeffs->nX; i++) {
+    //        sums[i] = (float *)psAlloc(coeffs->nY * sizeof(float));
+    //    }
+    //    for (i = 0; i < coeffs->nX; i++) {
+    //        for (j = 0; j < coeffs->nY; j++) {
+    //            sums[i][j] = 0.0;
+    //        }
+    //    }
 
     // We scale the pixel positions to values
@@ -375,8 +374,7 @@
                 for (j = 0; j < coeffs->nY; j++) {
                     polySum +=
-                        psPolynomial1DEval
-                        (rScalingFactors[x],
-                         chebPolys[i]) *
-                        psPolynomial1DEval(cScalingFactors[y], chebPolys[j]) * coeffs->coeff[i][j];
+                        psPolynomial1DEval(rScalingFactors[x], chebPolys[i]) *
+                        psPolynomial1DEval(cScalingFactors[y], chebPolys[j]) *
+                        coeffs->coeff[i][j];
 
                 }
@@ -394,8 +392,8 @@
 
     // Free some data
-    for (i = 0; i < coeffs->nX; i++) {
-        psFree(sums[i]);
-    }
-    psFree(sums);
+    //    for (i = 0; i < coeffs->nX; i++) {
+    //        psFree(sums[i]);
+    //    }
+    //    psFree(sums);
     psFree(cScalingFactors);
     psFree(rScalingFactors);