|
D.3.2.15 U_D_O
Procedure from library linalg.lib (see linalg_lib).
- Usage:
- U_D_O(A); constant invertible matrix A
- Return:
- list Z: Z[1]=P , Z[2]=U , Z[3]=D , Z[4]=O
gives a permutation matrix P,
a normalized lower triangular matrix U ,
a diagonal matrix D, and
a normalized upper triangular matrix O
with P*A=U*D*O
- Note:
- Z[1]=-1 means that A is not regular (proc uses gaussred)
Example:
|