Index: trunk/psLib/src/math/psBinaryOp.h
===================================================================
--- trunk/psLib/src/math/psBinaryOp.h	(revision 4160)
+++ trunk/psLib/src/math/psBinaryOp.h	(revision 4162)
@@ -1,3 +1,2 @@
-
 /** @file  psBinaryOp.h
  *
@@ -31,6 +30,6 @@
  *  @author Robert DeSonia, MHPCC
  *
- *  @version $Revision: 1.1 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2005-02-10 02:36:41 $
+ *  @version $Revision: 1.2 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2005-06-08 23:40:45 $
  *
  *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -65,3 +64,3 @@
 /// @}
 
-#endif
+#endif // #ifndef PSBINARY_OP_H
Index: trunk/psLib/src/math/psCompare.h
===================================================================
--- trunk/psLib/src/math/psCompare.h	(revision 4160)
+++ trunk/psLib/src/math/psCompare.h	(revision 4162)
@@ -6,11 +6,11 @@
  *  @ingroup Compare
  *
- *  @version $Revision: 1.6 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2005-02-17 19:26:23 $
+ *  @version $Revision: 1.7 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2005-06-08 23:40:45 $
  *
  *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
  */
 
-#if !defined(PS_COMPARE_H)
+#ifndef PS_COMPARE_H
 #define PS_COMPARE_H
 
@@ -24,7 +24,7 @@
  *  e.g., for psList of pointers to numeric values.
  *
