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)
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real b(3) integer i
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 ...
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j