Changeset 15820 for trunk/psLib/src/math/psMatrix.h

Timestamp:

Dec 13, 2007, 2:41:17 PM (18 years ago)

Author:

Paul Price

Message:

Removing Numerical Recipes code that was noted to be 'temporary': replaced the Gauss-Jordan solver with our LU Decomposition solver, since I couldn't find a GPL Gauss-Jordan solver code (none in GSL). NR says that the LU solver operation count is a factor of 3 better than GJ (p48 of NR ed 2), but we could be slowed down by converting to GSL vector/matrix representation; this could be sped up by doing something other than copying the elements individually that MHPCC originally implemented.

File:

• trunk/psLib/src/math/psMatrix.h (modified) (3 diffs)

trunk/psLib/src/math/psMatrix.h

@@ -19 +19 @@
  * @author Ross Harman, MHPCC
  *
- * @version $Revision: 1.27 $ $Name: not supported by cvs2svn $
- * @date $Date: 2007-03-07 20:48:55 $
+ * @version $Revision: 1.28 $ $Name: not supported by cvs2svn $
+ * @date $Date: 2007-12-14 00:41:17 $
  * Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
  */
@@ -61 +61 @@
 );
 
-/** Gauss-Jordan numerical solver.
+/** Used to be a Gauss-Jordan numerical solver, but now uses LU decomposition.
  *
  *  @return bool:   True if successful.
@@ -70 +70 @@
 );
 
-/** Gauss-Jordan numerical solver for F32 input data
- *
- *  @return bool:   True if successful.
- */
-bool psMatrixGJSolveF32(
-    psImage *A,                   ///< Matrix to be solved
-    psVector *b                   ///< Vector of values
-);
+/** Gauss-Jordan numerical solver for F32 input data */
+#define psMatrixGJSolveF32(A,B) psMatrixGJSolve(A,B)
+
 
 /** Invert psImage matrix.