Index: trunk/psLib/src/astronomy/psCoord.c
===================================================================
--- trunk/psLib/src/astronomy/psCoord.c	(revision 2650)
+++ trunk/psLib/src/astronomy/psCoord.c	(revision 2778)
@@ -10,6 +10,6 @@
 *  @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,5 +314,5 @@
         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,9 +354,9 @@
     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,5 +375,5 @@
         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);