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 )
subroutine zpttrf(N, D, E, INFO)
ZPTTRF
subroutine zpttrs(UPLO, N, NRHS, D, E, B, LDB, INFO)
ZPTTRS
subroutine zgttrf(N, DL, D, DU, DU2, IPIV, INFO)
ZGTTRF
subroutine zptcon(N, D, E, ANORM, RCOND, RWORK, INFO)
ZPTCON
subroutine zgttrs(TRANS, N, NRHS, DL, D, DU, DU2, IPIV, B, LDB, INFO)
ZGTTRS
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 zgtrfs(TRANS, N, NRHS, DL, D, DU, DLF, DF, DUF, DU2, IPIV, B, LDB, X, LDX, FERR, BERR, WORK, RWORK, INFO)
ZGTRFS
subroutine zptrfs(UPLO, N, NRHS, D, E, DF, EF, B, LDB, X, LDX, FERR, BERR, WORK, RWORK, INFO)
ZPTRFS
logical function lsamen(N, CA, CB)
LSAMEN
subroutine zerrgt(PATH, NUNIT)
ZERRGT
subroutine alaesm(PATH, OK, NOUT)
ALAESM
subroutine zgtcon(NORM, N, DL, D, DU, DU2, IPIV, ANORM, RCOND, WORK, INFO)
ZGTCON