Changeset 39018 for trunk/Ohana/src/addstar/src/loadgalphot_fit2d.c

Timestamp:

Nov 2, 2015, 7:49:22 AM (11 years ago)

Author:

eugene

Message:

add -forcedgalaxy mode to mkcmf; fix fitting error in loadgalphot

File:

• trunk/Ohana/src/addstar/src/loadgalphot_fit2d.c (modified) (3 diffs)

trunk/Ohana/src/addstar/src/loadgalphot_fit2d.c

@@ -183 +183 @@
       if (!finite(*x) || !finite(*y) || !finite(*z)) continue;
       *zf = 0;
-      double X = 1;
+      double Y = 1;
       for (i = ny = 0; ny < fit->nterm; ny++) {
-        double Y = X;
+        double X = Y;
         for (nx = 0; nx < fit->nterm - ny; nx++, i++) {
-          *zf += fit->b[i][0]*Y;
-          Y = Y * (*y);
-        }
-        X = X * (*x);
+          *zf += fit->b[i][0]*X;
+          X = X * (*x);
+        }
+        Y = Y * (*y);
       }
     }
@@ -213 +213 @@
     maxsigma = fit->ClipNsigma * sigma;
 
-    // if (VERBOSE) gprint (GP_ERR, "mean: %g, sigma: %g, maxsigma: %g\n", mean, sigma, maxsigma);
+    if (VERBOSE) gprint (GP_ERR, "mean: %g, sigma: %g, maxsigma: %g\n", mean, sigma, maxsigma);
 
     /* mask outlier points */
@@ -230 +230 @@
       } 
     }
-    // if (VERBOSE) gprint (GP_ERR, "pass: %d, Nmask: %d\n", N, Nmask);
+    if (VERBOSE) {
+      gprint (GP_ERR, "pass: %d, Nmask: %d\n", N, Nmask);
+    }
   }