Changeset 5294 for trunk/psModules/src/imsubtract

Timestamp:

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

Author:

gusciora

Message:

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

Location:

trunk/psModules/src/imsubtract

Files:

• pmSubtractBias.c (modified) (2 diffs)
• pmSubtractSky.c (modified) (6 diffs)

trunk/psModules/src/imsubtract/pmSubtractBias.c

@@ -6 +6 @@
  *  @author GLG, MHPCC
  *
- *  @version $Revision: 1.2 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2005-09-29 21:57:30 $
+ *  @version $Revision: 1.3 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2005-10-12 21:02:04 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -26 +26 @@
 {
     printf("-------------- PS_POLY1D_PRINT() --------------\n");
-    printf("poly->COOL_1D_n is %d\n", poly->COOL_1D_n);
-    for (psS32 i = 0 ; i < (1 + poly->COOL_1D_n) ; i++) {
+    printf("poly->nX is %d\n", poly->nX);
+    for (psS32 i = 0 ; i < (1 + poly->nX) ; i++) {
         printf("poly->coeff[%d] is %f\n", i, poly->coeff[i]);
     }

trunk/psModules/src/imsubtract/pmSubtractSky.c

@@ -6 +6 @@
  *  @author GLG, MHPCC
  *
- *  @version $Revision: 1.1 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2005-09-28 20:43:52 $
+ *  @version $Revision: 1.2 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2005-10-12 21:02:04 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -306 +306 @@
     // polynomial if there are not enough rows/columns in the input image.
 
-    if ((myPoly->COOL_2D_nX + 1) > dataImage->numRows) {
+    if ((myPoly->nX + 1) > dataImage->numRows) {
         psLogMsg(__func__, PS_LOG_WARN,
                  "WARNING: ImageFitPolynomial(): Reducing polynomial complexity in x-dimension.\n");
-        oldPolyX = myPoly->COOL_2D_nX;
-        myPoly->COOL_2D_nX = dataImage->numRows - 1;
-    }
-    if ((myPoly->COOL_2D_nY + 1) > dataImage->numCols) {
+        oldPolyX = myPoly->nX;
+        myPoly->nX = dataImage->numRows - 1;
+    }
+    if ((myPoly->nY + 1) > dataImage->numCols) {
         psLogMsg(__func__, PS_LOG_WARN,
                  "WARNING: ImageFitPolynomial(): Reducing polynomial complexity in y-dimension.\n");
-        oldPolyY = myPoly->COOL_2D_nY;
-        myPoly->COOL_2D_nY = dataImage->numCols - 1;
+        oldPolyY = myPoly->nY;
+        myPoly->nY = dataImage->numCols - 1;
     }
     psS32 i;
@@ -324 +324 @@
     psS32 aRow;
     psS32 aCol;
-    psS32 **polyTerms = buildPolyTerms(myPoly->COOL_2D_nX, myPoly->COOL_2D_nY);
+    psS32 **polyTerms = buildPolyTerms(myPoly->nX, myPoly->nY);
     // We determine how many coefficients will be in the polynomial that we
     // are fitting to this image.
-    psS32 localPolyTerms = CalculatePolyTerms(myPoly->COOL_2D_nX, myPoly->COOL_2D_nY);
+    psS32 localPolyTerms = CalculatePolyTerms(myPoly->nX, myPoly->nY);
     psImage *A = psImageAlloc(localPolyTerms, localPolyTerms, PS_TYPE_F64);
     psImage *Aout = psImageAlloc(localPolyTerms, localPolyTerms, PS_TYPE_F64);
@@ -345 +345 @@
         for (y=0;y<dataImage->numCols;y++) {
             if (maskImage->data.U8[x][y] == 0) {
-                buildSums((psF64) x, (psF64) y, myPoly->COOL_2D_nX, myPoly->COOL_2D_nY);
+                buildSums((psF64) x, (psF64) y, myPoly->nX, myPoly->nY);
 
                 /************************************************************
@@ -428 +428 @@
     //
     if (oldPolyX != -1) {
-        myPoly->COOL_2D_nX = oldPolyX;
-        for (i=oldPolyX ; i < (1 + myPoly->COOL_2D_nX) ; i++) {
-            for (j=0;j<(1 + myPoly->COOL_2D_nY) ; j++) {
+        myPoly->nX = oldPolyX;
+        for (i=oldPolyX ; i < (1 + myPoly->nX) ; i++) {
+            for (j=0;j<(1 + myPoly->nY) ; j++) {
                 myPoly->coeff[i][j] = 0.0;
             }
@@ -436 +436 @@
     }
     if (oldPolyY != -1) {
-        myPoly->COOL_2D_nY = oldPolyY;
-        for (i=0 ; i < (1 + myPoly->COOL_2D_nX) ; i++) {
-            for (j=oldPolyY;j < (1 + myPoly->COOL_2D_nY) ; j++) {
+        myPoly->nY = oldPolyY;
+        for (i=0 ; i < (1 + myPoly->nX) ; i++) {
+            for (j=oldPolyY;j < (1 + myPoly->nY) ; j++) {
                 myPoly->coeff[i][j] = 0.0;
             }