125 SUBROUTINE ctrt01( UPLO, DIAG, N, A, LDA, AINV, LDAINV, RCOND,
135 INTEGER lda, ldainv, n
140 COMPLEX a( lda, * ), ainv( ldainv, * )
147 parameter( zero = 0.0e+0, one = 1.0e+0 )
151 REAL ainvnm, anorm, eps
177 anorm =
clantr(
'1', uplo, diag, n, n, a, lda, rwork )
178 ainvnm =
clantr(
'1', uplo, diag, n, n, ainv, ldainv, rwork )
179 IF( anorm.LE.zero .OR. ainvnm.LE.zero )
THEN
184 rcond = ( one / anorm ) / ainvnm
188 IF(
lsame( diag,
'U' ) )
THEN
196 IF(
lsame( uplo,
'U' ) )
THEN
198 CALL
ctrmv(
'Upper',
'No transpose', diag,
j, a, lda,
203 CALL
ctrmv(
'Lower',
'No transpose', diag, n-
j+1, a(
j,
j ),
204 $ lda, ainv(
j,
j ), 1 )
211 ainv(
j,
j ) = ainv(
j,
j ) - one
216 resid =
clantr(
'1', uplo,
'Non-unit', n, n, ainv, ldainv, rwork )
218 resid = ( ( resid*rcond ) /
REAL( N ) ) / eps
logical function lsame(CA, CB)
LSAME
real function slamch(CMACH)
SLAMCH
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j
subroutine ctrmv(UPLO, TRANS, DIAG, N, A, LDA, X, INCX)
CTRMV
real function clantr(NORM, UPLO, DIAG, M, N, A, LDA, WORK)
CLANTR returns the value of the 1-norm, or the Frobenius norm, or the infinity norm, or the element of largest absolute value of a trapezoidal or triangular matrix.
subroutine ctrt01(UPLO, DIAG, N, A, LDA, AINV, LDAINV, RCOND, RWORK, RESID)
CTRT01