Index: trunk/psLib/src/math/psMatrix.h
===================================================================
--- trunk/psLib/src/math/psMatrix.h	(revision 15820)
+++ trunk/psLib/src/math/psMatrix.h	(revision 24084)
@@ -39,5 +39,5 @@
  *  @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,5 +54,5 @@
  *  @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,5 +61,24 @@
 );
 
-/** 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,8 +88,4 @@
     psVector *b                   ///< Vector of values
 );
-
-/** Gauss-Jordan numerical solver for F32 input data */
-#define psMatrixGJSolveF32(A,B) psMatrixGJSolve(A,B)
-
 
 /** Invert psImage matrix.