Changeset 6949 for trunk/psphot/src/psphotSourceFits.c

Timestamp:

Apr 21, 2006, 11:50:22 AM (20 years ago)

Author:

eugene

Message:

added tests for NaN results

File:

• trunk/psphot/src/psphotSourceFits.c (modified) (7 diffs)

trunk/psphot/src/psphotSourceFits.c

@@ -10 +10 @@
     // save the PSF model from the Ensemble fit
     pmModel *PSF = pmModelCopy (source->modelPSF);
+    if (isnan(PSF->params->data.F32[1])) psAbort ("psphot", "nan in psf fit");
 
     // extend source radius as needed
@@ -191 +192 @@
     PSF = source->modelPSF;
     psphotCheckRadiusPSFBlend (readout, source, PSF, 8.0);
+    if (isnan(PSF->params->data.F32[1])) psAbort ("psphot", "nan in dbl fit");
 
     modelSet = psArrayAlloc (2);
-    modelSet->n = 0;
+    // DROP modelSet->n = 0;
 
     DBL = pmModelCopy (PSF);
@@ -224 +226 @@
     // use the source moments, etc to guess basic model parameters
     pmModel *EXT = pmSourceModelGuess (source, modelTypeEXT); 
+    // if (isnan(EXT->params->data.F32[1])) psAbort ("psphot", "nan in ext fit");
 
     psphotCheckRadiusEXT (readout, source, EXT);
@@ -229 +232 @@
     x = EXT->params->data.F32[2];
     y = EXT->params->data.F32[3];
+
+    if ((source->moments->Sx < 1e-3) || (source->moments->Sx < 1e-3)) {
+        psTrace ("psphotSourceFits", 5, "problem source: moments: %f %f\n", source->moments->Sx, source->moments->Sy);
+    }
 
     // fit EXT (not PSF) model (set/unset the pixel mask)
@@ -251 +258 @@
     // save the PSF model from the Ensemble fit
     pmModel *PSF = pmModelCopy (source->modelPSF);
+    if (isnan(PSF->params->data.F32[1])) psAbort ("psphot", "nan in blend fit primary");
 
     x = PSF->params->data.F32[2];
@@ -256 +264 @@
 
     psArray *modelSet = psArrayAlloc (source->blends->n + 1);
-    modelSet->n = 0;
+    // DROP modelSet->n = 0;
     psArrayAdd (modelSet, 16, PSF);
     
     psArray *sourceSet = psArrayAlloc (source->blends->n + 1);
-    sourceSet->n = 0;
+    // DROP sourceSet->n = 0;
     psArrayAdd (sourceSet, 16, source);
 
@@ -282 +290 @@
         // XXX assume local sky is 0.0?
         model->params->data.F32[1] = blend->moments->Peak - blend->moments->Sky;
+        if (isnan(model->params->data.F32[1])) psAbort ("psphot", "nan in blend fit");
         model->params->data.F32[2] = blend->peak->x;
         model->params->data.F32[3] = blend->peak->y;