Changeset 1320 for branches/Release2.1.2/src/RowAccumulator.cpp

Timestamp:

03/02/07 11:02:27 (17 years ago)

Author:

mar637

Message:

merge from trunk, to get most recent bug fixes

File:

• branches/Release2.1.2/src/RowAccumulator.cpp (modified) (4 diffs)

branches/Release2.1.2/src/RowAccumulator.cpp

@@ -13 +13 @@
 
 #include <casa/Arrays/MaskArrMath.h>
+#include <casa/Arrays/MaskArrLogi.h>
 #include <casa/Arrays/ArrayMath.h>
+#include <casa/Arrays/ArrayLogical.h>
 #include "RowAccumulator.h"
 
@@ -40 +42 @@
   if (!initialized_) {
     Vector<Float> dummy(v.nelements(), 0.0);
-    spectrum_.setData(dummy, Vector<Bool>(m.nelements(), True));
-    n_.setData(Vector<Float>(v.nelements(), 0.0), m);
-    weightSum_.setData(Vector<Float>(v.nelements(), 0.0), m);
+    Vector<Bool> dummymsk(m.nelements(), True);
+    spectrum_.setData(dummy, dummymsk);
+    n_.setData(Vector<Float>(v.nelements(), 0.0), dummymsk);
+    weightSum_.setData(Vector<Float>(v.nelements(), 0.0), dummymsk);
     tsysSum_.resize(tsys.nelements()); tsysSum_=0.0;
   }
@@ -48 +51 @@
   Float totalweight = 1.0;
   totalweight *= addTsys(tsys);
-  totalweight *= addInterval(interval);
-  addTime(time);
+  // only add these if not everything masked
+  if ( !allEQ(m, False) ) {
+    totalweight *= addInterval(interval);
+    addTime(time);
+  }
   addSpectrum(v, m, totalweight);
   initialized_ = True;
@@ -134 +140 @@
 casa::Vector< casa::Float > asap::RowAccumulator::getTsys( ) const
 {
-  // @fixme this assummes tsys.nelements() == 1
+  // @fixme this assumes tsys.nelements() == 1
   return tsysSum_/max(n_);
 }