Changeset 2461 for trunk/src

Timestamp:

04/11/12 12:32:41 (13 years ago)

Author:

Takeshi Nakazato

Message:

New Development: No

JIRA Issue: No

Ready for Test: Yes

Interface Changes: No

What Interface Changed: Please list interface changes

Test Programs: sdgrid unit test

Put in Release Notes: No

Module(s): Module Names change impacts.

Description: Describe your changes here...

Fixed a bug which is caused by the fact that R.A. can be shifted by 2*pi
when it is represented in radian. Center coordinate is shifted by 2*pi
or -2*pi to match DIRECTION value if it is too far from observed position.

File:

• trunk/src/STGrid.cpp (modified) (1 diff)

trunk/src/STGrid.cpp

@@ -1058 +1058 @@
     center_(0) = xcen.getValue( "rad" ) ;
     center_(1) = ycen.getValue( "rad" ) ;
-  }
+    double base = 0.5 * (xmin + xmax) ;
+    int maxrotate = 1 ;
+    int nelem = 2 * maxrotate + 1 ;
+    double *sep = new double[nelem] ;
+    for ( int i = 0 ; i < nelem ; i++ )
+      sep[i] = abs(base - center_[0] - (i-maxrotate) * C::_2pi) ;
+//     os << "sep[0]=" << sep[0] << endl  
+//        << "sep[1]=" << sep[1] << endl
+//        << "sep[2]=" << sep[2] << LogIO::POST ;
+    int idx = 0 ;
+    base = sep[0] ;
+    int nrotate = 0 ;
+    while ( idx < nelem ) {
+      if ( base > sep[idx] ) {
+        base = sep[idx] ;
+        nrotate = idx ;
+      }
+      idx++ ;
+    }
+    delete sep ;
+    nrotate -= maxrotate ;
+//     os << "nrotate = " << nrotate << LogIO::POST ;
+    center_[0] += nrotate * C::_2pi ;
+  }
+//   os << "xmin=" << xmin << LogIO::POST ;
+//   os << "center_=" << center_ << LogIO::POST ;