66 INTEGER I, INFO, ITER, J
70 COMPLEX*16 A( nmax, nmax ), AF( nmax, nmax ), B( nmax ),
71 $ c( nmax ), r( nmax ), r1( nmax ), r2( nmax ),
72 $ w( 2*nmax ), x( nmax )
75 DOUBLE PRECISION RWORK(1)
88 COMMON / infoc / infot, nout, ok, lerr
89 COMMON / srnamc / srnamt
97 WRITE( nout, fmt = * )
103 a( i, j ) = 1.d0 / dble( i+j )
104 af( i, j ) = 1.d0 / dble( i+j )
119 CALL
zcgesv(-1,0,a,1,ip,b,1,x,1,work,swork,rwork,iter,info)
120 CALL
chkxer(
'ZCGESV', infot, nout, lerr, ok )
122 CALL
zcgesv(0,-1,a,1,ip,b,1,x,1,work,swork,rwork,iter,info)
123 CALL
chkxer(
'ZCGESV', infot, nout, lerr, ok )
125 CALL
zcgesv(2,1,a,1,ip,b,2,x,2,work,swork,rwork,iter,info)
126 CALL
chkxer(
'ZCGESV', infot, nout, lerr, ok )
128 CALL
zcgesv(2,1,a,2,ip,b,1,x,2,work,swork,rwork,iter,info)
129 CALL
chkxer(
'ZCGESV', infot, nout, lerr, ok )
131 CALL
zcgesv(2,1,a,2,ip,b,2,x,1,work,swork,rwork,iter,info)
132 CALL
chkxer(
'ZCGESV', infot, nout, lerr, ok )
137 WRITE( nout, fmt = 9999 )
'ZCGESV'
139 WRITE( nout, fmt = 9998 )
'ZCGESV'
142 9999
FORMAT( 1x, a6,
' drivers passed the tests of the error exits' )
143 9998
FORMAT(
' *** ', a6,
' drivers failed the tests of the error ',
subroutine zerrab(NUNIT)
ZERRAB
subroutine chkxer(SRNAMT, INFOT, NOUT, LERR, OK)
subroutine zcgesv(N, NRHS, A, LDA, IPIV, B, LDB, X, LDX, WORK, SWORK, RWORK, ITER, INFO)
ZCGESV computes the solution to system of linear equations A * X = B for GE matrices (mixed precisio...