72 parameter ( nmax = 2 )
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 cpttrf(N, D, E, INFO)
CPTTRF
subroutine cptcon(N, D, E, ANORM, RCOND, RWORK, INFO)
CPTCON
subroutine alaesm(PATH, OK, NOUT)
ALAESM
subroutine cgttrf(N, DL, D, DU, DU2, IPIV, INFO)
CGTTRF
subroutine cgtcon(NORM, N, DL, D, DU, DU2, IPIV, ANORM, RCOND, WORK, INFO)
CGTCON
subroutine cpttrs(UPLO, N, NRHS, D, E, B, LDB, INFO)
CPTTRS
subroutine cgttrs(TRANS, N, NRHS, DL, D, DU, DU2, IPIV, B, LDB, INFO)
CGTTRS
subroutine chkxer(SRNAMT, INFOT, NOUT, LERR, OK)
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 cgtrfs(TRANS, N, NRHS, DL, D, DU, DLF, DF, DUF, DU2, IPIV, B, LDB, X, LDX, FERR, BERR, WORK, RWORK, INFO)
CGTRFS