Changeset 1585 for trunk/archive/pslib/include/psMatrix.h

Timestamp:

Aug 18, 2004, 5:18:31 PM (22 years ago)

Author:

Paul Price

Message:

Updating from SDRS.

File:

• trunk/archive/pslib/include/psMatrix.h (modified) (4 diffs)

trunk/archive/pslib/include/psMatrix.h

@@ -15 +15 @@
 /* Linear Algebra */
 
+/** LU Decomposition of a matrix */
+psImage *
+psMatrixLUD(psImage *out,               ///< Matrix to return, or NULL
+            psVector *perm,             ///< Permutation vector
+            const psImage *in)          ///< Matrix to decompose
+;
+
+/** LU Solution.  Solves for and returns x in the equation Ax = b */
+psVector *
+psMatrixLUSolve(psVector *out,          ///< Vector to return, or NULL
+                const psImage *LU,      ///< LU-decomposed matrix
+                const psVector *RHS,    ///< right-hand-side of the equation
+                const psVector *perm    ///< Permutation vector
+;
+
 /** Invert matrix.  Not using restrict, to allow inversion to be done in-place */
 psImage *
 psMatrixInvert(psImage *out,            ///< Matrix to return, or NULL
                const psImage *myMatrix, ///< Matrix to be inverted
-               float *restrict determinant) ///< Determinant to return, or NULL
+               float *determinant)      ///< Determinant to return, or NULL
 ;
 
 /** Matrix determinant */
 float
-psMatrixDeterminant(const psImage *restrict myMatrix) ///< Matrix to get determinant for
+psMatrixDeterminant(const psImage *myMatrix) ///< Matrix to get determinant for
 ;
 
@@ -40 +55 @@
 ;
 
-/** LU Decomposition of a matrix */
-psImage *
-psMatrixLUD(psImage *out,               ///< Matrix to return, or NULL
-            psImage *in)                ///< Matrix to decompose
-;
-
-/** LU Solution.  Solves for and returns x in the equation Ax = b */
-psVector *
-psMatrixLUSolve(psVector *out,          ///< Vector to return, or NULL
-                const psImage *LU,      ///< LU-decomposed matrix
-                const psVector *RHS)    ///< right-hand-side of the equation
-;
-
 /** Eigenvectors of a matrix */
 psVector *
@@ -65 +67 @@
 psVector *
 psMatrixToVector(psVector *out,         ///< Vector to return, or NULL
-                 psImage *in)           ///< Matrix to convert
+                 const psImage *in)     ///< Matrix to convert
 ;
 
@@ -71 +73 @@
 psImage *
 psVectorToMatrix(psImage *out,          ///< Matrix to return, or NULL
-                 psVector *in)          ///< Vector to convert
+                 const psVector *in)    ///< Vector to convert
 ;