Index: trunk/archive/pslib/include/psMatrix.h
===================================================================
--- trunk/archive/pslib/include/psMatrix.h	(revision 753)
+++ trunk/archive/pslib/include/psMatrix.h	(revision 1585)
@@ -15,14 +15,29 @@
 /* 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,17 +55,4 @@
 ;
 
-/** 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,5 +67,5 @@
 psVector *
 psMatrixToVector(psVector *out,		///< Vector to return, or NULL
-		 psImage *in)		///< Matrix to convert
+		 const psImage *in)	///< Matrix to convert
 ;
 
@@ -71,5 +73,5 @@
 psImage *
 psVectorToMatrix(psImage *out,		///< Matrix to return, or NULL
-		 psVector *in)		///< Vector to convert
+		 const psVector *in)	///< Vector to convert
 ;