Quaternion Function Reference

lu

LU decomposition
(Quaternion and octonion overloadings of standard MATLAB® function)

Syntax

[L,U,P] = lu(A)

Description

The third output parameter is optional. Computes a decomposition such that P * A = L * U, where L and U are lower and upper triangular respectively.

If the third output parameter is omitted, the result in L will be permuted so that A = L * U. This behaviour matches that of the corresponding MATLAB® function.

The LU decomposition is not guaranteed to exist for all quaternion matrices. At present there is no checking for zero divisors or divisors of zero during the Gaussian elimination steps, and consequently the results returned may contain NaNs without warning. The same applies to octonion matrices.

The algorithm used includes partial pivoting and is based on that described by Golub and van Loan for real matrices.

See Also

MATLAB® function: lu
QTFM functions: qr, svd

References

  1. Gene H. Golub and Charles van Loan, Matrix Computations, 3rd edition, Johns Hopkins University Press, 1996.

© 2008-2013 Stephen J. Sangwine and Nicolas Le Bihan

License terms.