81 REAL d( nmax ), df( nmax ), r1( nmax ), r2( nmax ),
83 COMPLEX 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
cgttrf( -1, dl, e, du, du2, ip, info )
127 CALL
chkxer(
'CGTTRF', infot, nout, lerr, ok )
133 CALL
cgttrs(
'/', 0, 0, dl, e, du, du2, ip,
x, 1, info )
134 CALL
chkxer(
'CGTTRS', infot, nout, lerr, ok )
136 CALL
cgttrs(
'N', -1, 0, dl, e, du, du2, ip,
x, 1, info )
137 CALL
chkxer(
'CGTTRS', infot, nout, lerr, ok )
139 CALL
cgttrs(
'N', 0, -1, dl, e, du, du2, ip,
x, 1, info )
140 CALL
chkxer(
'CGTTRS', infot, nout, lerr, ok )
142 CALL
cgttrs(
'N', 2, 1, dl, e, du, du2, ip,
x, 1, info )
143 CALL
chkxer(
'CGTTRS', infot, nout, lerr, ok )
149 CALL
cgtrfs(
'/', 0, 0, dl, e, du, dlf, ef, duf, du2, ip,
b, 1,
150 $
x, 1, r1, r2, w, rw, info )
151 CALL
chkxer(
'CGTRFS', infot, nout, lerr, ok )
153 CALL
cgtrfs(
'N', -1, 0, dl, e, du, dlf, ef, duf, du2, ip,
b,
154 $ 1,
x, 1, r1, r2, w, rw, info )
155 CALL
chkxer(
'CGTRFS', infot, nout, lerr, ok )
157 CALL
cgtrfs(
'N', 0, -1, dl, e, du, dlf, ef, duf, du2, ip,
b,
158 $ 1,
x, 1, r1, r2, w, rw, info )
159 CALL
chkxer(
'CGTRFS', infot, nout, lerr, ok )
161 CALL
cgtrfs(
'N', 2, 1, dl, e, du, dlf, ef, duf, du2, ip,
b, 1,
162 $
x, 2, r1, r2, w, rw, info )
163 CALL
chkxer(
'CGTRFS', infot, nout, lerr, ok )
165 CALL
cgtrfs(
'N', 2, 1, dl, e, du, dlf, ef, duf, du2, ip,
b, 2,
166 $
x, 1, r1, r2, w, rw, info )
167 CALL
chkxer(
'CGTRFS', infot, nout, lerr, ok )
173 CALL
cgtcon(
'/', 0, dl, e, du, du2, ip, anorm, rcond, w,
175 CALL
chkxer(
'CGTCON', infot, nout, lerr, ok )
177 CALL
cgtcon(
'I', -1, dl, e, du, du2, ip, anorm, rcond, w,
179 CALL
chkxer(
'CGTCON', infot, nout, lerr, ok )
181 CALL
cgtcon(
'I', 0, dl, e, du, du2, ip, -anorm, rcond, w,
183 CALL
chkxer(
'CGTCON', infot, nout, lerr, ok )
185 ELSE IF(
lsamen( 2, c2,
'PT' ) )
THEN
194 CALL
cpttrf( -1, d, e, info )
195 CALL
chkxer(
'CPTTRF', infot, nout, lerr, ok )
201 CALL
cpttrs(
'/', 1, 0, d, e,
x, 1, info )
202 CALL
chkxer(
'CPTTRS', infot, nout, lerr, ok )
204 CALL
cpttrs(
'U', -1, 0, d, e,
x, 1, info )
205 CALL
chkxer(
'CPTTRS', infot, nout, lerr, ok )
207 CALL
cpttrs(
'U', 0, -1, d, e,
x, 1, info )
208 CALL
chkxer(
'CPTTRS', infot, nout, lerr, ok )
210 CALL
cpttrs(
'U', 2, 1, d, e,
x, 1, info )
211 CALL
chkxer(
'CPTTRS', infot, nout, lerr, ok )
217 CALL
cptrfs(
'/', 1, 0, d, e, df, ef,
b, 1,
x, 1, r1, r2, w,
219 CALL
chkxer(
'CPTRFS', infot, nout, lerr, ok )
221 CALL
cptrfs(
'U', -1, 0, d, e, df, ef,
b, 1,
x, 1, r1, r2, w,
223 CALL
chkxer(
'CPTRFS', infot, nout, lerr, ok )
225 CALL
cptrfs(
'U', 0, -1, d, e, df, ef,
b, 1,
x, 1, r1, r2, w,
227 CALL
chkxer(
'CPTRFS', infot, nout, lerr, ok )
229 CALL
cptrfs(
'U', 2, 1, d, e, df, ef,
b, 1,
x, 2, r1, r2, w,
231 CALL
chkxer(
'CPTRFS', infot, nout, lerr, ok )
233 CALL
cptrfs(
'U', 2, 1, d, e, df, ef,
b, 2,
x, 1, r1, r2, w,
235 CALL
chkxer(
'CPTRFS', infot, nout, lerr, ok )
241 CALL
cptcon( -1, d, e, anorm, rcond, rw, info )
242 CALL
chkxer(
'CPTCON', infot, nout, lerr, ok )
244 CALL
cptcon( 0, d, e, -anorm, rcond, rw, info )
245 CALL
chkxer(
'CPTCON', infot, nout, lerr, ok )
250 CALL
alaesm( path, ok, nout )
subroutine cpttrs(UPLO, N, NRHS, D, E, B, LDB, INFO)
CPTTRS
subroutine cgtcon(NORM, N, DL, D, DU, DU2, IPIV, ANORM, RCOND, WORK, INFO)
CGTCON
subroutine cptcon(N, D, E, ANORM, RCOND, RWORK, INFO)
CPTCON
subroutine cgttrs(TRANS, N, NRHS, DL, D, DU, DU2, IPIV, B, LDB, INFO)
CGTTRS
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 cpttrf(N, D, E, INFO)
CPTTRF
subroutine cgttrf(N, DL, D, DU, DU2, IPIV, INFO)
CGTTRF
subroutine cerrgt(PATH, NUNIT)
CERRGT
subroutine cptrfs(UPLO, N, NRHS, D, E, DF, EF, B, LDB, X, LDX, FERR, BERR, WORK, RWORK, INFO)
CPTRFS
LOGICAL function lsamen(N, CA, CB)
LSAMEN
subroutine alaesm(PATH, OK, NOUT)
ALAESM
subroutine cgtrfs(TRANS, N, NRHS, DL, D, DU, DLF, DF, DUF, DU2, IPIV, B, LDB, X, LDX, FERR, BERR, WORK, RWORK, INFO)
CGTRFS