LUDecomposition D

jet.spad line 6224 [edit on github]

• D: Field

LUDecomposition contains procedures to solve linear systems of equations or to compute inverses using a LU decomposition.

LUDecomp: Matrix D -> Record(LU: Matrix D, Perm: Vector Integer, Pivots: List D)

LUDecomp(A) computes a LU decomposition of A using the algorithm of Crout. LU contains both triangular matrices; Perm is the permutation used for partial pivoting and Pivots yields the used pivots.

LUInverse: Matrix D -> Record(Inv: Matrix D, Pivots: List D)

LUInverse(A) computes the inverse of A using a LU decomposition.

LUSolve: (Matrix D, Vector Integer, Vector D) -> Vector D

LUSolve(LU, Perm, B) uses a previously computed LU decomposition to solve a linear system with right hand side B. LU and Perm are as given by LUDecomp.