lugetP(A)
Returns matrix P of the LU factorization of matrix A. P is permutation matrix such that A = P*L*U. Matrix P is not the identity matrix in the case row exchanges were performed.
See also lugetL lugetU cholesky
Example
A = ( 1, 0, 4 ; 3, 2, 1 ; 3, -2, 1 )
1 0 4
3 2 1
3 -2 1
P = lugetP(A)
1 0 0
0 1 0
0 0 1