procedure replaces an
, with the
LU decomposition of a row-wise permutation of itself.
LUDC is based on the routine
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.
array of any type except string. Upon output,
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
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.
See the description of LUSOL for an example using this function.