SCULIB_INVERT_MATRIX
invert a square matrix
Description:
This routine inverts a matrix by
a method opaque to casual inspection but which has the advantage that it doesn’t require any more
space than the input matrix itself provides. The method is described in ‘Data Reduction and Error
Analysis for the Physical Sciences’ by Bevington and Robinson, and the code is adapted from the
Pascal version listed there.
Invocation
CALL SCULIB_INVERT_MATRIX (M, ARRAY, DET, IK, JK,
STATUS)
Arguments
M = INTEGER (Given)
the dimensions of the matrix
ARRAY (M,M) =
DOUBLE PRECISION (Given and returned)
the matrix to be inverted
DET = DOUBLE
PRECISION (Returned)
the determinant of the matrix
IK (M) = INTEGER (Scratch)
i
indices of largest array elements
JK (M) = INTEGER (Scratch)
j indices of largest array
elements
STATUS = INTEGER (Given and returned)
global status
Copyright
Copyright
©1995,1996,1997,1998,1999 Particle Physics and Astronomy Research Council. All Rights Reserved.