Changeset 31381

Timestamp:

Apr 25, 2011, 5:26:32 PM (15 years ago)

Author:

eugene

Message:

psphotEff needs to set the fit and aperture radii; do not report an error unless the flux is nan; add test code to dump the sparse matrix data

Location:

branches/eam_branches/ipp-20110404/psphot/src

Files:

• psphotEfficiency.c (modified) (1 diff)
• psphotFitSourcesLinear.c (modified) (1 diff)
• psphotMagnitudes.c (modified) (1 diff)

branches/eam_branches/ipp-20110404/psphot/src/psphotEfficiency.c

@@ -400 +400 @@
                 source->type = PM_SOURCE_TYPE_STAR;
 
+                source->modelPSF->fitRadius = sourceRadius;
+                source->apRadius = sourceRadius;
+
                 numFound++;
                 psArrayAdd(sources, sources->n, source);

branches/eam_branches/ipp-20110404/psphot/src/psphotFitSourcesLinear.c

@@ -267 +267 @@
         }
     }
+
+# if (0)
+    static int npass = 0;
+    char name[128];
+    FILE *f1 = NULL;
+    int fd = -1;
+
+    snprintf (name, 128, "sparse.Aij.%02d.dat", npass);
+    f1 = fopen (name, "w");
+    psAssert (f1, "failed to open file\n");
+    fd = fileno (f1);
+    p_psVectorPrint (fd, sparse->Aij, "Aij");
+    fclose (f1);
+
+    snprintf (name, 128, "sparse.Bfj.%02d.dat", npass);
+    f1 = fopen (name, "w");
+    psAssert (f1, "failed to open file\n");
+    fd = fileno (f1);
+    p_psVectorPrint (fd, sparse->Bfj, "Bfj");
+    fclose (f1);
+
+    snprintf (name, 128, "sparse.Qii.%02d.dat", npass);
+    f1 = fopen (name, "w");
+    psAssert (f1, "failed to open file\n");
+    fd = fileno (f1);
+    p_psVectorPrint (fd, sparse->Qii, "Qii");
+    fclose (f1);
+    npass ++;
+# endif
 
     psSparseResort (sparse);

branches/eam_branches/ipp-20110404/psphot/src/psphotMagnitudes.c

@@ -210 +210 @@
 
         status = pmSourceMagnitudes (source, psf, photMode, maskVal, markVal, source->apRadius);
-        if (status && isfinite(source->apMag)) {
+        if (status && isfinite(source->apFlux)) {
             Nap ++;
         } else {