Changeset 4716

Timestamp:

Aug 5, 2005, 9:56:18 AM (21 years ago)

Author:

eugene

Message:

fixed test for out-of-bounds AIT coords

File:

• trunk/Ohana/src/libohana/src/coordops.c (modified) (2 diffs)

trunk/Ohana/src/libohana/src/coordops.c

@@ -12 +12 @@
   int Zenith1, Zenith2, Zenithal, Polynomial, Cartesian, PseudoCyl;
   char *type;
-  double L, M, X, Y, T, Z;
+  double L, M, X, Y, T, Z, Z2;
   double R, sphi, cphi, stht, ctht;
   double alpha, delta, salp, calp, sdel, sdp, cdp;
@@ -113 +113 @@
   if (PseudoCyl) {
     if (!strcmp(type, "-AIT")) {
-      Z = sqrt (1.0 - SQ(RAD_DEG*0.25*L) - SQ(RAD_DEG*0.5*M));
-      alpha = 2.0 * DEG_RAD * atan2 (RAD_DEG*0.5*Z*L, 2.0*SQ(Z) - 1.0);
+      Z2 = (1.0 - SQ(RAD_DEG*0.25*L) - SQ(RAD_DEG*0.5*M));
+      if (Z2 < 0) return (FALSE);
+      Z = sqrt (Z2);
+      alpha = 2.0 * DEG_RAD * atan2 (RAD_DEG*0.5*Z*L, 2.0*Z2 - 1.0);
       delta = DEG_RAD * asin (RAD_DEG*M*Z);
       *ra  = alpha + coords[0].crval1;