Changeset 5294 for trunk/psLib/src/imageops/psImageStats.c

Timestamp:

Oct 12, 2005, 11:02:20 AM (21 years ago)

Author:

gusciora

Message:

Polynomial members are now nX, nY, nZ, nT.

File:

• trunk/psLib/src/imageops/psImageStats.c (modified) (12 diffs)

trunk/psLib/src/imageops/psImageStats.c

@@ -9 +9 @@
  *  @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 +233 @@
     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 +291 @@
     // 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 +303 @@
     // 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 +332 @@
     }
 
-    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 +355 @@
 
     // 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 +438 @@
     // 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 +449 @@
     // 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 +481 @@
     }
 
-    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 +504 @@
 
     // 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 +540 @@
     // 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 +550 @@
         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]) *