Changeset 31990 for branches/eam_branches/ipp-20110710/psphot/src/psphotSersicModelClass.c

Timestamp:

Aug 3, 2011, 10:30:02 AM (15 years ago)

Author:

eugene

Message:

add new window function using first radial moment; use new sersic guess function

File:

• branches/eam_branches/ipp-20110710/psphot/src/psphotSersicModelClass.c (modified) (5 diffs)

branches/eam_branches/ipp-20110710/psphot/src/psphotSersicModelClass.c

@@ -505 +505 @@
 }
 
-void psphotWriteGuess(pmModel *model, float dKronMag, float PSFratio, float ASPratio);
+# define TIMING false
 
 bool psphotSersicModelClassGuessPCM (pmPCMdata *pcm, pmSource *source) {
+
+    float t1, t2, t4, t5;
+    if (TIMING) { psTimerStart ("SersicGuess"); }
 
     psF32 *PAR = pcm->modelConv->params->data.F32;
@@ -549 +552 @@
 
     // find the containing model class:
+
+    if (TIMING) { t1 = psTimerMark ("SersicGuess"); }
 
     psphotSersicModelClass *class = NULL;    
@@ -561 +566 @@
     }
         
+    if (TIMING) { t2 = psTimerMark ("SersicGuess"); }
+
     psAssert (class, "PSFratio and ASPratio must be in range");
 
@@ -610 +617 @@
     PAR[PM_PAR_I0]  = 1.0;
 
+    if (TIMING) { t4 = psTimerMark ("SersicGuess"); }
+
     // set the normalization by linear fit between model and data
     psphotSersicModelNorm (pcm, source);
+
+    if (TIMING) { t5 = psTimerMark ("SersicGuess"); }
+
+    if (TIMING) {
+        fprintf (stderr, "guess, t1: %6.4f, t2: %6.4f, t4: %6.4f, t5: %6.4f\n", t1, t2, t4, t5);
+    }
 
     // float flux = pcm->modelConv->modelFlux(pcm->modelConv->params);
@@ -617 +632 @@
     // float Io = pow(10.0, -0.4*(totalMag - normMag));
     // PAR[PM_PAR_I0] = Io;
-
-    psphotWriteGuess(pcm->modelConv, dKronMag, PSFratio, ASPratio);
 
     return true;