compute the PLU or PLU1R decomposition of a Matrix
LUDecomposition(A, out, options)
(optional) equation of the form output = obj where obj is one of 'P', 'L', 'U', 'U1', or 'R', or a list consisting of one or more of these names; selects result objects to compute
(optional) parameters; for a complete list, see LinearAlgebra[LUDecomposition]
The LUDecomposition command computes a PLU decomposition or a modified PLU1R decomposition of the Matrix A.
Depending on what is included in the output option (out), an expression sequence containing one or more of the factors P, L, U, U1, and R can be returned. The objects are returned in the same order as specified in the output list.
Note: Either U or the pair U1 and R may be returned, but not both.
The LUDecomposition(A) calling sequence is equivalent to LUDecomposition(A, output=['P','L','U']). This LU decomposition generates a square unit lower triangular L factor and an upper triangular factor U with the same dimensions as A so that A=P·L·U. The Matrix P is a permutation Matrix.
The PLU1R decomposition is achieved by using LUDecomposition(A, output=['P','L','U1','R']). This further factors U into U1·R where U1 is square upper triangular factor and R is the unique reduced row echelon form of the Matrix A. In this case, A=P·L·U1·R.
The output option (out) determines the content of the returned expression sequence.
A ≔ 0,−2,0,3|1,3,0,1|1,1,0,0|−3,4,1,0
p,l,u ≔ LUDecomposition⁡A
B ≔ 1,0,2|3,1,4|6,1,4|1,1,2:
To reduce a Matrix using Gaussian elimination, specify the 'U' object:
To reduce a Matrix using Gauss-Jordan elimination, specify the 'R' object:
For more information on the P*L*U1*R decomposition see:
Corless, Robert M., and Jeffrey, David J. "The Turing Factorization of a Rectangular Matrix." Sigsam Bulletin, Vol. 31 no. 3. (September 1997): 20-28.
This paper names the U1 factor U.
Download Help Document