Note: If you are working with complex inputs, use the LA_LUDC procedure instead.
See the description of LUSOL for an example using this procedure.
An n by n array of any type except string. Upon output, A is replaced with its LU decomposition.
An output vector that records the row permutations which occurred as a result of partial pivoting.
Set this keyword if the input array A is in column-major format (composed of column vectors) rather than in row-major format (composed of row vectors).
Set this keyword to force the computation to be done in double-precision arithmetic.
An output variable that is set to positive 1 if the number of row interchanges was even, or to negative 1 if the number of interchanges was odd.
Resources and References
LUDC is based on the routine ludcmp described in section 2.3 of Numerical Recipes in C: The Art of Scientific Computing (Second Edition), published by Cambridge University Press, and is used by permission.