compute the PLU or PLU1R decomposition of a Matrix - Maple Help

Online Help

All Products    Maple    MapleSim


Home : Support : Online Help : Education : Student Package : Linear Algebra : Computation : Solvers : Student/LinearAlgebra/LUDecomposition

Student[LinearAlgebra][LUDecomposition] - compute the PLU or PLU1R decomposition of a Matrix

Calling Sequence

LUDecomposition(A, out, options)

Parameters

A

-

Matrix

out

-

(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

options

-

(optional) parameters; for a complete list, see LinearAlgebra[LUDecomposition]

Description

• 

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.

Examples

withStudent[LinearAlgebra]:

A:=0,2,0,3|1,3,0,1|1,1,0,0|3,4,1,0

A:=0113231400013100

(1)

p,l,u:=LUDecompositionA

p,l,u:=0100100000010010,1000010032112100001,231401130044520001

(2)

p.l.u

0113231400013100

(3)

B:=1,0,2|3,1,4|6,1,4|1,1,2:

To reduce a Matrix using Gaussian elimination, specify the 'U' object:

LUDecompositionA,output='U'

231401130044520001

(4)

To reduce a Matrix using Gauss-Jordan elimination, specify the 'R' object:

LUDecompositionB,output='R'

10010104300113

(5)

See Also

LinearAlgebra[LUDecomposition], Student[LinearAlgebra], Student[LinearAlgebra][GaussianElimination], Student[LinearAlgebra][ReducedRowEchelonForm]

References

  

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

Was this information helpful?



Please add your Comment (Optional)
E-mail Address (Optional)
What is ? This question helps us to combat spam