Changeset 2778 for trunk/psLib/src/astronomy/psCoord.c

Timestamp:

Dec 21, 2004, 10:42:07 AM (22 years ago)

Author:

gusciora

Message:

Added errors arguments to the VectorStats and VectorHistogram functions.
The arguments are ignored in the histogram and robust stats code.

File:

• trunk/psLib/src/astronomy/psCoord.c (modified) (4 diffs)

trunk/psLib/src/astronomy/psCoord.c

@@ -10 +10 @@
 *  @author GLG, MHPCC
 *
-*  @version $Revision: 1.40 $ $Name: not supported by cvs2svn $
-*  @date $Date: 2004-12-07 22:28:28 $
+*  @version $Revision: 1.41 $ $Name: not supported by cvs2svn $
+*  @date $Date: 2004-12-21 20:42:07 $
 *
 *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -314 +314 @@
         tmp->y = log(tan(DEG_TO_RAD(45.0) + (0.5 * coord->r)));
     } else if (projection->type == PS_PROJ_AIT) {
-        psF64 tmpF64  = sqrt(0.5 * (1.0 + (cos(coord->r) * cos(0.5 * coord->d))));
+        psF64 tmpF64  = PS_SQRT_F32(0.5 * (1.0 + (cos(coord->r) * cos(0.5 * coord->d))));
         tmpF64 = 1.0 / tmpF64;
         tmp->x = 2.0 * tmpF64 * cos(coord->r) * sin(0.5 * coord->d);
@@ -354 +354 @@
     if (projection->type == PS_PROJ_TAN) {
         tmp->r = atan2(-coord->y, coord->x);
-        R = sqrt((coord->x * coord->x) + (coord->y * coord->y));
+        R = PS_SQRT_F32((coord->x * coord->x) + (coord->y * coord->y));
         tmp->d = atan(1.0 / R);
 
     } else if (projection->type == PS_PROJ_SIN) {
-        R = sqrt((coord->x * coord->x) + (coord->y * coord->y));
+        R = PS_SQRT_F32((coord->x * coord->x) + (coord->y * coord->y));
         tmp->d = atan2(-coord->y, coord->x);
         tmp->r = acos((R * M_PI) / 180.0);
@@ -375 +375 @@
         chu2 = (coord->y * M_PI) / 360.0;
         chu2 *= chu2;
-        chu = sqrt(1.0 - chu1 - chu2);
+        chu = PS_SQRT_F32(1.0 - chu1 - chu2);
         tmp->d = 2.0 * atan2((2.0 * chu * chu) - 1.0, (coord->x * chu * M_PI) / 360.0);
         tmp->r = asin((coord->y * chu * M_PI) / 180.0);