Code covered by the BSD License

### Highlights from LU factorization with complete pivoting.

3.0
3.0 | 1 rating Rate this file 4 Downloads (last 30 days) File Size: 3.88 KB File ID: #27249 Version: 1.3

# LU factorization with complete pivoting.

### Nick Henderson (view profile)

14 Apr 2010 (Updated )

An implementation of algorithm 3.4.2 from Matrix Computations.

File Information
Description

To compute the LU factorization under default settings:

[L U p q] = lucp(A)

This produces a factorization such that L*U = A(p,q). Vectors p and q permute the rows and columns, respectively.

The pivot tolerance can be controlled:

[L U p q] = lucp(A,tol)

The algorithm will terminate if the absolute value of the pivot is less than tol.

Permutation matrices can be generated:

[L U P Q] = lucp(A,tol,'matrix')
[L U P Q] = lucp(A,tol,'sparse')

The first will generate full permutation matrices P and Q such that L*U = P*A*Q. The second generates sparse P and Q.

If A is sparse, L and U will be sparse. However, no effort is taken to reduce fill in.

This function works on non-square matrices.

Acknowledgements

Gauss Elimination With Complete Pivoting and Gaussian Elimination Using Complete Pivoting inspired this file.

MATLAB release MATLAB 7.8 (R2009a)
25 Apr 2011 cidem m