72 parameter ( nmax = 2 )
80 INTEGER ip( nmax ), iw( nmax )
81 REAL b( nmax ), c( nmax ), cf( nmax ), d( nmax ),
82 $ df( nmax ), e( nmax ), ef( nmax ), f( nmax ),
83 $ r1( nmax ), r2( nmax ), w( nmax ), x( nmax )
99 COMMON / infoc / infot, nout, ok, lerr
100 COMMON / srnamc / srnamt
105 WRITE( nout, fmt = * )
118 IF(
lsamen( 2, c2,
'GT' ) )
THEN
126 CALL sgttrf( -1, c, d, e, f, ip, info )
127 CALL chkxer(
'SGTTRF', infot, nout, lerr, ok )
133 CALL sgttrs(
'/', 0, 0, c, d, e, f, ip, x, 1, info )
134 CALL chkxer(
'SGTTRS', infot, nout, lerr, ok )
136 CALL sgttrs(
'N', -1, 0, c, d, e, f, ip, x, 1, info )
137 CALL chkxer(
'SGTTRS', infot, nout, lerr, ok )
139 CALL sgttrs(
'N', 0, -1, c, d, e, f, ip, x, 1, info )
140 CALL chkxer(
'SGTTRS', infot, nout, lerr, ok )
142 CALL sgttrs(
'N', 2, 1, c, d, e, f, ip, x, 1, info )
143 CALL chkxer(
'SGTTRS', infot, nout, lerr, ok )
149 CALL sgtrfs(
'/', 0, 0, c, d, e, cf, df, ef, f, ip, b, 1, x, 1,
150 $ r1, r2, w, iw, info )
151 CALL chkxer(
'SGTRFS', infot, nout, lerr, ok )
153 CALL sgtrfs(
'N', -1, 0, c, d, e, cf, df, ef, f, ip, b, 1, x,
154 $ 1, r1, r2, w, iw, info )
155 CALL chkxer(
'SGTRFS', infot, nout, lerr, ok )
157 CALL sgtrfs(
'N', 0, -1, c, d, e, cf, df, ef, f, ip, b, 1, x,
158 $ 1, r1, r2, w, iw, info )
159 CALL chkxer(
'SGTRFS', infot, nout, lerr, ok )
161 CALL sgtrfs(
'N', 2, 1, c, d, e, cf, df, ef, f, ip, b, 1, x, 2,
162 $ r1, r2, w, iw, info )
163 CALL chkxer(
'SGTRFS', infot, nout, lerr, ok )
165 CALL sgtrfs(
'N', 2, 1, c, d, e, cf, df, ef, f, ip, b, 2, x, 1,
166 $ r1, r2, w, iw, info )
167 CALL chkxer(
'SGTRFS', infot, nout, lerr, ok )
173 CALL sgtcon(
'/', 0, c, d, e, f, ip, anorm, rcond, w, iw,
175 CALL chkxer(
'SGTCON', infot, nout, lerr, ok )
177 CALL sgtcon(
'I', -1, c, d, e, f, ip, anorm, rcond, w, iw,
179 CALL chkxer(
'SGTCON', infot, nout, lerr, ok )
181 CALL sgtcon(
'I', 0, c, d, e, f, ip, -anorm, rcond, w, iw,
183 CALL chkxer(
'SGTCON', infot, nout, lerr, ok )
185 ELSE IF(
lsamen( 2, c2,
'PT' ) )
THEN
194 CALL spttrf( -1, d, e, info )
195 CALL chkxer(
'SPTTRF', infot, nout, lerr, ok )
201 CALL spttrs( -1, 0, d, e, x, 1, info )
202 CALL chkxer(
'SPTTRS', infot, nout, lerr, ok )
204 CALL spttrs( 0, -1, d, e, x, 1, info )
205 CALL chkxer(
'SPTTRS', infot, nout, lerr, ok )
207 CALL spttrs( 2, 1, d, e, x, 1, info )
208 CALL chkxer(
'SPTTRS', infot, nout, lerr, ok )
214 CALL sptrfs( -1, 0, d, e, df, ef, b, 1, x, 1, r1, r2, w, info )
215 CALL chkxer(
'SPTRFS', infot, nout, lerr, ok )
217 CALL sptrfs( 0, -1, d, e, df, ef, b, 1, x, 1, r1, r2, w, info )
218 CALL chkxer(
'SPTRFS', infot, nout, lerr, ok )
220 CALL sptrfs( 2, 1, d, e, df, ef, b, 1, x, 2, r1, r2, w, info )
221 CALL chkxer(
'SPTRFS', infot, nout, lerr, ok )
223 CALL sptrfs( 2, 1, d, e, df, ef, b, 2, x, 1, r1, r2, w, info )
224 CALL chkxer(
'SPTRFS', infot, nout, lerr, ok )
230 CALL sptcon( -1, d, e, anorm, rcond, w, info )
231 CALL chkxer(
'SPTCON', infot, nout, lerr, ok )
233 CALL sptcon( 0, d, e, -anorm, rcond, w, info )
234 CALL chkxer(
'SPTCON', infot, nout, lerr, ok )
239 CALL alaesm( path, ok, nout )
subroutine sgtrfs(TRANS, N, NRHS, DL, D, DU, DLF, DF, DUF, DU2, IPIV, B, LDB, X, LDX, FERR, BERR, WORK, IWORK, INFO)
SGTRFS
subroutine sgtcon(NORM, N, DL, D, DU, DU2, IPIV, ANORM, RCOND, WORK, IWORK, INFO)
SGTCON
subroutine spttrf(N, D, E, INFO)
SPTTRF
subroutine sptcon(N, D, E, ANORM, RCOND, WORK, INFO)
SPTCON
logical function lsamen(N, CA, CB)
LSAMEN
subroutine alaesm(PATH, OK, NOUT)
ALAESM
subroutine sgttrf(N, DL, D, DU, DU2, IPIV, INFO)
SGTTRF
subroutine chkxer(SRNAMT, INFOT, NOUT, LERR, OK)
subroutine sptrfs(N, NRHS, D, E, DF, EF, B, LDB, X, LDX, FERR, BERR, WORK, INFO)
SPTRFS
subroutine sgttrs(TRANS, N, NRHS, DL, D, DU, DU2, IPIV, B, LDB, INFO)
SGTTRS
subroutine spttrs(N, NRHS, D, E, B, LDB, INFO)
SPTTRS