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