Changeset 4901 for trunk/psphot/src/fit_galaxies.c

Timestamp:

Aug 29, 2005, 5:44:55 PM (21 years ago)

Author:

eugene

Message:

various psphot updates

File:

• trunk/psphot/src/fit_galaxies.c (modified) (3 diffs)

trunk/psphot/src/fit_galaxies.c

@@ -24 +24 @@
 
     psTimerStart ("psphot");
+
+    float dt = 0;
+    float dT = 0;
+    float dT2 = 0;
+    int dTn = 0;
+
     for (int i = 0; i < sources->n; i++) {
         psSource *source = sources->data[i];
@@ -72 +78 @@
         // fit as FLT, not PSF (skip poor fits)
         psImageKeepCircle (source->mask, x, y, model->radius, OR, 0x80);
+        psTimerStart ("galfit");
         status = pmSourceFitModel_EAM (source, model, false);
+        dt = psTimerMark ("galfit");
+        dT += dt;
+        dT2 += dt*dt;
+        dTn ++;
         psImageKeepCircle (source->mask, x, y, model->radius, AND, 0x7f);
         if (!status) {
@@ -105 +116 @@
 
     }
+    dt = dT / dTn;
+    fprintf (stderr, "dt: %f +/- %f\n", dt, sqrt(dT2/dTn - PS_SQR(dt)));
     psLogMsg ("psphot", 3, "fit galaxies: %f sec for %d galaxies (%d failures, %d total iterations)\n", psTimerMark ("psphot"), Nfit, Nfail, Niter);
     return (true);