66 INTEGER I, INFO, ITER, J
69 COMPLEX*16 A( nmax, nmax ), AF( nmax, nmax ), B( nmax ),
70 $ c( nmax ), r( nmax ), r1( nmax ), r2( nmax ),
71 $ w( 2*nmax ), x( nmax )
72 DOUBLE PRECISION RWORK( nmax )
73 COMPLEX*16 WORK(nmax*nmax)
74 COMPLEX SWORK(nmax*nmax)
85 COMMON / infoc / infot, nout, ok, lerr
86 COMMON / srnamc / srnamt
94 WRITE( nout, fmt = * )
100 a( i, j ) = 1.d0 / dble( i+j )
101 af( i, j ) = 1.d0 / dble( i+j )
115 CALL
zcposv(
'/',0,0,a,1,b,1,x,1,work,swork,rwork,iter,info)
116 CALL
chkxer(
'ZCPOSV', infot, nout, lerr, ok )
118 CALL
zcposv(
'U',-1,0,a,1,b,1,x,1,work,swork,rwork,iter,info)
119 CALL
chkxer(
'ZCPOSV', infot, nout, lerr, ok )
121 CALL
zcposv(
'U',0,-1,a,1,b,1,x,1,work,swork,rwork,iter,info)
122 CALL
chkxer(
'ZCPOSV', infot, nout, lerr, ok )
124 CALL
zcposv(
'U',2,1,a,1,b,2,x,2,work,swork,rwork,iter,info)
125 CALL
chkxer(
'ZCPOSV', infot, nout, lerr, ok )
127 CALL
zcposv(
'U',2,1,a,2,b,1,x,2,work,swork,rwork,iter,info)
128 CALL
chkxer(
'ZCPOSV', infot, nout, lerr, ok )
130 CALL
zcposv(
'U',2,1,a,2,b,2,x,1,work,swork,rwork,iter,info)
131 CALL
chkxer(
'ZCPOSV', infot, nout, lerr, ok )
136 WRITE( nout, fmt = 9999 )
'ZCPOSV'
138 WRITE( nout, fmt = 9998 )
'ZCPOSV'
141 9999
FORMAT( 1x, a6,
' drivers passed the tests of the error exits' )
142 9998
FORMAT(
' *** ', a6,
' drivers failed the tests of the error ',
subroutine chkxer(SRNAMT, INFOT, NOUT, LERR, OK)
subroutine zerrac(NUNIT)
ZERRAC
subroutine zcposv(UPLO, N, NRHS, A, LDA, B, LDB, X, LDX, WORK, SWORK, RWORK, ITER, INFO)
ZCPOSV computes the solution to system of linear equations A * X = B for PO matrices ...