Index: /trunk/psLib/src/math/psMinimizePolyFit.c
===================================================================
--- /trunk/psLib/src/math/psMinimizePolyFit.c	(revision 6938)
+++ /trunk/psLib/src/math/psMinimizePolyFit.c	(revision 6939)
@@ -10,6 +10,6 @@
  *  @author EAM, IfA
  *
- *  @version $Revision: 1.9 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2006-03-06 20:57:42 $
+ *  @version $Revision: 1.10 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2006-04-21 20:59:51 $
  *
  *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -38,4 +38,6 @@
 /* DEFINE STATEMENTS                                                         */
 /*****************************************************************************/
+
+# define USE_GAUSS_JORDAN 1
 
 #define PS_VECTOR_GEN_CHEBY_INDEX(VEC, SIZE, TYPE) \
@@ -391,5 +393,5 @@
     }
 
-    if (0) {
+    if (USE_GAUSS_JORDAN) {
         // GaussJordan version
         if (false == psGaussJordan(A, B)) {
@@ -555,5 +557,5 @@
     }
 
-    if (0) {
+    if (USE_GAUSS_JORDAN) {
         // GaussJordan version
         if (false == psGaussJordan(A, B)) {
@@ -846,5 +848,5 @@
 
     // XXX: rel10_ifa used psGaussJordan().  However, this failed tests.  So, I'm using psMatrixLUD().
-    if (0) {
+    if (USE_GAUSS_JORDAN) {
         // GaussJordan version
         if (false == psGaussJordan(A, B)) {
@@ -1766,5 +1768,5 @@
 
     // XXX: rel10_ifa used psGaussJordan().  However, this failed tests.  So, I'm using psMatrixLUD().
-    if (0) {
+    if (USE_GAUSS_JORDAN) {
         // does the solution in place
         // The matrices were overflowing, so I switched to LUD.
@@ -2323,5 +2325,5 @@
 
     // XXX: rel10_ifa used psGaussJordan().  However, this failed tests.  So, I'm using psMatrixLUD().
-    if (0) {
+    if (USE_GAUSS_JORDAN) {
         // does the solution in place
         // The GaussJordan version was overflowing, so I'm using LUD.