Changeset 36210

Timestamp:

Oct 13, 2013, 1:32:54 PM (13 years ago)

Author:

eugene

Message:

some error handling, fixed error messages, call psphotInitRadiusEXT

File:

• branches/eam_branches/ipp-20130904/psphot/src/psphotGalaxyShape.c (modified) (4 diffs)

branches/eam_branches/ipp-20130904/psphot/src/psphotGalaxyShape.c

@@ -6 +6 @@
 
     fprintf (stdout, "\n");
-    psLogMsg ("psphot", PS_LOG_INFO, "--- psphot Kron Fluxes ---");
+    psLogMsg ("psphot", PS_LOG_INFO, "--- psphot Galaxy Shapes ---");
 
     // select the appropriate recipe information
@@ -58 +58 @@
 
     if (!sources->n) {
-        psLogMsg ("psphot", PS_LOG_INFO, "no sources, skipping kron fluxes");
+        psLogMsg ("psphot", PS_LOG_INFO, "no sources, skipping galaxy shapes");
         return true;
     }
 
     psTimerStart ("psphot.galaxy");
+
+    psphotInitRadiusEXT (recipe, readout);
 
     // determine the number of allowed threads
@@ -218 +220 @@
     // we are using fitOptions->mode : be sure this makes sense
     pmPCMdata *pcm = pmPCMinit (source, fitOptions, model, maskVal, psfSize);
+    if (!pcm) return false;
 
     // we are fitting only PM_PAR_I0; the shape elements are generated from a grid
@@ -302 +305 @@
     float Chisq = YY - 2 * Io * YM + Io * Io * MM;
 
-    fprintf (stderr, "Io: %f, Chisq: %f\n", Io, Chisq);
+    // fprintf (stderr, "Io: %f, Chisq: %f\n", Io, Chisq);
 
     psVectorAppend (source->galaxyFits->Io, Io);