108 SUBROUTINE cget03( N, A, LDA, AINV, LDAINV, WORK, LDWORK, RWORK,
116 INTEGER LDA, LDAINV, LDWORK, N
121 COMPLEX A( LDA, * ), AINV( LDAINV, * ),
129 parameter( zero = 0.0e+0, one = 1.0e+0 )
131 parameter( czero = ( 0.0e+0, 0.0e+0 ),
132 $ cone = ( 1.0e+0, 0.0e+0 ) )
136 REAL AINVNM, ANORM, EPS
140 EXTERNAL clange, slamch
160 eps = slamch(
'Epsilon' )
161 anorm = clange(
'1', n, n, a, lda, rwork )
162 ainvnm = clange(
'1', n, n, ainv, ldainv, rwork )
163 IF( anorm.LE.zero .OR. ainvnm.LE.zero )
THEN
168 rcond = ( one/anorm ) / ainvnm
172 CALL cgemm(
'No transpose',
'No transpose', n, n, n, -cone,
173 $ ainv, ldainv, a, lda, czero, work, ldwork )
175 work( i, i ) = cone + work( i, i )
180 resid = clange(
'1', n, n, work, ldwork, rwork )
182 resid = ( ( resid*rcond )/eps ) / real( n )
subroutine cget03(n, a, lda, ainv, ldainv, work, ldwork, rwork, rcond, resid)
CGET03
subroutine cgemm(transa, transb, m, n, k, alpha, a, lda, b, ldb, beta, c, ldc)
CGEMM