126 $ ldc, rwork, resid )
135 INTEGER lda, ldafac, ldc, n
136 DOUBLE PRECISION resid
140 DOUBLE PRECISION rwork( * )
141 COMPLEX*16 a( lda, * ), afac( ldafac, * ), c( ldc, * )
147 DOUBLE PRECISION zero, one
148 parameter( zero = 0.0d+0, one = 1.0d+0 )
149 COMPLEX*16 czero, cone
150 parameter( czero = ( 0.0d+0, 0.0d+0 ),
151 $ cone = ( 1.0d+0, 0.0d+0 ) )
155 DOUBLE PRECISION anorm, eps
166 INTRINSIC dimag, dble
180 anorm =
zlanhe(
'1', uplo, n, a, lda, rwork )
186 IF( dimag( afac(
j,
j ) ).NE.zero )
THEN
194 CALL
zlaset(
'Full', n, n, czero, cone, c, ldc )
198 CALL
zlavhe_rook( uplo,
'Conjugate',
'Non-unit', n, n, afac,
199 $ ldafac, ipiv, c, ldc, info )
203 CALL
zlavhe_rook( uplo,
'No transpose',
'Unit', n, n, afac,
204 $ ldafac, ipiv, c, ldc, info )
208 IF(
lsame( uplo,
'U' ) )
THEN
211 c( i,
j ) = c( i,
j ) - a( i,
j )
213 c(
j,
j ) = c(
j,
j ) - dble( a(
j,
j ) )
217 c(
j,
j ) = c(
j,
j ) - dble( a(
j,
j ) )
219 c( i,
j ) = c( i,
j ) - a( i,
j )
226 resid =
zlanhe(
'1', uplo, n, c, ldc, rwork )
228 IF( anorm.LE.zero )
THEN
232 resid = ( ( resid/dble( n ) )/anorm ) / eps