Changeset 12987

Timestamp:

Apr 24, 2007, 8:40:45 AM (19 years ago)

Author:

eugene

Message:

changed pmModelSub to pmSourceSub in non-linear fits

File:

• branches/eam_02_branch/psphot/src/psphotSourceFits.c (modified) (6 diffs)

branches/eam_02_branch/psphot/src/psphotSourceFits.c

@@ -108 +108 @@
         blend->modelPSF = psMemIncrRefCounter (model);
         psTrace ("psphot", 5, "fitted blend as PSF\n");
-        pmModelSub (source->pixels, source->maskObj, model, PM_MODEL_OP_FULL);
+        pmSourceCacheModel (source);
+        pmSourceSub (source, PM_MODEL_OP_FULL);
         blend->mode |=  PM_SOURCE_MODE_SUBTRACTED;
         blend->mode &= ~PM_SOURCE_MODE_TEMPSUB;
@@ -126 +127 @@
 
     psTrace ("psphot", 5, "fitted primary as PSF\n");
-    pmModelSub (source->pixels, source->maskObj, PSF, PM_MODEL_OP_FULL);
+    pmSourceCacheModel (source);
+    pmSourceSub (source, PM_MODEL_OP_FULL);
     psFree (source->modelPSF);
     source->modelPSF = PSF;
@@ -161 +163 @@
 
     psTrace ("psphot", 5, "fitted as PSF\n");
-    // should this be pmSourceSub?
-    pmModelSub (source->pixels, source->maskObj, PSF, PM_MODEL_OP_FULL);
+    pmSourceCacheModel (source);
+    pmSourceSub (source, PM_MODEL_OP_FULL);
 
     // free old model, save new model
@@ -257 +259 @@
     // sub EXT
     psFree (DBL);
-    pmModelSub (source->pixels, source->maskObj, EXT, PM_MODEL_OP_FULL);
-    psTrace ("psphot", 5, "blob as EXT: %f %f\n", EXT->params->data.F32[PM_PAR_XPOS], EXT->params->data.F32[PM_PAR_YPOS]);
 
     // save new model
     source->modelEXT = EXT;
+
+    // XXX need to set model before subtracting
+    pmSourceCacheModel (source);
+    pmSourceSub (source, PM_MODEL_OP_FULL);
+    psTrace ("psphot", 5, "blob as EXT: %f %f\n", EXT->params->data.F32[PM_PAR_XPOS], EXT->params->data.F32[PM_PAR_YPOS]);
+
     source->mode |=  PM_SOURCE_MODE_SUBTRACTED;
     source->mode &= ~PM_SOURCE_MODE_TEMPSUB;
@@ -269 +275 @@
     // sub DLB
     psFree (EXT);
-    pmModelSub (source->pixels, source->maskObj, (pmModel *) DBL->data[0], PM_MODEL_OP_FULL);
-    pmModelSub (source->pixels, source->maskObj, (pmModel *) DBL->data[1], PM_MODEL_OP_FULL);
-    psTrace ("psphot", 5, "blob as DBL: %f %f\n", ONE->params->data.F32[PM_PAR_XPOS], ONE->params->data.F32[PM_PAR_YPOS]);
+
+    // XXX double check this
+    // pmModelSub (source->pixels, source->maskObj, (pmModel *) DBL->data[0], PM_MODEL_OP_FULL);
+    // pmModelSub (source->pixels, source->maskObj, (pmModel *) DBL->data[1], PM_MODEL_OP_FULL);
 
     // drop old model, save new second model...
@@ -290 +297 @@
     newSrc->type     = source->type;
     newSrc->mode     = source->mode;
+
+    replace with:
+      pmSourceSub (source, PM_MODEL_OP_FULL);
+      pmSourceSub (newSrc, PM_MODEL_OP_FULL);
+
+    psTrace ("psphot", 5, "blob as DBL: %f %f\n", ONE->params->data.F32[PM_PAR_XPOS], ONE->params->data.F32[PM_PAR_YPOS]);
+
     psArrayAdd (sources, 100, newSrc);
     psFree (newSrc);