1 SUBROUTINE dpot01( UPLO, N, A, LDA, AFAC, LDAFAC, RWORK, RESID )
10 DOUBLE PRECISION resid
13 DOUBLE PRECISION a( lda, * ), afac( ldafac, * ), rwork( * )
61 DOUBLE PRECISION zero, one
62 parameter( zero = 0.0d+0, one = 1.0d+0 )
66 DOUBLE PRECISION anorm, eps, t
70 DOUBLE PRECISION ddot, dlamch, dlansy
71 EXTERNAL lsame, ddot, dlamch, dlansy
74 EXTERNAL dscal, dsyr, dtrmv
90 eps = dlamch(
'Epsilon' )
91 anorm = dlansy(
'1',
uplo, n, a, lda, rwork )
92 IF( anorm.LE.zero )
THEN
99 IF( lsame(
uplo,
'U' ) )
THEN
104 t = ddot( k, afac( 1, k ), 1, afac( 1, k ), 1 )
109 CALL dtrmv(
'Upper',
'Transpose',
'Non-unit', k-1, afac,
110 $ ldafac, afac( 1, k ), 1 )
123 $ CALL dsyr(
'Lower', n-k, one, afac( k+1, k ), 1,
124 $ afac( k+1, k+1 ), ldafac )
129 CALL dscal( n-k+1, t, afac( k, k ), 1 )
136 IF( lsame(
uplo,
'U' ) )
THEN
139 afac( i, j ) = afac( i, j ) - a( i, j )
145 afac( i, j ) = afac( i, j ) - a( i, j )
152 resid = dlansy(
'1',
uplo, n, afac, ldafac, rwork )
154 resid = ( ( resid / dble( n ) ) / anorm ) / eps