124 SUBROUTINE dsyt01( UPLO, N, A, LDA, AFAC, LDAFAC, IPIV, C, LDC,
134 INTEGER lda, ldafac, ldc, n
135 DOUBLE PRECISION resid
139 DOUBLE PRECISION a( lda, * ), afac( ldafac, * ), c( ldc, * ),
146 DOUBLE PRECISION zero, one
147 parameter( zero = 0.0d+0, one = 1.0d+0 )
151 DOUBLE PRECISION anorm, eps
176 anorm =
dlansy(
'1', uplo, n, a, lda, rwork )
180 CALL
dlaset(
'Full', n, n, zero, one, c, ldc )
184 CALL
dlavsy( uplo,
'Transpose',
'Non-unit', n, n, afac, ldafac,
185 $ ipiv, c, ldc, info )
189 CALL
dlavsy( uplo,
'No transpose',
'Unit', n, n, afac, ldafac,
190 $ ipiv, c, ldc, info )
194 IF(
lsame( uplo,
'U' ) )
THEN
197 c( i,
j ) = c( i,
j ) - a( i,
j )
203 c( i,
j ) = c( i,
j ) - a( i,
j )
210 resid =
dlansy(
'1', uplo, n, c, ldc, rwork )
212 IF( anorm.LE.zero )
THEN
216 resid = ( ( resid / dble( n ) ) / anorm ) / eps
LOGICAL function lsame(CA, CB)
LSAME
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.
subroutine dlavsy(UPLO, TRANS, DIAG, N, NRHS, A, LDA, IPIV, B, LDB, INFO)
DLAVSY
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j
DOUBLE PRECISION function dlamch(CMACH)
DLAMCH
subroutine dlaset(UPLO, M, N, ALPHA, BETA, A, LDA)
DLASET initializes the off-diagonal elements and the diagonal elements of a matrix to given values...
subroutine dsyt01(UPLO, N, A, LDA, AFAC, LDAFAC, IPIV, C, LDC, RWORK, RESID)
DSYT01