Changeset 36240

Timestamp:

Oct 26, 2013, 2:56:36 PM (13 years ago)

Author:

eugene

Message:

add pcm->poissonErrors; if not(poissonErrors), set variance to a constant in pmPCM_Minimize.c

Location:

branches/eam_branches/ipp-20130904/psModules/src/objects

Files:

• pmPCM_MinimizeChisq.c (modified) (1 diff)
• pmPCMdata.c (modified) (2 diffs)
• pmPCMdata.h (modified) (1 diff)

branches/eam_branches/ipp-20130904/psModules/src/objects/pmPCM_MinimizeChisq.c

@@ -580 +580 @@
 
             // XXXX note this point here:::
-            float yweight = 1.0 / source->variance->data.F32[i][j];
+            float yweight = pcm->poissonErrors ? 1.0 / source->variance->data.F32[i][j] : 1.0;
             float delta = ymodel - source->pixels->data.F32[i][j];
 

branches/eam_branches/ipp-20130904/psModules/src/objects/pmPCMdata.c

@@ -95 +95 @@
     pcm->nDOF = 0;
 
+    pcm->poissonErrors = true;
+
     // full convolution with the PSF is expensive.  if we have to save time, we can do a 1D
     // convolution with a Gaussian approximation to the kernel
@@ -365 +367 @@
     pcm->modelConv = psMemIncrRefCounter(model);
     pcm->constraint = constraint;
+
+    pcm->poissonErrors = fitOpions->poissonErrors;
 
     pcm->nPix = nPix;

branches/eam_branches/ipp-20130904/psModules/src/objects/pmPCMdata.h

@@ -35 +35 @@
     int nPar;
     int nDOF;
+
+    bool poissonErrors;
 
     bool use1Dgauss;