121 SUBROUTINE cspt02( UPLO, N, NRHS, A, X, LDX, B, LDB, RWORK,
130 INTEGER LDB, LDX, N, NRHS
135 COMPLEX A( * ), B( LDB, * ), X( LDX, * )
142 parameter( zero = 0.0e+0, one = 1.0e+0 )
144 parameter( cone = ( 1.0e+0, 0.0e+0 ) )
148 REAL ANORM, BNORM, EPS, XNORM
151 REAL CLANSP, SCASUM, SLAMCH
152 EXTERNAL clansp, scasum, slamch
164 IF( n.LE.0 .OR. nrhs.LE.0 )
THEN
171 eps = slamch(
'Epsilon' )
172 anorm = clansp(
'1', uplo, n, a, rwork )
173 IF( anorm.LE.zero )
THEN
181 CALL cspmv( uplo, n, -cone, a, x( 1, j ), 1, cone, b( 1, j ),
190 bnorm = scasum( n, b( 1, j ), 1 )
191 xnorm = scasum( n, x( 1, j ), 1 )
192 IF( xnorm.LE.zero )
THEN
195 resid = max( resid, ( ( bnorm/anorm )/xnorm )/eps )
subroutine cspt02(uplo, n, nrhs, a, x, ldx, b, ldb, rwork, resid)
CSPT02
subroutine cspmv(uplo, n, alpha, ap, x, incx, beta, y, incy)
CSPMV computes a matrix-vector product for complex vectors using a complex symmetric packed matrix