Changeset 12695 for branches/rel-1_0/psLib/src/math/psMatrix.c

Timestamp:

Mar 30, 2007, 10:25:50 AM (19 years ago)

Author:

rhl

Message:

singular matrix test was checking wrong vector

File:

• branches/rel-1_0/psLib/src/math/psMatrix.c (modified) (2 diffs)

branches/rel-1_0/psLib/src/math/psMatrix.c

@@ -22 +22 @@
  *  @author Andy Becker, University of Washington (SVD).
  *
- *  @version $Revision: 1.47.2.2 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2007-03-08 23:21:39 $
+ *  @version $Revision: 1.47.2.3 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2007-03-30 20:25:50 $
  *
  *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -348 +348 @@
     int solver_failed = gsl_linalg_LU_solve(a_gsl, perm_gsl, b_gsl, x_gsl);
 
-    for (int i = 0; i < b->n && !solver_failed; i++) { // gsl_linalg_LU_solve doesn't always catch errors
-        if (!gsl_finite(b_gsl->data[i])) {
+    for (int i = 0; i < x_gsl->size && !solver_failed; i++) { // gsl_linalg_LU_solve doesn't always catch errors
+        if (!gsl_finite(x_gsl->data[i])) {
             solver_failed = 1;
         }