Changeset 20582

Timestamp:

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

Author:

eugene

Message:

adding pmSourceVisual

Location:

trunk/psModules/src

Files:

• objects/Makefile.am (modified) (2 diffs)
• objects/pmPSFtry.c (modified) (8 diffs)
• objects/pmSourceVisual.c (added)
• objects/pmSourceVisual.h (added)
• psmodules.h (modified) (1 diff)

trunk/psModules/src/objects/Makefile.am

@@ -39 +39 @@
      pmSourcePlotMoments.c \
      pmSourcePlotApResid.c \
+     pmSourceVisual.c \
      pmResiduals.c \
      pmPSF.c \
@@ -74 +75 @@
      pmSourceIO.h \
      pmSourcePlots.h \
+     pmSourceVisual.h \
      pmResiduals.h \
      pmPSF.h \

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

trunk/psModules/src/psmodules.h

@@ -121 +121 @@
 #include <pmModelUtils.h>
 #include <pmSourcePhotometry.h>
+#include <pmSourceVisual.h>
 
 // The following headers are from random locations, here because they cross bounds