Changeset 31381 for branches/eam_branches/ipp-20110404/psphot/src/psphotFitSourcesLinear.c

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

File:

• branches/eam_branches/ipp-20110404/psphot/src/psphotFitSourcesLinear.c (modified) (1 diff)

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);