Index: trunk/psLib/src/imageops/psImageStats.c
===================================================================
--- trunk/psLib/src/imageops/psImageStats.c	(revision 5213)
+++ trunk/psLib/src/imageops/psImageStats.c	(revision 5294)
@@ -9,6 +9,6 @@
  *  @author GLG, MHPCC
  *
- *  @version $Revision: 1.82 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2005-09-30 23:09:13 $
+ *  @version $Revision: 1.83 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2005-10-12 21:02:20 $
  *
  *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -233,8 +233,8 @@
     chebPolys[1]->coeff[1] = 1;
     for (i = 2; i < maxChebyPoly; i++) {
-        for (j = 0; j < (1 + chebPolys[i - 1]->COOL_1D_n); j++) {
+        for (j = 0; j < (1 + chebPolys[i - 1]->nX); j++) {
             chebPolys[i]->coeff[j + 1] = 2 * chebPolys[i - 1]->coeff[j];
         }
-        for (j = 0; j < (1 + chebPolys[i - 2]->COOL_1D_n); j++) {
+        for (j = 0; j < (1 + chebPolys[i - 2]->nX); j++) {
             chebPolys[i]->coeff[j] -= chebPolys[i - 2]->coeff[j];
         }
@@ -291,7 +291,7 @@
     // 29 in the ADD: sums[k][l] = SUM {
     // image(x,y) * Tk(x) * Tl(y) }
-    sums = (double **)psAlloc((1 + coeffs->COOL_2D_nX) * sizeof(double *));
-    for (i = 0; i < (1 + coeffs->COOL_2D_nX); i++) {
-        sums[i] = (double *)psAlloc((1 + coeffs->COOL_2D_nY) * sizeof(double));
+    sums = (double **)psAlloc((1 + coeffs->nX) * sizeof(double *));
+    for (i = 0; i < (1 + coeffs->nX); i++) {
+        sums[i] = (double *)psAlloc((1 + coeffs->nY) * sizeof(double));
     }
     // We scale the pixel positions to values
@@ -303,13 +303,13 @@
     // are needed, then create them.
     // XXX: recorde or verify the poly order/nterm change
-    maxChebyPoly = coeffs->COOL_2D_nX;
-    if (coeffs->COOL_2D_nY > coeffs->COOL_2D_nX) {
-        maxChebyPoly = coeffs->COOL_2D_nY;
+    maxChebyPoly = coeffs->nX;
+    if (coeffs->nY > coeffs->nX) {
+        maxChebyPoly = coeffs->nY;
     }
     chebPolys = p_psCreateChebyshevPolys(maxChebyPoly + 1);
 
     // Compute the sums[][] data structure.
-    for (i = 0; i < (1 + coeffs->COOL_2D_nX); i++) {
-        for (j = 0; j < (1 + coeffs->COOL_2D_nY); j++) {
+    for (i = 0; i < (1 + coeffs->nX); i++) {
+        for (j = 0; j < (1 + coeffs->nY); j++) {
             sums[i][j] = 0.0;
             for (x = 0; x < input->numRows; x++) {
@@ -332,6 +332,6 @@
     }
 
-    for (i = 0; i < (1 + coeffs->COOL_2D_nX); i++) {
-        for (j = 0; j < (1 + coeffs->COOL_2D_nY); j++) {
+    for (i = 0; i < (1 + coeffs->nX); i++) {
+        for (j = 0; j < (1 + coeffs->nY); j++) {
             coeffs->coeff[i][j] = sums[i][j];
             coeffs->coeff[i][j] /= (double)(input->numRows * input->numCols);
@@ -355,5 +355,5 @@
 
     // Free some data
-    for (i = 0; i < (1 + coeffs->COOL_2D_nX); i++) {
+    for (i = 0; i < (1 + coeffs->nX); i++) {
         psFree(sums[i]);
     }
@@ -438,7 +438,7 @@
     // 29 in the ADD: sums[k][l] = SUM {
     // image(x,y) * Tk(x) * Tl(y) }
-    sums = (double **)psAlloc((1 + coeffs->COOL_2D_nX) * sizeof(double *));
-    for (i = 0; i < (1 + coeffs->COOL_2D_nX); i++) {
-        sums[i] = (double *)psAlloc((1 + coeffs->COOL_2D_nY) * sizeof(double));
+    sums = (double **)psAlloc((1 + coeffs->nX) * sizeof(double *));
+    for (i = 0; i < (1 + coeffs->nX); i++) {
+        sums[i] = (double *)psAlloc((1 + coeffs->nY) * sizeof(double));
     }
     // We scale the pixel positions to values
@@ -449,13 +449,13 @@
     // Determine how many Chebyshev polynomials
     // are needed, then create them.
-    maxChebyPoly = coeffs->COOL_2D_nX;
-    if (coeffs->COOL_2D_nY > coeffs->COOL_2D_nX) {
-        maxChebyPoly = coeffs->COOL_2D_nY;
+    maxChebyPoly = coeffs->nX;
+    if (coeffs->nY > coeffs->nX) {
+        maxChebyPoly = coeffs->nY;
     }
     chebPolys = p_psCreateChebyshevPolys(maxChebyPoly + 1);
 
     // Compute the sums[][] data structure.
-    for (i = 0; i < (1 + coeffs->COOL_2D_nX); i++) {
-        for (j = 0; j < (1 + coeffs->COOL_2D_nY); j++) {
+    for (i = 0; i < (1 + coeffs->nX); i++) {
+        for (j = 0; j < (1 + coeffs->nY); j++) {
             sums[i][j] = 0.0;
             for (x = 0; x < input->numRows; x++) {
@@ -481,6 +481,6 @@
     }
 
-    for (i = 0; i < (1 + coeffs->COOL_2D_nX); i++) {
-        for (j = 0; j < (1 + coeffs->COOL_2D_nY); j++) {
+    for (i = 0; i < (1 + coeffs->nX); i++) {
+        for (j = 0; j < (1 + coeffs->nY); j++) {
             coeffs->coeff[i][j] = sums[i][j];
             coeffs->coeff[i][j] /= (double)(input->numRows * input->numCols);
@@ -504,5 +504,5 @@
 
     // Free some data
-    for (i = 0; i < (1 + coeffs->COOL_2D_nX); i++) {
+    for (i = 0; i < (1 + coeffs->nX); i++) {
         psFree(sums[i]);
     }
@@ -540,7 +540,7 @@
     // Determine how many Chebyshev polynomials
     // are needed, then create them.
-    maxChebyPoly = coeffs->COOL_2D_nX;
-    if (coeffs->COOL_2D_nY > coeffs->COOL_2D_nX) {
-        maxChebyPoly = coeffs->COOL_2D_nY;
+    maxChebyPoly = coeffs->nX;
+    if (coeffs->nY > coeffs->nX) {
+        maxChebyPoly = coeffs->nY;
     }
 
@@ -550,6 +550,6 @@
         for (y = 0; y < input->numCols; y++) {
             polySum = 0.0;
-            for (i = 0; i < (1 + coeffs->COOL_2D_nX); i++) {
-                for (j = 0; j < (1 + coeffs->COOL_2D_nY); j++) {
+            for (i = 0; i < (1 + coeffs->nX); i++) {
+                for (j = 0; j < (1 + coeffs->nY); j++) {
                     polySum +=
                         psPolynomial1DEval(chebPolys[i], rScalingFactors[x]) *