77 DOUBLE PRECISION ANORM, RCOND
81 DOUBLE PRECISION D( NMAX ), DF( NMAX ), R1( NMAX ), R2( NMAX ),
83 COMPLEX*16 B( NMAX ), DL( NMAX ), DLF( NMAX ), DU( NMAX ),
84 $ DU2( NMAX ), DUF( NMAX ), E( NMAX ),
85 $ EF( NMAX ), W( NMAX ), X( NMAX )
101 COMMON / infoc / infot, nout, ok, lerr
102 COMMON / srnamc / srnamt
107 WRITE( nout, fmt = * )
118 IF(
lsamen( 2, c2,
'GT' ) )
THEN
126 CALL zgttrf( -1, dl, e, du, du2, ip, info )
127 CALL chkxer(
'ZGTTRF', infot, nout, lerr, ok )
133 CALL zgttrs(
'/', 0, 0, dl, e, du, du2, ip, x, 1, info )
134 CALL chkxer(
'ZGTTRS', infot, nout, lerr, ok )
136 CALL zgttrs(
'N', -1, 0, dl, e, du, du2, ip, x, 1, info )
137 CALL chkxer(
'ZGTTRS', infot, nout, lerr, ok )
139 CALL zgttrs(
'N', 0, -1, dl, e, du, du2, ip, x, 1, info )
140 CALL chkxer(
'ZGTTRS', infot, nout, lerr, ok )
142 CALL zgttrs(
'N', 2, 1, dl, e, du, du2, ip, x, 1, info )
143 CALL chkxer(
'ZGTTRS', infot, nout, lerr, ok )
149 CALL zgtrfs(
'/', 0, 0, dl, e, du, dlf, ef, duf, du2, ip, b, 1,
150 $ x, 1, r1, r2, w, rw, info )
151 CALL chkxer(
'ZGTRFS', infot, nout, lerr, ok )
153 CALL zgtrfs(
'N', -1, 0, dl, e, du, dlf, ef, duf, du2, ip, b,
154 $ 1, x, 1, r1, r2, w, rw, info )
155 CALL chkxer(
'ZGTRFS', infot, nout, lerr, ok )
157 CALL zgtrfs(
'N', 0, -1, dl, e, du, dlf, ef, duf, du2, ip, b,
158 $ 1, x, 1, r1, r2, w, rw, info )
159 CALL chkxer(
'ZGTRFS', infot, nout, lerr, ok )
161 CALL zgtrfs(
'N', 2, 1, dl, e, du, dlf, ef, duf, du2, ip, b, 1,
162 $ x, 2, r1, r2, w, rw, info )
163 CALL chkxer(
'ZGTRFS', infot, nout, lerr, ok )
165 CALL zgtrfs(
'N', 2, 1, dl, e, du, dlf, ef, duf, du2, ip, b, 2,
166 $ x, 1, r1, r2, w, rw, info )
167 CALL chkxer(
'ZGTRFS', infot, nout, lerr, ok )
173 CALL zgtcon(
'/', 0, dl, e, du, du2, ip, anorm, rcond, w,
175 CALL chkxer(
'ZGTCON', infot, nout, lerr, ok )
177 CALL zgtcon(
'I', -1, dl, e, du, du2, ip, anorm, rcond, w,
179 CALL chkxer(
'ZGTCON', infot, nout, lerr, ok )
181 CALL zgtcon(
'I', 0, dl, e, du, du2, ip, -anorm, rcond, w,
183 CALL chkxer(
'ZGTCON', infot, nout, lerr, ok )
185 ELSE IF(
lsamen( 2, c2,
'PT' ) )
THEN
194 CALL zpttrf( -1, d, e, info )
195 CALL chkxer(
'ZPTTRF', infot, nout, lerr, ok )
201 CALL zpttrs(
'/', 1, 0, d, e, x, 1, info )
202 CALL chkxer(
'ZPTTRS', infot, nout, lerr, ok )
204 CALL zpttrs(
'U', -1, 0, d, e, x, 1, info )
205 CALL chkxer(
'ZPTTRS', infot, nout, lerr, ok )
207 CALL zpttrs(
'U', 0, -1, d, e, x, 1, info )
208 CALL chkxer(
'ZPTTRS', infot, nout, lerr, ok )
210 CALL zpttrs(
'U', 2, 1, d, e, x, 1, info )
211 CALL chkxer(
'ZPTTRS', infot, nout, lerr, ok )
217 CALL zptrfs(
'/', 1, 0, d, e, df, ef, b, 1, x, 1, r1, r2, w,
219 CALL chkxer(
'ZPTRFS', infot, nout, lerr, ok )
221 CALL zptrfs(
'U', -1, 0, d, e, df, ef, b, 1, x, 1, r1, r2, w,
223 CALL chkxer(
'ZPTRFS', infot, nout, lerr, ok )
225 CALL zptrfs(
'U', 0, -1, d, e, df, ef, b, 1, x, 1, r1, r2, w,
227 CALL chkxer(
'ZPTRFS', infot, nout, lerr, ok )
229 CALL zptrfs(
'U', 2, 1, d, e, df, ef, b, 1, x, 2, r1, r2, w,
231 CALL chkxer(
'ZPTRFS', infot, nout, lerr, ok )
233 CALL zptrfs(
'U', 2, 1, d, e, df, ef, b, 2, x, 1, r1, r2, w,
235 CALL chkxer(
'ZPTRFS', infot, nout, lerr, ok )
241 CALL zptcon( -1, d, e, anorm, rcond, rw, info )
242 CALL chkxer(
'ZPTCON', infot, nout, lerr, ok )
244 CALL zptcon( 0, d, e, -anorm, rcond, rw, info )
245 CALL chkxer(
'ZPTCON', infot, nout, lerr, ok )
250 CALL alaesm( path, ok, nout )