Computes an estimate of the reciprocal condition number of a matrix.
USAGE:
rcond.Matrix(x, one.norm=T)
REQUIRED ARGUMENTS:
x:
A numeric or complex Matrix inheriting from class "Matrix".
OPTIONAL ARGUMENTS:
one.norm:
A logical variable which indicateds computation of the 1-norm condition
estimate if TRUE and the infinity-norm condition estimate otherwise.
VALUE:
A numeric value of class "rcond", representing the reciprocal one or infinity
norm condition estimate.
A copy of the call to "rcond" is returned as an attribute.
DETAILS:
A QR factorization is computed and the condition estimate returned is that of
the triangular factor, using the functions dgeqpf, dtzrqf, dtrcon, zgeqpf,
ztzrqf, ztrcon from Lapack (Anderson et al. (1994)).
The condition number of a square matrix is the product of the norm of that
matrix and the norm of its inverse. Its values fall in the range [1, Inf),
where a value of Inf would imply a singular matrix. A matrix is said to
be ill-conditioned if its has a large condition number. Another way to
view a condition number of a matrix is as a factor by which errors for
solutions to systems of equations with that matrix as coefficient matrix
can be multiplied. Condition numbers usually are estimated rather than
computed exactly for reasons of efficiency.
REFERENCES:
Anderson, E., et al. (1994).
LAPACK User's Guide,
2nd edition, SIAM, Philadelphia.
Golub, G., and Van Loan, C. F. (1989).
Matrix Computations,
2nd edition, Johns Hopkins, Baltimore.