132 SUBROUTINE zget02( TRANS, M, N, NRHS, A, LDA, X, LDX, B, LDB,
141 INTEGER LDA, LDB, LDX, M, N, NRHS
142 DOUBLE PRECISION RESID
145 DOUBLE PRECISION RWORK( * )
146 COMPLEX*16 A( LDA, * ), B( LDB, * ), X( LDX, * )
152 DOUBLE PRECISION ZERO, ONE
153 parameter( zero = 0.0d+0, one = 1.0d+0 )
155 parameter( cone = ( 1.0d+0, 0.0d+0 ) )
159 DOUBLE PRECISION ANORM, BNORM, EPS, XNORM
163 DOUBLE PRECISION DLAMCH, DZASUM, ZLANGE
164 EXTERNAL lsame, dlamch, dzasum, zlange
176 IF( m.LE.0 .OR. n.LE.0 .OR. nrhs.EQ.0 )
THEN
181 IF( lsame( trans,
'T' ) .OR. lsame( trans,
'C' ) )
THEN
191 eps = dlamch(
'Epsilon' )
192 IF( lsame( trans,
'N' ) )
THEN
193 anorm = zlange(
'1', m, n, a, lda, rwork )
195 anorm = zlange(
'I', m, n, a, lda, rwork )
197 IF( anorm.LE.zero )
THEN
204 CALL zgemm( trans,
'No transpose', n1, nrhs, n2, -cone, a, lda, x,
205 $ ldx, cone, b, ldb )
212 bnorm = dzasum( n1, b( 1, j ), 1 )
213 xnorm = dzasum( n2, x( 1, j ), 1 )
214 IF( xnorm.LE.zero )
THEN
217 resid = max( resid, ( ( bnorm / anorm ) / xnorm ) / eps )
subroutine zgemm(TRANSA, TRANSB, M, N, K, ALPHA, A, LDA, B, LDB, BETA, C, LDC)
ZGEMM
subroutine zget02(TRANS, M, N, NRHS, A, LDA, X, LDX, B, LDB, RWORK, RESID)
ZGET02