Changeset 34057 for branches/eam_branches/ipp-20120601/ppSim/src/ppSimInsertStars.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/ppSimInsertStars.c (modified) (3 diffs)

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

@@ -6 +6 @@
     PTR = NULL;
 
-float imageSum (psImage *image);
+double imageSum (psImage *image);
 
 bool ppSimInsertStars (pmReadout *readout, psImage *expCorr, psArray *stars, pmConfig *config) {
@@ -111 +111 @@
 
         // XXX test:
-        flux = star->flux * expCorr->data.F32[0][0];
+        // flux = star->flux * expCorr->data.F32[0][0];
 
         // instantiate a model for the PSF at this location, set desired flux
@@ -145 +145 @@
             readout->image->data.F32[(int)(yCell)][(int)(xCell)] += flux;
         } else {
-            // insert the source flux in the image
+            // insert the source flux in the image and measure inserted flux
+            // this is 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 sum1 = imageSum(source->pixels);
             pmSourceAddWithOffset (source, PM_MODEL_OP_FULL, 0xff, dX, dY);