- *  @return int      an integer less than, equal to, or greater than zero if 
- *                   the first argument is considered to be respectively less 
- *                   than, equal to, or greater than the second. 
+ *  @return int      an integer less than, equal to, or greater than zero if
+ *                   the first argument is considered to be respectively less
+ *                   than, equal to, or greater than the second.
  */
 typedef int (*psComparePtrFcn) (
@@ -35,7 +35,7 @@
 /** A comparison function for sorting.
  *
- *  @return int      an integer less than, equal to, or greater than zero if 
- *                   the first argument is considered to be respectively less 
- *                   than, equal to, or greater than the second. 
+ *  @return int      an integer less than, equal to, or greater than zero if
+ *                   the first argument is considered to be respectively less
+ *                   than, equal to, or greater than the second.
  */
 typedef int (*psCompareFcn) (
@@ -46,7 +46,7 @@
 /** Compare function of psS8 data.  For use with psListSort.
  *
- *  @return int      an integer less than, equal to, or greater than zero if 
- *                   the first argument is considered to be respectively less 
- *                   than, equal to, or greater than the second. 
+ *  @return int      an integer less than, equal to, or greater than zero if
+ *                   the first argument is considered to be respectively less
+ *                   than, equal to, or greater than the second.
  */
 int psCompareS8Ptr(
@@ -57,7 +57,7 @@
 /** Compare function of psS16 data.  For use with psListSort.
  *
- *  @return int      an integer less than, equal to, or greater than zero if 
- *                   the first argument is considered to be respectively less 
- *                   than, equal to, or greater than the second. 
+ *  @return int      an integer less than, equal to, or greater than zero if
+ *                   the first argument is considered to be respectively less
+ *                   than, equal to, or greater than the second.
  */
 int psCompareS16Ptr(
@@ -68,7 +68,7 @@
 /** Compare function of psS32 data.  For use with psListSort.
  *
- *  @return int      an integer less than, equal to, or greater than zero if 
- *                   the first argument is considered to be respectively less 
- *                   than, equal to, or greater than the second. 
+ *  @return int      an integer less than, equal to, or greater than zero if
+ *                   the first argument is considered to be respectively less
+ *                   than, equal to, or greater than the second.
  */
 int psCompareS32Ptr(
@@ -79,7 +79,7 @@
 /** Compare function of psS64 data.  For use with psListSort.
  *
- *  @return int      an integer less than, equal to, or greater than zero if 
- *                   the first argument is considered to be respectively less 
- *                   than, equal to, or greater than the second. 
+ *  @return int      an integer less than, equal to, or greater than zero if
+ *                   the first argument is considered to be respectively less
+ *                   than, equal to, or greater than the second.
  */
 int psCompareS64Ptr(
@@ -90,7 +90,7 @@
 /** Compare function of psU8 data.  For use with psListSort.
  *
- *  @return int      an integer less than, equal to, or greater than zero if 
- *                   the first argument is considered to be respectively less 
- *                   than, equal to, or greater than the second. 
+ *  @return int      an integer less than, equal to, or greater than zero if
+ *                   the first argument is considered to be respectively less
+ *                   than, equal to, or greater than the second.
  */
 int psCompareU8Ptr(
@@ -101,7 +101,7 @@
 /** Compare function of psU16 data.  For use with psListSort.
  *
- *  @return int      an integer less than, equal to, or greater than zero if 
- *                   the first argument is considered to be respectively less 
- *                   than, equal to, or greater than the second. 
+ *  @return int      an integer less than, equal to, or greater than zero if
+ *                   the first argument is considered to be respectively less
+ *                   than, equal to, or greater than the second.
  */
 int psCompareU16Ptr(
@@ -112,7 +112,7 @@
 /** Compare function of psU32 data.  For use with psListSort.
  *
- *  @return int      an integer less than, equal to, or greater than zero if 
- *                   the first argument is considered to be respectively less 
- *                   than, equal to, or greater than the second. 
+ *  @return int      an integer less than, equal to, or greater than zero if
+ *                   the first argument is considered to be respectively less
+ *                   than, equal to, or greater than the second.
  */
 int psCompareU32Ptr(
@@ -123,7 +123,7 @@
 /** Compare function of psU64 data.  For use with psListSort.
  *
- *  @return int      an integer less than, equal to, or greater than zero if 
- *                   the first argument is considered to be respectively less 
- *                   than, equal to, or greater than the second. 
+ *  @return int      an integer less than, equal to, or greater than zero if
+ *                   the first argument is considered to be respectively less
+ *                   than, equal to, or greater than the second.
  */
 int psCompareU64Ptr(
@@ -134,7 +134,7 @@
 /** Compare function of psF32 data.  For use with psListSort.
  *
- *  @return int      an integer less than, equal to, or greater than zero if 
- *                   the first argument is considered to be respectively less 
- *                   than, equal to, or greater than the second. 
+ *  @return int      an integer less than, equal to, or greater than zero if
+ *                   the first argument is considered to be respectively less
+ *                   than, equal to, or greater than the second.
  */
 int psCompareF32Ptr(
@@ -145,7 +145,7 @@
 /** Compare function of psF64 data.  For use with psListSort.
  *
- *  @return int      an integer less than, equal to, or greater than zero if 
- *                   the first argument is considered to be respectively less 
- *                   than, equal to, or greater than the second. 
+ *  @return int      an integer less than, equal to, or greater than zero if
+ *                   the first argument is considered to be respectively less
+ *                   than, equal to, or greater than the second.
  */
 int psCompareF64Ptr(
@@ -156,7 +156,7 @@
 /** Compare function of psS8 data.  For use with psListSort for descending ordering.
  *
- *  @return int      an integer less than, equal to, or greater than zero if 
- *                   the first argument is considered to be respectively greater 
- *                   than, equal to, or less than the second. 
+ *  @return int      an integer less than, equal to, or greater than zero if
+ *                   the first argument is considered to be respectively greater
+ *                   than, equal to, or less than the second.
  */
 int psCompareDescendingS8Ptr(
@@ -167,7 +167,7 @@
 /** Compare function of psS16 data.  For use with psListSort for descending ordering.
  *
- *  @return int      an integer less than, equal to, or greater than zero if 
- *                   the first argument is considered to be respectively greater 
- *                   than, equal to, or less than the second. 
+ *  @return int      an integer less than, equal to, or greater than zero if
+ *                   the first argument is considered to be respectively greater
+ *                   than, equal to, or less than the second.
  */
 int psCompareDescendingS16Ptr(
@@ -178,7 +178,7 @@
 /** Compare function of psS32 data.  For use with psListSort for descending ordering.
  *
- *  @return int      an integer less than, equal to, or greater than zero if 
- *                   the first argument is considered to be respectively greater 
- *                   than, equal to, or less than the second. 
+ *  @return int      an integer less than, equal to, or greater than zero if
+ *                   the first argument is considered to be respectively greater
+ *                   than, equal to, or less than the second.
  */
 int psCompareDescendingS32Ptr(
@@ -189,7 +189,7 @@
 /** Compare function of psS64 data.  For use with psListSort for descending ordering.
  *
- *  @return int      an integer less than, equal to, or greater than zero if 
- *                   the first argument is considered to be respectively greater 
- *                   than, equal to, or less than the second. 
+ *  @return int      an integer less than, equal to, or greater than zero if
+ *                   the first argument is considered to be respectively greater
+ *                   than, equal to, or less than the second.
  */
 int psCompareDescendingS64Ptr(
@@ -200,7 +200,7 @@
 /** Compare function of psU8 data.  For use with psListSort for descending ordering.
  *
- *  @return int      an integer less than, equal to, or greater than zero if 
- *                   the first argument is considered to be respectively greater 
- *                   than, equal to, or less than the second. 
+ *  @return int      an integer less than, equal to, or greater than zero if
+ *                   the first argument is considered to be respectively greater
+ *                   than, equal to, or less than the second.
  */
 int psCompareDescendingU8Ptr(
@@ -211,7 +211,7 @@
 /** Compare function of psU16 data.  For use with psListSort for descending ordering.
  *
- *  @return int      an integer less than, equal to, or greater than zero if 
- *                   the first argument is considered to be respectively greater 
- *                   than, equal to, or less than the second. 
+ *  @return int      an integer less than, equal to, or greater than zero if
+ *                   the first argument is considered to be respectively greater
+ *                   than, equal to, or less than the second.
  */
 int psCompareDescendingU16Ptr(
@@ -222,7 +222,7 @@
 /** Compare function of psU32 data.  For use with psListSort for descending ordering.
  *
- *  @return int      an integer less than, equal to, or greater than zero if 
- *                   the first argument is considered to be respectively greater 
- *                   than, equal to, or lessg than the second. 
+ *  @return int      an integer less than, equal to, or greater than zero if
+ *                   the first argument is considered to be respectively greater
+ *                   than, equal to, or lessg than the second.
  */
 int psCompareDescendingU32Ptr(
@@ -233,7 +233,7 @@
 /** Compare function of psU64 data.  For use with psListSort for descending ordering.
  *
- *  @return int      an integer less than, equal to, or greater than zero if 
- *                   the first argument is considered to be respectively greater 
- *                   than, equal to, or lessg than the second. 
+ *  @return int      an integer less than, equal to, or greater than zero if
+ *                   the first argument is considered to be respectively greater
+ *                   than, equal to, or lessg than the second.
  */
 int psCompareDescendingU64Ptr(
@@ -244,7 +244,7 @@
 /** Compare function of psF32 data.  For use with psListSort for descending ordering.
  *
- *  @return int      an integer less than, equal to, or greater than zero if 
- *                   the first argument is considered to be respectively greater 
- *                   than, equal to, or lessg than the second. 
+ *  @return int      an integer less than, equal to, or greater than zero if
+ *                   the first argument is considered to be respectively greater
+ *                   than, equal to, or lessg than the second.
  */
 int psCompareDescendingF32Ptr(
@@ -255,7 +255,7 @@
 /** Compare function of psF64 data.  For use with psListSort for descending ordering.
  *
- *  @return int      an integer less than, equal to, or greater than zero if 
- *                   the first argument is considered to be respectively greater 
- *                   than, equal to, or lessg than the second. 
+ *  @return int      an integer less than, equal to, or greater than zero if
+ *                   the first argument is considered to be respectively greater
+ *                   than, equal to, or lessg than the second.
  */
 int psCompareDescendingF64Ptr(
@@ -266,7 +266,7 @@
 /** Compare function of psS8 data.
  *
- *  @return int      an integer less than, equal to, or greater than zero if 
- *                   the first argument is considered to be respectively less 
- *                   than, equal to, or greater than the second. 
+ *  @return int      an integer less than, equal to, or greater than zero if
+ *                   the first argument is considered to be respectively less
+ *                   than, equal to, or greater than the second.
  */
 int psCompareS8(
@@ -278,7 +278,7 @@
 /** Compare function of psS16 data.
  *
- *  @return int      an integer less than, equal to, or greater than zero if 
- *                   the first argument is considered to be respectively less 
- *                   than, equal to, or greater than the second. 
+ *  @return int      an integer less than, equal to, or greater than zero if
+ *                   the first argument is considered to be respectively less
+ *                   than, equal to, or greater than the second.
  */
 int psCompareS16(
@@ -290,7 +290,7 @@
 /** Compare function of psS32 data.
  *
- *  @return int      an integer less than, equal to, or greater than zero if 
- *                   the first argument is considered to be respectively less 
- *                   than, equal to, or greater than the second. 
+ *  @return int      an integer less than, equal to, or greater than zero if
+ *                   the first argument is considered to be respectively less
+ *                   than, equal to, or greater than the second.
  */
 int psCompareS32(
@@ -302,7 +302,7 @@
 /** Compare function of psS64 data.
  *
- *  @return int      an integer less than, equal to, or greater than zero if 
- *                   the first argument is considered to be respectively less 
- *                   than, equal to, or greater than the second. 
+ *  @return int      an integer less than, equal to, or greater than zero if
+ *                   the first argument is considered to be respectively less
+ *                   than, equal to, or greater than the second.
  */
 int psCompareS64(
@@ -314,7 +314,7 @@
 /** Compare function of psU8 data.
  *
- *  @return int      an integer less than, equal to, or greater than zero if 
- *                   the first argument is considered to be respectively less 
- *                   than, equal to, or greater than the second. 
+ *  @return int      an integer less than, equal to, or greater than zero if
+ *                   the first argument is considered to be respectively less
+ *                   than, equal to, or greater than the second.
  */
 int psCompareU8(
@@ -326,7 +326,7 @@
 /** Compare function of psU16 data.
  *
- *  @return int      an integer less than, equal to, or greater than zero if 
- *                   the first argument is considered to be respectively less 
- *                   than, equal to, or greater than the second. 
+ *  @return int      an integer less than, equal to, or greater than zero if
+ *                   the first argument is considered to be respectively less
+ *                   than, equal to, or greater than the second.
  */
 int psCompareU16(
@@ -338,7 +338,7 @@
 /** Compare function of psU32 data.
  *
- *  @return int      an integer less than, equal to, or greater than zero if 
- *                   the first argument is considered to be respectively less 
- *                   than, equal to, or greater than the second. 
+ *  @return int      an integer less than, equal to, or greater than zero if
+ *                   the first argument is considered to be respectively less
+ *                   than, equal to, or greater than the second.
  */
 int psCompareU32(
@@ -350,7 +350,7 @@
 /** Compare function of psU64 data.
  *
- *  @return int      an integer less than, equal to, or greater than zero if 
- *                   the first argument is considered to be respectively less 
- *                   than, equal to, or greater than the second. 
+ *  @return int      an integer less than, equal to, or greater than zero if
+ *                   the first argument is considered to be respectively less
+ *                   than, equal to, or greater than the second.
  */
 int psCompareU64(
@@ -362,7 +362,7 @@
 /** Compare function of psF32 data.
  *
- *  @return int      an integer less than, equal to, or greater than zero if 
- *                   the first argument is considered to be respectively less 
- *                   than, equal to, or greater than the second. 
+ *  @return int      an integer less than, equal to, or greater than zero if
+ *                   the first argument is considered to be respectively less
+ *                   than, equal to, or greater than the second.
  */
 int psCompareF32(
@@ -374,7 +374,7 @@
 /** Compare function of psF64 data.
  *
- *  @return int      an integer less than, equal to, or greater than zero if 
- *                   the first argument is considered to be respectively less 
- *                   than, equal to, or greater than the second. 
+ *  @return int      an integer less than, equal to, or greater than zero if
+ *                   the first argument is considered to be respectively less
+ *                   than, equal to, or greater than the second.
  */
 int psCompareF64(
@@ -386,7 +386,7 @@
 /** Compare function of psS8 data.
  *
- *  @return int      an integer less than, equal to, or greater than zero if 
- *                   the first argument is considered to be respectively greater 
- *                   than, equal to, or less than the second. 
+ *  @return int      an integer less than, equal to, or greater than zero if
+ *                   the first argument is considered to be respectively greater
+ *                   than, equal to, or less than the second.
  */
 int psCompareDescendingS8(
@@ -398,7 +398,7 @@
 /** Compare function of psS16 data.
  *
- *  @return int      an integer less than, equal to, or greater than zero if 
- *                   the first argument is considered to be respectively greater 
- *                   than, equal to, or less than the second. 
+ *  @return int      an integer less than, equal to, or greater than zero if
+ *                   the first argument is considered to be respectively greater
+ *                   than, equal to, or less than the second.
  */
 int psCompareDescendingS16(
@@ -410,7 +410,7 @@
 /** Compare function of psS32 data.
  *
- *  @return int      an integer less than, equal to, or greater than zero if 
- *                   the first argument is considered to be respectively greater 
- *                   than, equal to, or less than the second. 
+ *  @return int      an integer less than, equal to, or greater than zero if
+ *                   the first argument is considered to be respectively greater
+ *                   than, equal to, or less than the second.
  */
 int psCompareDescendingS32(
@@ -422,7 +422,7 @@
 /** Compare function of psS64 data.
  *
- *  @return int      an integer less than, equal to, or greater than zero if 
- *                   the first argument is considered to be respectively greater 
- *                   than, equal to, or less than the second. 
+ *  @return int      an integer less than, equal to, or greater than zero if
+ *                   the first argument is considered to be respectively greater
+ *                   than, equal to, or less than the second.
  */
 int psCompareDescendingS64(
@@ -434,7 +434,7 @@
 /** Compare function of psU8 data.
  *
- *  @return int      an integer less than, equal to, or greater than zero if 
- *                   the first argument is considered to be respectively greater 
- *                   than, equal to, or less than the second. 
+ *  @return int      an integer less than, equal to, or greater than zero if
+ *                   the first argument is considered to be respectively greater
+ *                   than, equal to, or less than the second.
  */
 int psCompareDescendingU8(
@@ -446,7 +446,7 @@
 /** Compare function of psU16 data.
  *
- *  @return int      an integer less than, equal to, or greater than zero if 
- *                   the first argument is considered to be respectively greater 
- *                   than, equal to, or less than the second. 
+ *  @return int      an integer less than, equal to, or greater than zero if
+ *                   the first argument is considered to be respectively greater
+ *                   than, equal to, or less than the second.
  */
 int psCompareDescendingU16(
@@ -458,7 +458,7 @@
 /** Compare function of psU32 data.
  *
- *  @return int      an integer less than, equal to, or greater than zero if 
- *                   the first argument is considered to be respectively greater 
- *                   than, equal to, or less than the second. 
+ *  @return int      an integer less than, equal to, or greater than zero if
+ *                   the first argument is considered to be respectively greater
+ *                   than, equal to, or less than the second.
  */
 int psCompareDescendingU32(
@@ -470,7 +470,7 @@
 /** Compare function of psU64 data.
  *
- *  @return int      an integer less than, equal to, or greater than zero if 
- *                   the first argument is considered to be respectively greater 
- *                   than, equal to, or less than the second. 
+ *  @return int      an integer less than, equal to, or greater than zero if
+ *                   the first argument is considered to be respectively greater
+ *                   than, equal to, or less than the second.
  */
 int psCompareDescendingU64(
@@ -482,7 +482,7 @@
 /** Compare function of psF32 data.
  *
- *  @return int      an integer less than, equal to, or greater than zero if 
- *                   the first argument is considered to be respectively greater 
- *                   than, equal to, or less than the second. 
+ *  @return int      an integer less than, equal to, or greater than zero if
+ *                   the first argument is considered to be respectively greater
+ *                   than, equal to, or less than the second.
  */
 int psCompareDescendingF32(
@@ -494,7 +494,7 @@
 /** Compare function of psF64 data.
  *
- *  @return int      an integer less than, equal to, or greater than zero if 
- *                   the first argument is considered to be respectively greater 
- *                   than, equal to, or less than the second. 
+ *  @return int      an integer less than, equal to, or greater than zero if
+ *                   the first argument is considered to be respectively greater
+ *                   than, equal to, or less than the second.
  */
 int psCompareDescendingF64(
@@ -506,3 +506,3 @@
 /// @}
 
-#endif
+#endif  // #ifndef PS_COMPARE_H
Index: trunk/psLib/src/math/psConstants.h
===================================================================
--- trunk/psLib/src/math/psConstants.h	(revision 4160)
+++ trunk/psLib/src/math/psConstants.h	(revision 4162)
@@ -6,6 +6,6 @@
  *  @author GLG, MHPCC
  *
- *  @version $Revision: 1.72 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2005-06-04 00:49:48 $
+ *  @version $Revision: 1.73 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2005-06-08 23:40:45 $
  *
  *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -45,6 +45,5 @@
 #define M_1_PI 0.3183098861837906715377675267450287  /* 1/pi */
 #define M_2_PI 0.6366197723675813430755350534900574  /* 2/pi */
-#endif
-#define PS_PI   3.1415926535897932384626433832795029  /* pi */
+#endif // #ifndef M_PI
 
 #define DEG_TO_RAD(DEGREES) ((DEGREES) * M_PI / 180.0)
@@ -699,15 +698,9 @@
  *****************************************************************************/
 #define PS_MAX(A, B) \
-(((A) > (B)) ? (A) : (B)) \
+(((A) > (B)) ? (A) : (B))
 
 #define PS_MIN(A, B) \
-(((A) < (B)) ? (A) : (B)) \
+(((A) < (B)) ? (A) : (B))
 
 #define PS_SQR(A) \
-((A) * (A)) \
-
-#ifdef DARWIN
-#define PS_SQRT_F32(A) ((float) sqrt(A))
-#else
-#define PS_SQRT_F32(A) (sqrtf(A))
-#endif
+((A) * (A))
Index: trunk/psLib/src/math/psMatrix.h
===================================================================
--- trunk/psLib/src/math/psMatrix.h	(revision 4160)
+++ trunk/psLib/src/math/psMatrix.h	(revision 4162)
@@ -1,3 +1,2 @@
-
 /** @file  psMatrix.h
  *
@@ -22,6 +21,6 @@
  *  @author Ross Harman, MHPCC
  *
- *  @version $Revision: 1.15 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2005-02-17 19:26:23 $
+ *  @version $Revision: 1.16 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2005-06-08 23:40:45 $
  *
  *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -164,3 +163,3 @@
 /// @}
 
-#endif
+#endif // #ifndef PSMATRIX_H
Index: trunk/psLib/src/math/psMinimize.h
===================================================================
--- trunk/psLib/src/math/psMinimize.h	(revision 4160)
+++ trunk/psLib/src/math/psMinimize.h	(revision 4162)
@@ -8,6 +8,6 @@
  *  @author GLG, MHPCC
  *
- *  @version $Revision: 1.42 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2005-05-05 22:23:29 $
+ *  @version $Revision: 1.43 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2005-06-08 23:40:45 $
  *
  *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -15,5 +15,5 @@
  */
 
-#if !defined(PS_MINIMIZE_H)
+#ifndef PS_MINIMIZE_H
 #define PS_MINIMIZE_H
 
@@ -53,11 +53,11 @@
 /** Derive a polynomial fit.
  *
- *  psVectorFitPolynomial1d returns the polynomial that best fits the 
- *  observations. The input parameters are a polynomial that specifies the 
- *  fit order, myPoly, which will be altered and returned with the best-fit 
- *  coefficients; and the observations, x, y and yErr. The independent 
- *  variable list, x may be NULL, in which case the vector index is used. 
- *  The dependent variable error, yErr may be null, in which case the solution 
- *  is determined in the assumption that all data errors are equal. This 
+ *  psVectorFitPolynomial1d returns the polynomial that best fits the
+ *  observations. The input parameters are a polynomial that specifies the
+ *  fit order, myPoly, which will be altered and returned with the best-fit
+ *  coefficients; and the observations, x, y and yErr. The independent
+ *  variable list, x may be NULL, in which case the vector index is used.
+ *  The dependent variable error, yErr may be null, in which case the solution
+ *  is determined in the assumption that all data errors are equal. This
  *  function must be valid only for types psF32, psF64.
  *
@@ -142,4 +142,4 @@
 /* \} */// End of MathGroup Functions
 
-#endif
+#endif // #ifndef PS_MINIMIZE_H
 
Index: trunk/psLib/src/math/psPolynomial.c
===================================================================
--- trunk/psLib/src/math/psPolynomial.c	(revision 4160)
+++ trunk/psLib/src/math/psPolynomial.c	(revision 4162)
@@ -7,6 +7,6 @@
  *  polynomials.  It also contains a Gaussian functions.
  *
- *  @version $Revision: 1.104 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2005-05-25 20:26:55 $
+ *  @version $Revision: 1.105 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2005-06-08 23:40:45 $
  *
  *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -1024,5 +1024,5 @@
 
     if (normal == true) {
-        tmp = 1.0 / PS_SQRT_F32(2.0 * M_PI * (sigma * sigma));
+        tmp = 1.0 / sqrtf(2.0 * M_PI * (sigma * sigma));
     }
 
Index: trunk/psLib/src/math/psPolynomial.h
===================================================================
--- trunk/psLib/src/math/psPolynomial.h	(revision 4160)
+++ trunk/psLib/src/math/psPolynomial.h	(revision 4162)
@@ -12,11 +12,11 @@
  *  @author GLG, MHPCC
  *
- *  @version $Revision: 1.44 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2005-03-31 23:01:46 $
+ *  @version $Revision: 1.45 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2005-06-08 23:40:45 $
  *
  *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
  */
 
-#if !defined(PS_FUNCTIONS_H)
+#ifndef PS_FUNCTIONS_H
 #define PS_FUNCTIONS_H
 
@@ -34,8 +34,8 @@
 
 /** Evaluate a non-normalized Gaussian with the given mean and sigma at the
- *  given coordianate.  
- *
- *  Note that this is not a Gaussian deviate.  The evaluated Gaussian is: 
- *        \f[ exp(-\frac{(x-mean)^2}{2\sigma^2}) \f] 
+ *  given coordianate.
+ *
+ *  Note that this is not a Gaussian deviate.  The evaluated Gaussian is:
+ *        \f[ exp(-\frac{(x-mean)^2}{2\sigma^2}) \f]
  *
  *  @return psF32      value on the gaussian curve given the input parameters
@@ -49,8 +49,8 @@
 
 /** Produce a vector of random numbers from a Gaussian distribution with
- *  the specified mean and sigma 
- *  
+ *  the specified mean and sigma
+ *
  *  @return psVector*    vector of random numbers
- *  
+ *
  */
 psVector* p_psGaussianDev(
@@ -433,8 +433,8 @@
                          const psVector* y,
                          psF32 X);
-#endif
-
-/* \} */// End of MathGroup Functions
-
-#endif
-
+#endif // #if 0
+
+/** \} */ // End of MathGroup Functions
+
+#endif // #ifndef PS_FUNCTIONS_H
+
Index: trunk/psLib/src/math/psRandom.h
===================================================================
--- trunk/psLib/src/math/psRandom.h	(revision 4160)
+++ trunk/psLib/src/math/psRandom.h	(revision 4162)
@@ -10,11 +10,11 @@
 *  @author GLG, MHPCC
 *
-*  @version $Revision: 1.3 $ $Name: not supported by cvs2svn $
-*  @date $Date: 2005-02-17 19:26:23 $
+*  @version $Revision: 1.4 $ $Name: not supported by cvs2svn $
+*  @date $Date: 2005-06-08 23:40:45 $
 *
 *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
 */
 
-#if !defined(PS_RANDOM_H)
+#ifndef PS_RANDOM_H
 #define PS_RANDOM_H
 
@@ -56,4 +56,3 @@
 /* \} */// End of MathGroup Functions
 
-#endif
-
+#endif // #ifndef PS_RANDOM_H
Index: trunk/psLib/src/math/psSpline.c
===================================================================
--- trunk/psLib/src/math/psSpline.c	(revision 4160)
+++ trunk/psLib/src/math/psSpline.c	(revision 4162)
@@ -7,6 +7,6 @@
  *  polynomials.  It also contains a Gaussian functions.
  *
- *  @version $Revision: 1.104 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2005-05-25 20:26:55 $
+ *  @version $Revision: 1.105 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2005-06-08 23:40:45 $
  *
  *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -1024,5 +1024,5 @@
 
     if (normal == true) {
-        tmp = 1.0 / PS_SQRT_F32(2.0 * M_PI * (sigma * sigma));
+        tmp = 1.0 / sqrtf(2.0 * M_PI * (sigma * sigma));
     }
 
Index: trunk/psLib/src/math/psSpline.h
===================================================================
--- trunk/psLib/src/math/psSpline.h	(revision 4160)
+++ trunk/psLib/src/math/psSpline.h	(revision 4162)
@@ -12,11 +12,11 @@
  *  @author GLG, MHPCC
  *
- *  @version $Revision: 1.44 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2005-03-31 23:01:46 $
+ *  @version $Revision: 1.45 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2005-06-08 23:40:45 $
  *
  *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
  */
 
-#if !defined(PS_FUNCTIONS_H)
+#ifndef PS_FUNCTIONS_H
 #define PS_FUNCTIONS_H
 
@@ -34,8 +34,8 @@
 
 /** Evaluate a non-normalized Gaussian with the given mean and sigma at the
- *  given coordianate.  
- *
- *  Note that this is not a Gaussian deviate.  The evaluated Gaussian is: 
- *        \f[ exp(-\frac{(x-mean)^2}{2\sigma^2}) \f] 
+ *  given coordianate.
+ *
+ *  Note that this is not a Gaussian deviate.  The evaluated Gaussian is:
+ *        \f[ exp(-\frac{(x-mean)^2}{2\sigma^2}) \f]
  *
  *  @return psF32      value on the gaussian curve given the input parameters
@@ -49,8 +49,8 @@
 
 /** Produce a vector of random numbers from a Gaussian distribution with
- *  the specified mean and sigma 
- *  
+ *  the specified mean and sigma
+ *
  *  @return psVector*    vector of random numbers
- *  
+ *
  */
 psVector* p_psGaussianDev(
@@ -433,8 +433,8 @@
                          const psVector* y,
                          psF32 X);
-#endif
-
-/* \} */// End of MathGroup Functions
-
-#endif
-
+#endif // #if 0
+
+/** \} */ // End of MathGroup Functions
+
+#endif // #ifndef PS_FUNCTIONS_H
+
Index: trunk/psLib/src/math/psStats.c
===================================================================
--- trunk/psLib/src/math/psStats.c	(revision 4160)
+++ trunk/psLib/src/math/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)));
Index: trunk/psLib/src/math/psStats.h
===================================================================
--- trunk/psLib/src/math/psStats.h	(revision 4160)
+++ trunk/psLib/src/math/psStats.h	(revision 4162)
@@ -14,10 +14,10 @@
  *  @author GLG, MHPCC
  *
- *  @version $Revision: 1.41 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2005-06-07 23:11:54 $
+ *  @version $Revision: 1.42 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2005-06-08 23:40:45 $
  *
  *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
  */
-#if !defined(PS_STATS_H)
+#ifndef PS_STATS_H
 #define PS_STATS_H
 
@@ -112,9 +112,9 @@
 /** The basic histogram structure which contains bounds and bins.
  *
- *  In this structure, the vector bounds specifies the boundaries of the 
- *  histogram bins, and must of type psF32, while nums specifies the number 
- *  of entries in the bin, and must of type psU32. The value of bounds.n must 
- *  therefore be 1 greater than than nums.n. The two values minNum and maxNum 
- *  are the number of data values which fell below the lower limit bound or 
+ *  In this structure, the vector bounds specifies the boundaries of the
+ *  histogram bins, and must of type psF32, while nums specifies the number
+ *  of entries in the bin, and must of type psU32. The value of bounds.n must
+ *  therefore be 1 greater than than nums.n. The two values minNum and maxNum
+ *  are the number of data values which fell below the lower limit bound or
  *  above the upper limit bound, respectively.
  */
@@ -130,6 +130,6 @@
 
 /** Allocator for psHistogram where the bounds of the bins are implicitly
- *  specified through simply specifying an upper and lower limit along with 
- *  the size of the bins. 
+ *  specified through simply specifying an upper and lower limit along with
+ *  the size of the bins.
  *
  *  @return psHistogram*    Newly allocated psHistogram
@@ -142,5 +142,5 @@
 
 /** Allocator for psHistogram where the bounds of the bins are explicitly
- *  specified. 
+ *  specified.
  *
  *  @return psHistogram*    Newly allocated psHistogram
@@ -152,5 +152,5 @@
 /** Calculate a histogram
  *
- *  The following function populates the histogram bins from the specified 
+ *  The following function populates the histogram bins from the specified
  *  vector (in). It alters and returns the histogram out structure. The input
  *  vector may be of types psU8, psU16, psF32, psF64.
@@ -169,5 +169,5 @@
  *
  *  @return psBool    If more than one statistic result is set in stats->options,
- *                    false is returned and the value parameter is not set, 
+ *                    false is returned and the value parameter is not set,
  *                    otherwise true is returned.
  */
@@ -193,3 +193,3 @@
 /// @}
 
-#endif
+#endif // #ifndef PS_STATS_H
Index: trunk/psLib/src/math/psUnaryOp.h
===================================================================
--- trunk/psLib/src/math/psUnaryOp.h	(revision 4160)
+++ trunk/psLib/src/math/psUnaryOp.h	(revision 4162)
@@ -30,6 +30,6 @@
  *  @author Robert DeSonia, MHPCC
  *
- *  @version $Revision: 1.1 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2005-02-10 02:36:41 $
+ *  @version $Revision: 1.2 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2005-06-08 23:40:45 $
  *
  *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -67,3 +67,3 @@
 /// @}
 
-#endif
+#endif // #ifndef PSUNARY_OP_H