125 $ ldc, rwork, resid )
134 INTEGER lda, ldafac, ldc, n
139 REAL a( lda, * ), afac( ldafac, * ), c( ldc, * ),
147 parameter( zero = 0.0e+0, one = 1.0e+0 )
176 anorm =
slansy(
'1', uplo, n, a, lda, rwork )
180 CALL
slaset(
'Full', n, n, zero, one, c, ldc )
184 CALL
slavsy_rook( uplo,
'Transpose',
'Non-unit', n, n, afac,
185 $ ldafac, ipiv, c, ldc, info )
189 CALL
slavsy_rook( uplo,
'No transpose',
'Unit', n, n, afac,
190 $ ldafac, 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 =
slansy(
'1', uplo, n, c, ldc, rwork )
212 IF( anorm.LE.zero )
THEN
216 resid = ( ( resid /
REAL( N ) ) / anorm ) / eps