Index: trunk/psLib/src/image/psImageStats.c
===================================================================
--- trunk/psLib/src/image/psImageStats.c	(revision 1898)
+++ trunk/psLib/src/image/psImageStats.c	(revision 1939)
@@ -10,6 +10,6 @@
 *  @author George Gusciora, MHPCC
 *
-*  @version $Revision: 1.37 $ $Name: not supported by cvs2svn $
-*  @date $Date: 2004-09-25 19:11:01 $
+*  @version $Revision: 1.38 $ $Name: not supported by cvs2svn $
+*  @date $Date: 2004-09-30 23:06:36 $
 *
 *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -195,5 +195,5 @@
 }
 
-psPolynomial1D* *p_psCreateChebyshevPolys(int maxChebyPoly)
+psPolynomial1D **p_psCreateChebyshevPolys(int maxChebyPoly)
 {
     psPolynomial1D **chebPolys = NULL;
@@ -247,5 +247,4 @@
     float *cScalingFactors = NULL;
     float *rScalingFactors = NULL;
-    float tmp = 0.0;
 
     // Create the sums[][] data structure.  This
@@ -258,5 +257,4 @@
         sums[i] = (float *)psAlloc(coeffs->nY * sizeof(float));
     }
-
     // We scale the pixel positions to values
     // between -1.0 and 1.0
@@ -272,18 +270,21 @@
     chebPolys = p_psCreateChebyshevPolys(maxChebyPoly);
 
-    // Sanity check for the Chebyshevs.
-    for (i = 0; i < coeffs->nX; i++) {
-        for (j = 0; j < coeffs->nY; j++) {
-            tmp = 0.0;
-            for (x = 0; x < input->numRows; x++) {
-                tmp +=
-                    psPolynomial1DEval
-                    (rScalingFactors[x], chebPolys[i]) * psPolynomial1DEval(rScalingFactors[x], chebPolys[j]);
-
+
+    /*
+        // Sanity check for the Chebyshevs.
+        for (i = 0; i < coeffs->nX; i++) {
+            for (j = 0; j < coeffs->nY; j++) {
+                tmp = 0.0;
+                for (x = 0; x < input->numRows; x++) {
+                    tmp +=
+                        psPolynomial1DEval
+                        (rScalingFactors[x], chebPolys[i]) * psPolynomial1DEval(rScalingFactors[x], chebPolys[j]);
+     
+                }
+                // printf("SUM(Cheby(%d) * Cheby(%d))
+                // is %f\n", i, j, tmp);
             }
-            // printf("SUM(Cheby(%d) * Cheby(%d))
-            // is %f\n", i, j, tmp);
-        }
-    }
+        }
+    */
 
     // Compute the sums[][] data structure.
@@ -295,7 +296,6 @@
                     sums[i][j] +=
                         input->data.F32[x][y] *
-                        psPolynomial1DEval
-                        (rScalingFactors[x],
-                         chebPolys[i]) * psPolynomial1DEval(cScalingFactors[y], chebPolys[j]);
+                        psPolynomial1DEval(rScalingFactors[x], chebPolys[i]) *
+                        psPolynomial1DEval(cScalingFactors[y], chebPolys[j]);
                 }
             }