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

Timestamp:

May 6, 2009, 2:18:37 PM (17 years ago)

Author:

eugene

Message:

rename psMatrixLUSolve to psMatrixLUSolution, psMatrixLUD to psMatrixLUDecomposition, add psMatrixLUInvert and distinct psMatrixLUSolve vs psMatrixGJSolve

File:

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

trunk/psLib/src/math/psMatrix.h

@@ -39 +39 @@
  *  @return  psImage* : Pointer to LU decomposed psImage.
  */
-psImage* psMatrixLUD(
+psImage *psMatrixLUDecomposition(
     psImage* out,                      ///< Image to return, or NULL.
     psVector** perm,                   ///< Output permutation vector used by psMatrixLUSolve.
@@ -54 +54 @@
  *  @return  psVector* : Pointer to psVector solution of matrix equation.
  */
-psVector* psMatrixLUSolve(
+psVector *psMatrixLUSolution(
     psVector* out,                     ///< Vector to return, or NULL.
     const psImage* LU,                 ///< LU-decomposed matrix.
@@ -61 +61 @@
 );
 
-/** Used to be a Gauss-Jordan numerical solver, but now uses LU decomposition.
+/** LU Decomposition-based Matrix inversion
+ *
+ *  @return  psImage * : Pointer to psImage inverse of matrix
+ */
+psImage *psMatrixLUInvert(
+    psImage *out,                  ///< place result here if not NULL
+    const psImage* LU,             ///< LU-decomposed matrix.
+    const psVector* perm           ///< Permutation vector resulting from psMatrixLUD function.
+);
+
+/** LU Decomposition-based solver for Ax = B.
+ *
+ *  @return bool:   True if successful.
+ */
+bool psMatrixLUSolve(
+    psImage *A,                   ///< Matrix to be solved
+    psVector *b                   ///< Vector of values
+);
+
+/** Gauss-Jordan-based solver for Ax = B.
  *
  *  @return bool:   True if successful.
@@ -69 +88 @@
     psVector *b                   ///< Vector of values
 );
-
-/** Gauss-Jordan numerical solver for F32 input data */
-#define psMatrixGJSolveF32(A,B) psMatrixGJSolve(A,B)
-
 
 /** Invert psImage matrix.