Changeset 19823 for trunk/Ohana/src/opihi/cmd.astro

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/opihi/cmd.astro

Files:

• cplot.c (modified) (1 diff)
• czplot.c (modified) (1 diff)
• drizzle.c (modified) (1 diff)
• getvel.c (modified) (1 diff)
• polar.c (modified) (1 diff)
• rotcurve.c (modified) (1 diff)

trunk/Ohana/src/opihi/cmd.astro/cplot.c

@@ -41 +41 @@
   Npts = 0;
   for (i = 0; i < Xvec.Nelements; i++, r++, d++) {
+      *r = ohana_normalize_angle (*r);
     while (*r < Rmin) *r += 360.0;
     while (*r > Rmax) *r -= 360.0;

trunk/Ohana/src/opihi/cmd.astro/czplot.c

@@ -49 +49 @@
   for (i = 0; i < Zvec.Nelements; i++, in++, out++, r++, d++, x++, y++) {
     *out = MIN (1.0, MAX (0.01, (*in - min) / range));
-    while (*r < Rmin) *r += 360.0;
-    while (*r > Rmax) *r -= 360.0;
+    *r = ohana_normalize_angle (*r);
     fRD_to_XY (x, y, *r, *d, &graphmode.coords);
   }

trunk/Ohana/src/opihi/cmd.astro/drizzle.c

@@ -221 +221 @@
   *phi   = DEG_RAD*atan2(y,x) + rot_phi;
   
-  while (*phi <   0.0) *phi += 360.0;
-  while (*phi > 360.0) *phi -= 360.0;
+  *phi = ohana_normalize_angle(*phi);
   return (TRUE);
 }

trunk/Ohana/src/opihi/cmd.astro/getvel.c

@@ -43 +43 @@
   dV *= 0.001;
 
-  while (L >= 360) {L -= 360.0;}
-  while (L < 0.0)  {L += 360.0;}
+  L = ohana_normalize_angle (L);
   gprint (GP_ERR, "L: %f\n", L);
 

trunk/Ohana/src/opihi/cmd.astro/polar.c

@@ -55 +55 @@
       y = (j - Yo) * dN + No;
       r = hypot(x, y);
-      t = DEG_RAD*atan2 (y, x);
-      while (t < 360.0) {t += 360.0;}
-      while (t > 360.0) {t -= 360.0;}
+      t = ohana_normalize_angle(DEG_RAD*atan2 (y, x));
       X = (t - Lo) / dL + xo;
       Y = (r - Do) / dD + yo;

trunk/Ohana/src/opihi/cmd.astro/rotcurve.c

@@ -70 +70 @@
   gfits_scan (&out[0].header, "CRPIX2", "%lf", 1, &Yo);
 
-  while (L >= 360) {L -= 360.0;}
-  while (L < 0.0)  {L += 360.0;}
+  L = ohana_normalize_angle (L);
+
   X = (L - Lo) / dL + Xo;
   if ((X >= Nx) || (X < 0)) {