Changeset 19823 for trunk/Ohana/src/libohana

Timestamp:

Oct 2, 2008, 8:08:31 AM (18 years ago)

Author:

eugene

Message:

replace poor code with ohana_normalize_angle

Location:

trunk/Ohana/src/libohana

Files:

• include/ohana.h (modified) (1 diff)
• src/time.c (modified) (1 diff)

trunk/Ohana/src/libohana/include/ohana.h

@@ -203 +203 @@
 double  ohana_lst              PROTO((double jd, double longitude));
 int     ohana_str_to_radec     PROTO((double *ra, double *dec, char *str1, char *str2));
+double  ohana_normalize_angle  PROTO((double angle));
 
 int     hstgsc_hms_to_deg      PROTO((double *h0, double *h1, double *d0, double *d1, char *string));

trunk/Ohana/src/libohana/src/time.c

@@ -434 +434 @@
   return (sid_g);
 }
+
+double ohana_normalize_angle (double angle) {
+
+    double result;
+
+    // take an input angle and force the domain to be 0.0 - 360.0
+    // there are a few ways to do this:  
+
+    // option 1: This is a potentially very slow method, also subject to round off errors
+# if (0)
+    while (angle > 360.0) angle -= 360.0;
+    while (angle <   0.0) angle += 360.0;
+# endif
+    
+    // option 2: take sin & cos, apply atan2 (y, x) 
+# if (1)    
+    double x, y;
+
+    x = cos(angle*RAD_DEG);
+    y = sin(angle*RAD_DEG);
+
+    result = atan2 (y, x);
+    if (result < 0.0) result += 360.0;
+# endif
+
+# if (0)
+    // option 3:
+    int nCircle = angle / 360.0;
+    result -= 360.0*nCircle;
+    if (result < 0.0) result += 360.0;
+# endif
+
+    return (result);
+}