Changeset 13335 for trunk/psModules/src/objects/pmSource.c

Timestamp:

May 10, 2007, 9:48:35 AM (19 years ago)

Author:

magnier

Message:

adding PM_MODEL_OP_NOISE

File:

• trunk/psModules/src/objects/pmSource.c (modified) (4 diffs)

trunk/psModules/src/objects/pmSource.c

@@ -6 +6 @@
  *  @author EAM, IfA: significant modifications.
  *
- *  @version $Revision: 1.27 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2007-04-26 01:20:29 $
+ *  @version $Revision: 1.28 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2007-05-10 19:48:35 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -799 +799 @@
         }
 
+        psF32 **target = source->pixels->data.F32;
+        if (mode & PM_MODEL_OP_NOISE) {
+          target = source->noise->data.F32;
+        }
+
         // XXX need to respect the source and model masks?
         for (int iy = 0; iy < source->modelFlux->numRows; iy++) {
@@ -807 +812 @@
                 float value = Io*source->modelFlux->data.F32[iy][ix];
                 if (add) {
-                    source->pixels->data.F32[oy][ox] += value;
+                    target[oy][ox] += value;
                 } else {
-                    source->pixels->data.F32[oy][ox] -= value;
+                    target[oy][ox] -= value;
                 }
             }
@@ -816 +821 @@
     } 
 
+    psImage *target = source->pixels;
+    if (mode & PM_MODEL_OP_NOISE) {
+        target = source->noise;
     if (add) {
-        status = pmModelAdd (source->pixels, source->maskObj, model, PM_MODEL_OP_FULL);
+        status = pmModelAdd (target, source->maskObj, model, PM_MODEL_OP_FULL);
     } else {
-        status = pmModelSub (source->pixels, source->maskObj, model, PM_MODEL_OP_FULL);
+        status = pmModelSub (target, source->maskObj, model, PM_MODEL_OP_FULL);
     }