Changeset 139

Timestamp:

12/24/04 13:21:30 (20 years ago)

Author:

mar637

Message:

Added inSitu version of multiply/scale

Location:

trunk/src

Files:

• SDMath.cc (modified) (1 diff)
• SDMath.h (modified) (1 diff)
• SDMathWrapper.h (modified) (1 diff)

trunk/src/SDMath.cc

@@ -204 +204 @@
 }
 
+void SDMath::multiplyInSitu(SDMemTable* in, Float factor) {
+  SDMemTable* sdmt = new SDMemTable(*in);
+  Table t = sdmt->table();
+  ArrayColumn<Float> spec(t,"SPECTRA");  
+  for (uInt i=0; i < t.nrow(); i++) {
+    MaskedArray<Float> marr(sdmt->rowAsMaskedArray(i));
+    marr *= factor;
+    spec.put(i, marr.getArray());
+  }
+  in = sdmt;
+  delete sdmt;sdmt=0;
+}
+
 CountedPtr<SDMemTable>
 SDMath::multiply(const CountedPtr<SDMemTable>& in, Float factor) 

trunk/src/SDMath.h

@@ -46 +46 @@
   casa::CountedPtr<SDMemTable> quotient(const casa::CountedPtr<SDMemTable>& on, 
                                          const casa::CountedPtr<SDMemTable>& off);
+  void multiplyInSitu(SDMemTable* in, casa::Float factor);
+
   casa::CountedPtr<SDMemTable> multiply(const casa::CountedPtr<SDMemTable>& in, 
                                   casa::Float factor);

trunk/src/SDMathWrapper.h

@@ -57 +57 @@
     return SDMemTableWrapper(SDMath::multiply(in.getCP(),factor));
   }
+  void scaleInSitu(SDMemTableWrapper& in, casa::Float factor) {
+    SDMemTable* sdmt = in.getPtr();
+    SDMath::multiplyInSitu(in.getPtr(),factor);
+  }
 
   SDMemTableWrapper add(const SDMemTableWrapper& in,