127 SUBROUTINE dpot06( UPLO, N, NRHS, A, LDA, X, LDX, B, LDB,
137 INTEGER lda, ldb, ldx, n, nrhs
138 DOUBLE PRECISION resid
141 DOUBLE PRECISION a( lda, * ),
b( ldb, * ), rwork( * ),
148 DOUBLE PRECISION zero, one, negone
149 parameter( zero = 0.0d+0, one = 1.0d+0 )
150 parameter( negone = -1.0d+0 )
154 DOUBLE PRECISION anorm, bnorm, eps, xnorm
171 IF( n.LE.0 .OR. nrhs.EQ.0 )
THEN
179 anorm =
dlansy(
'I', uplo, n, a, lda, rwork )
180 IF( anorm.LE.zero )
THEN
188 CALL
dsymm(
'Left', uplo, n, nrhs, negone, a, lda, x,
197 xnorm = abs(x(
idamax( n, x( 1,
j ), 1 ),
j))
198 IF( xnorm.LE.zero )
THEN
201 resid = max( resid, ( ( bnorm / anorm ) / xnorm ) / eps )
subroutine dsymm(SIDE, UPLO, M, N, ALPHA, A, LDA, B, LDB, BETA, C, LDC)
DSYMM
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real b(3) integer i
double precision function dlamch(CMACH)
DLAMCH
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j
double precision function dlansy(NORM, UPLO, N, A, LDA, WORK)
DLANSY returns the value of the 1-norm, or the Frobenius norm, or the infinity norm, or the element of largest absolute value of a real symmetric matrix.
integer function idamax(N, DX, INCX)
IDAMAX
subroutine dpot06(UPLO, N, NRHS, A, LDA, X, LDX, B, LDB, RWORK, RESID)
DPOT06