Index: trunk/psLib/src/dataManip/psStats.c
===================================================================
--- trunk/psLib/src/dataManip/psStats.c	(revision 4138)
+++ trunk/psLib/src/dataManip/psStats.c	(revision 4162)
@@ -14,6 +14,6 @@
  *      stats->binsize
  *
- *  @version $Revision: 1.131 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2005-06-07 23:11:54 $
+ *  @version $Revision: 1.132 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2005-06-08 23:40:45 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -878,5 +878,5 @@
     } else {
         countFloat = (psF32)countInt;
-        stats->sampleStdev = PS_SQRT_F32((sumSquares - (sumDiffs * sumDiffs / countFloat)) / (countFloat - 1));
+        stats->sampleStdev = sqrtf((sumSquares - (sumDiffs * sumDiffs / countFloat)) / (countFloat - 1));
     }
 }
@@ -992,8 +992,8 @@
         // data ranges are used correctly.
         if (errors != NULL) {
-            stats->sampleStdev = (1.0 / PS_SQRT_F32(errorDivisor));
+            stats->sampleStdev = (1.0 / sqrtf(errorDivisor));
         } else {
             countFloat = (psF32)countInt;
-            stats->sampleStdev = PS_SQRT_F32((sumSquares - (sumDiffs * sumDiffs / countFloat)) / (countFloat - 1));
+            stats->sampleStdev = sqrtf((sumSquares - (sumDiffs * sumDiffs / countFloat)) / (countFloat - 1));
 
         }
@@ -1627,5 +1627,5 @@
         }
     }
-    myStdev = PS_SQRT_F32((sumSquares - (sumDiffs * sumDiffs / countFloat)) / (countFloat - 1));
+    myStdev = sqrtf((sumSquares - (sumDiffs * sumDiffs / countFloat)) / (countFloat - 1));
 
     p_psNormalizeVectorRangeF32(robustHistogramVector, 0.0, 1.0);
@@ -1657,5 +1657,5 @@
                      "WARNING: failed fit a 1D polynomial.\n");
         }
-        psF32 polyFitSigma = PS_SQRT_F32(-0.5 / tmpPoly->coeff[2]);
+        psF32 polyFitSigma = sqrtf(-0.5 / tmpPoly->coeff[2]);
         psF32 polyFitMean = tmpPoly->coeff[1] * PS_SQR(polyFitSigma);
         // psF32 polyFitNorm = exp(tmpPoly->coedd[0] + PS_SQR(polyFitMean) / (2.0 * PS_SQR(polyFitSigma)));
