163 SUBROUTINE cpot05( UPLO, N, NRHS, A, LDA, B, LDB, X, LDX, XACT,
164 $ LDXACT, FERR, BERR, RESLTS )
172 INTEGER LDA, LDB, LDX, LDXACT, N, NRHS
175 REAL BERR( * ), FERR( * ), RESLTS( * )
176 COMPLEX A( LDA, * ), B( LDB, * ), X( LDX, * ),
184 parameter( zero = 0.0e+0, one = 1.0e+0 )
188 INTEGER I, IMAX, J, K
189 REAL AXBI, DIFF, EPS, ERRBND, OVFL, TMP, UNFL, XNORM
196 EXTERNAL lsame, icamax, slamch
199 INTRINSIC abs, aimag, max, min, real
205 cabs1( zdum ) = abs( real( zdum ) ) + abs( aimag( zdum ) )
211 IF( n.LE.0 .OR. nrhs.LE.0 )
THEN
217 eps = slamch(
'Epsilon' )
218 unfl = slamch(
'Safe minimum' )
220 upper = lsame( uplo,
'U' )
228 imax = icamax( n, x( 1, j ), 1 )
229 xnorm = max( cabs1( x( imax, j ) ), unfl )
232 diff = max( diff, cabs1( x( i, j )-xact( i, j ) ) )
235 IF( xnorm.GT.one )
THEN
237 ELSE IF( diff.LE.ovfl*xnorm )
THEN
245 IF( diff / xnorm.LE.ferr( j ) )
THEN
246 errbnd = max( errbnd, ( diff / xnorm ) / ferr( j ) )
258 tmp = cabs1( b( i, k ) )
261 tmp = tmp + cabs1( a( j, i ) )*cabs1( x( j, k ) )
263 tmp = tmp + abs( real( a( i, i ) ) )*cabs1( x( i, k ) )
265 tmp = tmp + cabs1( a( i, j ) )*cabs1( x( j, k ) )
269 tmp = tmp + cabs1( a( i, j ) )*cabs1( x( j, k ) )
271 tmp = tmp + abs( real( a( i, i ) ) )*cabs1( x( i, k ) )
273 tmp = tmp + cabs1( a( j, i ) )*cabs1( x( j, k ) )
279 axbi = min( axbi, tmp )
282 tmp = berr( k ) / ( ( n+1 )*eps+( n+1 )*unfl /
283 $ max( axbi, ( n+1 )*unfl ) )
287 reslts( 2 ) = max( reslts( 2 ), tmp )
subroutine cpot05(uplo, n, nrhs, a, lda, b, ldb, x, ldx, xact, ldxact, ferr, berr, reslts)
CPOT05