126 $ ldc, rwork, resid )
135 INTEGER lda, ldafac, ldc, n
141 COMPLEX a( lda, * ), afac( ldafac, * ), c( ldc, * )
148 parameter( zero = 0.0e+0, one = 1.0e+0 )
150 parameter( czero = ( 0.0e+0, 0.0e+0 ),
151 $ cone = ( 1.0e+0, 0.0e+0 ) )
180 anorm =
clansy(
'1', uplo, n, a, lda, rwork )
184 CALL
claset(
'Full', n, n, czero, cone, c, ldc )
188 CALL
clavsy_rook( uplo,
'Transpose',
'Non-unit', n, n, afac,
189 $ ldafac, ipiv, c, ldc, info )
193 CALL
clavsy_rook( uplo,
'No transpose',
'Unit', n, n, afac,
194 $ ldafac, ipiv, c, ldc, info )
198 IF(
lsame( uplo,
'U' ) )
THEN
201 c( i,
j ) = c( i,
j ) - a( i,
j )
207 c( i,
j ) = c( i,
j ) - a( i,
j )
214 resid =
clansy(
'1', uplo, n, c, ldc, rwork )
216 IF( anorm.LE.zero )
THEN
220 resid = ( ( resid/
REAL( N ) )/anorm ) / eps