66 INTEGER I, INFO, ITER, J
70 DOUBLE PRECISION A( nmax, nmax ), AF( nmax, nmax ), B( nmax ),
71 $ c( nmax ), r( nmax ), r1( nmax ), r2( nmax ),
72 $ w( 2*nmax ), x( nmax )
73 DOUBLE PRECISION WORK(1)
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 )
116 CALL
dsgesv(-1,0,a,1,ip,b,1,x,1,work,swork,iter,info)
117 CALL
chkxer(
'DSGESV', infot, nout, lerr, ok )
119 CALL
dsgesv(0,-1,a,1,ip,b,1,x,1,work,swork,iter,info)
120 CALL
chkxer(
'DSGESV', infot, nout, lerr, ok )
122 CALL
dsgesv(2,1,a,1,ip,b,2,x,2,work,swork,iter,info)
123 CALL
chkxer(
'DSGESV', infot, nout, lerr, ok )
125 CALL
dsgesv(2,1,a,2,ip,b,1,x,2,work,swork,iter,info)
126 CALL
chkxer(
'DSGESV', infot, nout, lerr, ok )
128 CALL
dsgesv(2,1,a,2,ip,b,2,x,1,work,swork,iter,info)
129 CALL
chkxer(
'DSGESV', infot, nout, lerr, ok )
134 WRITE( nout, fmt = 9999 )
'DSGESV'
136 WRITE( nout, fmt = 9998 )
'DSGESV'
139 9999
FORMAT( 1x, a6,
' drivers passed the tests of the error exits' )
140 9998
FORMAT(
' *** ', a6,
' drivers failed the tests of the error ',
subroutine chkxer(SRNAMT, INFOT, NOUT, LERR, OK)
subroutine derrab(NUNIT)
DERRAB
subroutine dsgesv(N, NRHS, A, LDA, IPIV, B, LDB, X, LDX, WORK, SWORK, ITER, INFO)
DSGESV computes the solution to system of linear equations A * X = B for GE matrices (mixed precisio...