Changeset 1199 for tags/Release2.1.0b/src/STMath.cpp

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

File:

• tags/Release2.1.0b/src/STMath.cpp (modified) (3 diffs)

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);