Changeset 6939

Timestamp:

Apr 21, 2006, 10:59:51 AM (20 years ago)

Author:

magnier

Message:

restore the GaussJordan inversion

File:

• trunk/psLib/src/math/psMinimizePolyFit.c (modified) (7 diffs)

trunk/psLib/src/math/psMinimizePolyFit.c

@@ -10 +10 @@
  *  @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 +38 @@
 /* DEFINE STATEMENTS                                                         */
 /*****************************************************************************/
+
+# define USE_GAUSS_JORDAN 1
 
 #define PS_VECTOR_GEN_CHEBY_INDEX(VEC, SIZE, TYPE) \
@@ -391 +393 @@
     }
 
-    if (0) {
+    if (USE_GAUSS_JORDAN) {
         // GaussJordan version
         if (false == psGaussJordan(A, B)) {
@@ -555 +557 @@
     }
 
-    if (0) {
+    if (USE_GAUSS_JORDAN) {
         // GaussJordan version
         if (false == psGaussJordan(A, B)) {
@@ -846 +848 @@
 
     // 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 +1768 @@
 
     // 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 +2325 @@
 
     // 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.