Index: /trunk/psLib/src/dataManip/psMinimize.c
===================================================================
--- /trunk/psLib/src/dataManip/psMinimize.c	(revision 3852)
+++ /trunk/psLib/src/dataManip/psMinimize.c	(revision 3853)
@@ -9,6 +9,6 @@
  *  @author GLG, MHPCC
  *
- *  @version $Revision: 1.111 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2005-05-05 22:10:12 $
+ *  @version $Revision: 1.112 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2005-05-05 22:15:16 $
  *
  *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -808,5 +808,6 @@
 }
 
-# define SWAP(X,Y) {double tmp=(X); (X) = (Y); (Y) = tmp;}
+// XXX" put this in constants.h
+# define PS_SWAP(X,Y) {double tmp=(X); (X) = (Y); (Y) = tmp;}
 
 // XXX EAM : temporary gauss-jordan solver based on gene's
@@ -829,6 +830,7 @@
     indxr = psAlloc (Nx*sizeof(int));
     ipiv  = psAlloc (Nx*sizeof(int));
-    for (j = 0; j < Nx; j++)
+    for (j = 0; j < Nx; j++) {
         ipiv[j] = 0;
+    }
 
     irow = icol = 0;
@@ -864,7 +866,7 @@
         if (irow != icol) {
             for (l = 0; l < Nx; l++) {
-                SWAP (matrix[irow][l], matrix[icol][l]);
-            }
-            SWAP (vector[irow], vector[icol]);
+                PS_SWAP (matrix[irow][l], matrix[icol][l]);
+            }
+            PS_SWAP (vector[irow], vector[icol]);
         }
         indxr[i] = irow;
@@ -886,6 +888,7 @@
                 dum = matrix[ll][icol];
                 matrix[ll][icol] = 0.0;
-                for (l = 0; l < Nx; l++)
+                for (l = 0; l < Nx; l++) {
                     matrix[ll][l] -= matrix[icol][l]*dum;
+                }
                 vector[ll] -= vector[icol]*dum;
             }
@@ -894,7 +897,9 @@
 
     for (l = Nx - 1; l >= 0; l--) {
-        if (indxr[l] != indxc[l])
-            for (k = 0; k < Nx; k++)
-                SWAP (matrix[k][indxr[l]], matrix[k][indxc[l]]);
+        if (indxr[l] != indxc[l]) {
+            for (k = 0; k < Nx; k++) {
+                PS_SWAP (matrix[k][indxr[l]], matrix[k][indxc[l]]);
+            }
+        }
     }
     psFree (ipiv);
Index: /trunk/psLib/src/dataManip/psMinimize.h
===================================================================
--- /trunk/psLib/src/dataManip/psMinimize.h	(revision 3852)
+++ /trunk/psLib/src/dataManip/psMinimize.h	(revision 3853)
@@ -8,6 +8,6 @@
  *  @author GLG, MHPCC
  *
- *  @version $Revision: 1.40 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2005-05-05 22:10:12 $
+ *  @version $Revision: 1.41 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2005-05-05 22:15:16 $
  *
  *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -76,4 +76,5 @@
                                );
 
+// XXX: Add Doxygen comments here.
 typedef
 psF64 (*psMinimizeLMChi2Func)(psVector *deriv,
Index: /trunk/psLib/src/math/psMinimize.c
===================================================================
--- /trunk/psLib/src/math/psMinimize.c	(revision 3852)
+++ /trunk/psLib/src/math/psMinimize.c	(revision 3853)
@@ -9,6 +9,6 @@
  *  @author GLG, MHPCC
  *
- *  @version $Revision: 1.111 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2005-05-05 22:10:12 $
+ *  @version $Revision: 1.112 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2005-05-05 22:15:16 $
  *
  *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -808,5 +808,6 @@
 }
 
-# define SWAP(X,Y) {double tmp=(X); (X) = (Y); (Y) = tmp;}
+// XXX" put this in constants.h
+# define PS_SWAP(X,Y) {double tmp=(X); (X) = (Y); (Y) = tmp;}
 
 // XXX EAM : temporary gauss-jordan solver based on gene's
@@ -829,6 +830,7 @@
     indxr = psAlloc (Nx*sizeof(int));
     ipiv  = psAlloc (Nx*sizeof(int));
-    for (j = 0; j < Nx; j++)
+    for (j = 0; j < Nx; j++) {
         ipiv[j] = 0;
+    }
 
     irow = icol = 0;
@@ -864,7 +866,7 @@
         if (irow != icol) {
             for (l = 0; l < Nx; l++) {
-                SWAP (matrix[irow][l], matrix[icol][l]);
-            }
-            SWAP (vector[irow], vector[icol]);
+                PS_SWAP (matrix[irow][l], matrix[icol][l]);
+            }
+            PS_SWAP (vector[irow], vector[icol]);
         }
         indxr[i] = irow;
@@ -886,6 +888,7 @@
                 dum = matrix[ll][icol];
                 matrix[ll][icol] = 0.0;
-                for (l = 0; l < Nx; l++)
+                for (l = 0; l < Nx; l++) {
                     matrix[ll][l] -= matrix[icol][l]*dum;
+                }
                 vector[ll] -= vector[icol]*dum;
             }
@@ -894,7 +897,9 @@
 
     for (l = Nx - 1; l >= 0; l--) {
-        if (indxr[l] != indxc[l])
-            for (k = 0; k < Nx; k++)
-                SWAP (matrix[k][indxr[l]], matrix[k][indxc[l]]);
+        if (indxr[l] != indxc[l]) {
+            for (k = 0; k < Nx; k++) {
+                PS_SWAP (matrix[k][indxr[l]], matrix[k][indxc[l]]);
+            }
+        }
     }
     psFree (ipiv);
Index: /trunk/psLib/src/math/psMinimize.h
===================================================================
--- /trunk/psLib/src/math/psMinimize.h	(revision 3852)
+++ /trunk/psLib/src/math/psMinimize.h	(revision 3853)
@@ -8,6 +8,6 @@
  *  @author GLG, MHPCC
  *
- *  @version $Revision: 1.40 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2005-05-05 22:10:12 $
+ *  @version $Revision: 1.41 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2005-05-05 22:15:16 $
  *
  *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -76,4 +76,5 @@
                                );
 
+// XXX: Add Doxygen comments here.
 typedef
 psF64 (*psMinimizeLMChi2Func)(psVector *deriv,