Changeset 1199 for tags/Release2.1.0b/src

Timestamp:

08/29/06 11:00:24 (18 years ago)

Author:

mar637

Message:

changed lag_flag interface to use frequency width intead of lag width

Location:

tags/Release2.1.0b/src

Files:

• STMath.cpp (modified) (3 diffs)
• STMath.h (modified) (1 diff)
• STMathWrapper.h (modified) (1 diff)

tags/Release2.1.0b/src/STMath.cpp

@@ -1381 +1381 @@
 CountedPtr< Scantable >
   asap::STMath::lagFlag( const CountedPtr< Scantable > & in,
-                          double frequency, int width )
+                          double frequency, double width )
 {
   CountedPtr< Scantable > out = getScantable(in, false);
@@ -1399 +1399 @@
       Vector<Float> spec = specCol(i);
       Vector<uChar> flag = flagCol(i);
-      Int lag = Int(spec.nelements()*abs(inc)/frequency);
+      Int lag0 = Int(spec.nelements()*abs(inc)/(frequency+width)+0.5);
+      Int lag1 = Int(spec.nelements()*abs(inc)/(frequency-width)+0.5);
       for (int k=0; k < flag.nelements(); ++k ) {
         if (flag[k] > 0) {
@@ -1407 +1408 @@
       Vector<Complex> lags;
       ffts.fft(lags, spec);
-      Int start =  max(0, lag-width);
-      Int end =  min(Int(lags.nelements()-1), lag+width);
+      Int start =  max(0, lag0);
+      Int end =  min(Int(lags.nelements()-1), lag1);
       if (start == end) {
         lags[start] = Complex(0.0);

tags/Release2.1.0b/src/STMath.h

@@ -191 +191 @@
   casa::CountedPtr<Scantable>
     lagFlag( const casa::CountedPtr<Scantable>& in, double frequency,
-              int width);
+              double width);
 
 private:

tags/Release2.1.0b/src/STMathWrapper.h

@@ -157 +157 @@
 
   ScantableWrapper lagFlag( const ScantableWrapper& in,
-                            double frequency, int width )
+                            double frequency, double width )
   { return ScantableWrapper(STMath::lagFlag(in.getCP(), frequency, width)); }