38 DOUBLE PRECISION anrm, ccond, rcond
42 INTEGER hl( 2 ), hpiv( 2 )
43 DOUBLE PRECISION r( nmax ), r1( nmax ), r2( nmax )
44 COMPLEX*16 a( nmax, nmax ), af( nmax, nmax ), b( nmax ),
45 $ w( 2*nmax ), x( nmax )
53 $ zgbtrf, zgbtrs, zgecon,
zgeequ, zgerfs, zgetf2,
54 $ zgetrf, zgetri, zgetrs
62 common / infoc / infot, nout, ok, lerr
63 common / srnamc / srnamt
66 INTRINSIC dble, dcmplx
71 WRITE( nout, fmt = * )
83 a( i, j ) = dcmplx( 1.d0 / dble( i+j ),
84 $ -1.d0 / dble( i+j ) )
85 af( i, j ) = dcmplx( 1.d0 / dble( i+j ),
86 $ -1.d0 / dble( i+j ) )
100 IF(
lsamen( 2, c2,
'GE' ) )
THEN
108 $ 2, 1, hl, hpiv, info )
115 CALL
chkxer(
'ZGETRF', infot, nout, info, ok )
118 CALL
chkxer(
'ZGETRF', infot, nout, info, ok )
121 CALL
chkxer(
'ZGETRF', infot, nout, info, ok )
129 CALL
chkxer(
'ZGETRS', infot, nout, info, ok )
133 CALL
chkxer(
'ZGETRS', infot, nout, info, ok )
137 CALL
chkxer(
'ZGETRS', infot, nout, info, ok )
141 CALL
chkxer(
'ZGETRS', infot, nout, info, ok )
145 CALL
chkxer(
'ZGETRS', infot, nout, info, ok )
149 CALL plasma_dealloc_handle( hl, info )
150 CALL plasma_dealloc_handle( hpiv, info )
158 CALL
chkxer(
'ZGETRF', infot, nout, info, ok )
161 CALL
chkxer(
'ZGETRF', infot, nout, info, ok )
164 CALL
chkxer(
'ZGETRF', infot, nout, info, ok )
172 CALL
chkxer(
'ZGETRS', infot, nout, info, ok )
176 CALL
chkxer(
'ZGETRS', infot, nout, info, ok )
180 CALL
chkxer(
'ZGETRS', infot, nout, info, ok )
184 CALL
chkxer(
'ZGETRS', infot, nout, info, ok )
188 CALL
chkxer(
'ZGETRS', infot, nout, info, ok )
194 CALL
alaesm( path, ok, nout )