133 SUBROUTINE zget02( TRANS, M, N, NRHS, A, LDA, X, LDX, B, LDB,
143 INTEGER lda, ldb, ldx, m, n, nrhs
144 DOUBLE PRECISION resid
147 DOUBLE PRECISION rwork( * )
148 COMPLEX*16 a( lda, * ),
b( ldb, * ),
x( ldx, * )
154 DOUBLE PRECISION zero, one
155 parameter( zero = 0.0d+0, one = 1.0d+0 )
157 parameter( cone = 1.0d+0 )
161 DOUBLE PRECISION anorm, bnorm, eps, xnorm
178 IF( m.LE.0 .OR. n.LE.0 .OR. nrhs.EQ.0 )
THEN
183 IF(
lsame( trans,
'T' ) .OR.
lsame( trans,
'C' ) )
THEN
194 anorm =
zlange(
'1', n1, n2, a, lda, rwork )
195 IF( anorm.LE.zero )
THEN
202 CALL
zgemm( trans,
'No transpose', n1, nrhs, n2, -cone, a, lda,
x,
203 $ ldx, cone,
b, ldb )
212 IF( xnorm.LE.zero )
THEN
215 resid = max( resid, ( ( bnorm / anorm ) / xnorm ) / eps )
LOGICAL function lsame(CA, CB)
LSAME
DOUBLE PRECISION function dzasum(N, ZX, INCX)
DZASUM
DOUBLE PRECISION function zlange(NORM, M, N, A, LDA, WORK)
ZLANGE returns the value of the 1-norm, Frobenius norm, infinity-norm, or the largest absolute value ...
subroutine zgemm(TRANSA, TRANSB, M, N, K, ALPHA, A, LDA, B, LDB, BETA, C, LDC)
ZGEMM
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real b(3) integer i
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j
DOUBLE PRECISION function dlamch(CMACH)
DLAMCH
subroutine zget02(TRANS, M, N, NRHS, A, LDA, X, LDX, B, LDB, RWORK, RESID)
ZGET02