Changeset 20582 for trunk/psModules/src/objects/pmPSFtry.c

Timestamp:

Nov 7, 2008, 3:52:34 PM (18 years ago)

Author:

eugene

Message:

adding pmSourceVisual

File:

• trunk/psModules/src/objects/pmPSFtry.c (modified) (8 diffs)

trunk/psModules/src/objects/pmPSFtry.c

@@ -5 +5 @@
  *  @author EAM, IfA
  *
- *  @version $Revision: 1.65 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2008-10-29 01:04:27 $
+ *  @version $Revision: 1.66 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2008-11-08 01:52:34 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -34 +34 @@
 #include "pmSourceFitModel.h"
 #include "pmSourcePhotometry.h"
+#include "pmSourceVisual.h"
 
 bool printTrendMap (pmTrend2D *trend) {
@@ -683 +684 @@
             stdev = psStatsGetValue (trend->stats, stdevOption);
             psTrace ("psModules.objects", 4, "clipped E0 : %f +/- %f keeping %ld of %ld\n", mean, stdev, psf->psfTrendStats->clippedNvalues, e0->n);
+            pmSourceVisualPSFModelResid (trend, x, y, e0, srcMask);
 
             trend = psf->params->data[PM_PAR_E1];
@@ -689 +691 @@
             stdev = psStatsGetValue (trend->stats, stdevOption);
             psTrace ("psModules.objects", 4, "clipped E1 : %f +/- %f keeping %ld of %ld\n", mean, stdev, psf->psfTrendStats->clippedNvalues, e1->n);
+            pmSourceVisualPSFModelResid (trend, x, y, e1, srcMask);
 
             trend = psf->params->data[PM_PAR_E2];
@@ -695 +698 @@
             stdev = psStatsGetValue (trend->stats, stdevOption);
             psTrace ("psModules.objects", 4, "clipped E2 : %f +/- %f keeping %ld of %ld\n", mean, stdev, psf->psfTrendStats->clippedNvalues, e2->n);
+            pmSourceVisualPSFModelResid (trend, x, y, e2, srcMask);
 
             if (!status) {
@@ -853 +857 @@
         psImageMapCleanup (trend->map);
         // printTrendMap (trend);
+        pmSourceVisualPSFModelResid (trend, x, y, e0obs, mask);
 
         trend = psf->params->data[PM_PAR_E1];
@@ -862 +867 @@
         psImageMapCleanup (trend->map);
         // printTrendMap (trend);
+        pmSourceVisualPSFModelResid (trend, x, y, e1obs, mask);
 
         trend = psf->params->data[PM_PAR_E2];
@@ -871 +877 @@
         psImageMapCleanup (trend->map);
         // printTrendMap (trend);
+        pmSourceVisualPSFModelResid (trend, x, y, e2obs, mask);
     }
     psf->psfTrendStats->clipIter = nIter; // restore default setting