Changeset 30776 for branches/eam_branches/ipp-20110213/psphot/src

Timestamp:

Mar 2, 2011, 5:44:48 PM (15 years ago)

Author:

eugene

Message:

free some data as early as possible (petProfile, source pixels); added memory dump tests

Location:

branches/eam_branches/ipp-20110213/psphot/src

Files:

• psphotExtendedSourceAnalysisByObject.c (modified) (1 diff)
• psphotMergeSources.c (modified) (1 diff)
• psphotRadialAperturesByObject.c (modified) (1 diff)
• psphotRadialProfile.c (modified) (5 diffs)
• psphotStack.c (modified) (2 diffs)
• psphotStackImageLoop.c (modified) (3 diffs)
• psphotStackReadout.c (modified) (4 diffs)

branches/eam_branches/ipp-20110213/psphot/src/psphotExtendedSourceAnalysisByObject.c

@@ -148 +148 @@
                 psFree(source->extpars->radFlux);
                 psFree(source->extpars->ellipticalFlux);
+                psFree(source->extpars->petProfile);
             }
         }

branches/eam_branches/ipp-20110213/psphot/src/psphotMergeSources.c

@@ -699 +699 @@
             sourceOut->parent = sourceSrc;
 
-            // keep the original source flags
+            // keep the original source flags and sequence ID (if set)
+            sourceOut->seq      = sourceSrc->seq;
             sourceOut->type     = sourceSrc->type;
             sourceOut->mode     = sourceSrc->mode;

branches/eam_branches/ipp-20110213/psphot/src/psphotRadialAperturesByObject.c

@@ -150 +150 @@
             } else {
                 source->mode |= PM_SOURCE_MODE_RADIAL_FLUX;
+                if (source->parent) {
+                    source->parent->mode |= PM_SOURCE_MODE_RADIAL_FLUX;
+                }
             }
 

branches/eam_branches/ipp-20110213/psphot/src/psphotRadialProfile.c

@@ -25 +25 @@
         return false;
     }
+    // allocate: extpars->radFlux->radii,fluxes,theta
 
     // use the radial profiles to determine the radius of a given isophote.  this isophote
@@ -33 +34 @@
         return false;
     }
+    // allocate : extpars->radFlux->isophotalRadii (use profile->radii,fluxes)
+
 
     // convert the isophotal radius vs angle measurements to an elliptical contour
@@ -39 +42 @@
         return false;
     }
-  
+    // use extpars->radFlux->isophotalRadii,theta (result in extpars->axes)
+
     // generate a single, normalized radial profile following the elliptical contours.
     // the radius is normalized by the axis ratio so that on the major axis, 1 pixel = 1 pixel
@@ -46 +50 @@
         return false;
     }
+    // allocate extpars->ellipticalFlux->radiusElliptical,fluxElliptical (use axes to scale raw pixels)
   
     // generated profile in averaged bins
@@ -52 +57 @@
         return false;
     }
+    // allocate extpars->radProfile->binSB, binSBstdv, binSum, binFill, radialBins, area (small lengths)
+    // use radiusElliptical, fluxElliptical, 
   
     return true;

branches/eam_branches/ipp-20110213/psphot/src/psphotStack.c

@@ -2 +2 @@
 
 int main (int argc, char **argv) {
+
+    // XXX add memory dumps
+    psMemDumpSetState(true);
 
     psTimerStart ("complete");
@@ -12 +15 @@
 
     psphotVersionPrint();
+
+    psMemDump("start");
 
     // load input data (config and images (signal, noise, mask)

branches/eam_branches/ipp-20110213/psphot/src/psphotStackImageLoop.c

@@ -18 +18 @@
     pmCell *cell;
     pmReadout *readout;
+
+    psMemDump("startloop");
 
     pmFPAfile *inputRaw = psMetadataLookupPtr (&status, config->files, "PSPHOT.STACK.INPUT.RAW");
@@ -86 +88 @@
         if (!pmFPAfileIOChecks (config, view, PM_FPA_AFTER)) ESCAPE ("failed output for Chip in psphot.");
     }
+    psMemDump("doneloop");
 
     SetAstrometryFPA(config, view, bilevelAstrometry);
@@ -96 +99 @@
 
     psFree (view);
+
+    psMemDump("doneoutput");
     return true;
 }

branches/eam_branches/ipp-20110213/psphot/src/psphotStackReadout.c

@@ -211 +211 @@
     }
 
+    psMemDump("psfstats");
+
     psphotStackObjectsUnifyPosition (objects);
 
@@ -230 +232 @@
     }
 
+    psMemDump("extmeas");
+
     bool smoothAgain = true;
     for (int nMatchedPSF = 0; smoothAgain; nMatchedPSF++) {
@@ -250 +254 @@
         // smooth to the next FWHM, or set 'smoothAgain' to false if no more 
         psphotStackMatchPSFsNext(&smoothAgain, config, view, STACK_OUT, nMatchedPSF);
+        psMemDump("matched");
     }
 
@@ -264 +269 @@
 
     // drop the references to the image pixels held by each source
-    // psphotSourceFreePixels (config, view, STACK_OUT);
+    psphotSourceFreePixels (config, view, STACK_OUT);
     psphotSourceFreePixels (config, view, STACK_SRC);