Changeset 34057 for branches/eam_branches/ipp-20120601/ppSim/src/ppSimInsertGalaxies.c

Timestamp:

Jun 23, 2012, 6:56:57 AM (14 years ago)

Author:

eugene

Message:

need to use double for the imageSum

File:

• branches/eam_branches/ipp-20120601/ppSim/src/ppSimInsertGalaxies.c (modified) (2 diffs)

branches/eam_branches/ipp-20120601/ppSim/src/ppSimInsertGalaxies.c

@@ -2 +2 @@
 static char *defaultModel = "PS_MODEL_SERSIC";
 
+// this is used to measure the inserted flux as a difference of 2 large numbers : must use
+// 'double' or we get floating point errors and trends on a scale of ~0.005 mags at > -10 mags
+// instrumental
 double imageSum (psImage *image) {
     double sum = 0.0;
@@ -165 +168 @@
         
         // insert the source flux in the image
-        float sum1 = imageSum(source->pixels);
+        double sum1 = imageSum(source->pixels);
         pmSourceAddWithOffset (source, PM_MODEL_OP_FULL, 0xff, dX, dY);
-        float sum2 = imageSum(source->pixels);
+        double sum2 = imageSum(source->pixels);
         float flux = sum2 - sum1;