Computes the infinity norm condition number of op(A)*diag(x) for general matrices.
call cla_gercond_x( trans, n, a, lda, af, ldaf, ipiv, x, info, work, rwork )
call zla_gercond_x( trans, n, a, lda, af, ldaf, ipiv, x, info, work, rwork )
The function computes the infinity norm condition number of
op(A) * diag(x)
where the x is a COMPLEX vector for cla_gercond_x and a DOUBLE COMPLEX vector for zla_gercond_x.
CHARACTER*1. Must be 'N' or 'T' or 'C'.
Specifies the form of the system of equations:
If trans = 'N', the system has the form A*X = B (No transpose)
If trans = 'T', the system has the form AT*X = B (Transpose)
If trans = 'C', the system has the form AH*X = B (Conjugate Transpose = Transpose)
INTEGER. The number of linear equations, that is, the order of the matrix A; n≥ 0.
- a, af, x, work
COMPLEX for cla_gercond_x
DOUBLE COMPLEX for zla_gercond_x
a(lda,*) contains the original general n-by-n matrix A.
af(ldaf,*) contains the factors L and U from the factorization A=P*L*U as returned by ?getrf.
x, DIMENSIONn. The vector x in the formula op(A) * diag(x).
work is a workspace array of DIMENSION (2*n).
The second dimension of a and af must be at least max(1, n).
INTEGER. The leading dimension of the array a. lda≥max(1,n).
INTEGER. The leading dimension of af. ldaf≥max(1,n).
Array with DIMENSIONn. The pivot indices from the factorization A = P*L*U as computed by ?getrf. Row i of the matrix was interchanged with row ipiv(i).
REAL for cla_gercond_x
DOUBLE PRECISION for zla_gercond_x
Array rwork with DIMENSIONn is a workspace.
If info = 0, the execution is successful.
If i > 0, the i-th parameter is invalid.