Changeset 15824 for trunk/psLib/src/math/psMatrix.c

Timestamp:

Dec 13, 2007, 4:49:44 PM (18 years ago)

Author:

Paul Price

Message:

LU decomp doesn't allow in-place processing.

File:

• trunk/psLib/src/math/psMatrix.c (modified) (2 diffs)

trunk/psLib/src/math/psMatrix.c

@@ -22 +22 @@
  *  @author Andy Becker, University of Washington (SVD).
  *
- *  @version $Revision: 1.54 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2007-12-14 02:46:45 $
+ *  @version $Revision: 1.55 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2007-12-14 02:49:44 $
  *
  *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -363 +363 @@
         return false;
     }
-    b = psMatrixLUSolve(b, lu, b, perm);
+    psVector *ans = psMatrixLUSolve(NULL, lu, b, perm); // Answer
     psFree(lu);
     psFree(perm);
-    if (!b) {
+    if (!ans) {
         psError(PS_ERR_UNKNOWN, false, "Unable to solve matrix equation.");
         return false;
     }
+
+    memcpy(b->data.U8, ans->data.U8, b->n * PSELEMTYPE_SIZEOF(ans->type.type));
+    psFree(ans);
 
     return true;