Changeset 25638 for branches/eam_branches/20090715/psModules/src/objects/pmSourcePhotometry.c

Timestamp:

Sep 28, 2009, 5:51:13 PM (17 years ago)

Author:

eugene

Message:

remove old test code

File:

• branches/eam_branches/20090715/psModules/src/objects/pmSourcePhotometry.c (modified) (1 diff)

branches/eam_branches/20090715/psModules/src/objects/pmSourcePhotometry.c

@@ -750 +750 @@
     return flux;
 }
-
-// XXX this is test code to verify the shift is doing the right thing (seems to be)
-# if (0)
-// measure centroid of unshifted gaussian (should be 16.0,16.0)
-        {
-          psImage *image = source->pixels;
-          float xo = 0.0;
-          float yo = 0.0;
-          float xo2 = 0.0;
-          float yo2 = 0.0;
-          float no = 0.0;
-          for (int j = 0; j < image->numRows; j++)
-          {
-            for (int i = 0; i < image->numCols; i++) {
-              xo += i*image->data.F32[j][i];
-              yo += j*image->data.F32[j][i];
-              xo2 += i*i*image->data.F32[j][i];
-              yo2 += j*j*image->data.F32[j][i];
-              no += image->data.F32[j][i];
-            }
-          }
-          xo /= no;
-          yo /= no;
-          xo2 = sqrt (xo2/no - xo*xo);
-          yo2 = sqrt (yo2/no - yo*yo);
-          fprintf (stderr, "pre-shift centroid: %f,%f, sigma: %f,%f: flux: %f\n", xo, yo, xo2, yo2, no);
-        }
-
-// measure centroid of unshifted gaussian (should be 16.0,16.0)
-        {
-          psImage *image = flux;
-          float xo = 0.0;
-          float yo = 0.0;
-          float xo2 = 0.0;
-          float yo2 = 0.0;
-          float no = 0.0;
-          for (int j = 0; j < image->numRows; j++)
-          {
-            for (int i = 0; i < image->numCols; i++) {
-              xo += i*image->data.F32[j][i];
-              yo += j*image->data.F32[j][i];
-              xo2 += i*i*image->data.F32[j][i];
-              yo2 += j*j*image->data.F32[j][i];
-              no += image->data.F32[j][i];
-            }
-          }
-          xo /= no;
-          yo /= no;
-          xo2 = sqrt (xo2/no - xo*xo);
-          yo2 = sqrt (yo2/no - yo*yo);
-          fprintf (stderr, "pre-shift centroid: %f,%f, sigma: %f,%f: flux: %f\n", xo, yo, xo2, yo2, no);
-        }
-# endif