74 DOUBLE PRECISION ANORM, RCOND
77 INTEGER IP( NMAX ), IW( NMAX )
78 DOUBLE PRECISION B( NMAX ), C( NMAX ), CF( NMAX ), D( NMAX ),
79 $ DF( NMAX ), E( NMAX ), EF( NMAX ), F( NMAX ),
80 $ R1( NMAX ), R2( NMAX ), W( NMAX ), X( NMAX )
96 COMMON / infoc / infot, nout, ok, lerr
97 COMMON / srnamc / srnamt
102 WRITE( nout, fmt = * )
115 IF( lsamen( 2, c2,
'GT' ) )
THEN
123 CALL dgttrf( -1, c, d, e, f, ip, info )
124 CALL chkxer(
'DGTTRF', infot, nout, lerr, ok )
130 CALL dgttrs(
'/', 0, 0, c, d, e, f, ip, x, 1, info )
131 CALL chkxer(
'DGTTRS', infot, nout, lerr, ok )
133 CALL dgttrs(
'N', -1, 0, c, d, e, f, ip, x, 1, info )
134 CALL chkxer(
'DGTTRS', infot, nout, lerr, ok )
136 CALL dgttrs(
'N', 0, -1, c, d, e, f, ip, x, 1, info )
137 CALL chkxer(
'DGTTRS', infot, nout, lerr, ok )
139 CALL dgttrs(
'N', 2, 1, c, d, e, f, ip, x, 1, info )
140 CALL chkxer(
'DGTTRS', infot, nout, lerr, ok )
146 CALL dgtrfs(
'/', 0, 0, c, d, e, cf, df, ef, f, ip, b, 1, x, 1,
147 $ r1, r2, w, iw, info )
148 CALL chkxer(
'DGTRFS', infot, nout, lerr, ok )
150 CALL dgtrfs(
'N', -1, 0, c, d, e, cf, df, ef, f, ip, b, 1, x,
151 $ 1, r1, r2, w, iw, info )
152 CALL chkxer(
'DGTRFS', infot, nout, lerr, ok )
154 CALL dgtrfs(
'N', 0, -1, c, d, e, cf, df, ef, f, ip, b, 1, x,
155 $ 1, r1, r2, w, iw, info )
156 CALL chkxer(
'DGTRFS', infot, nout, lerr, ok )
158 CALL dgtrfs(
'N', 2, 1, c, d, e, cf, df, ef, f, ip, b, 1, x, 2,
159 $ r1, r2, w, iw, info )
160 CALL chkxer(
'DGTRFS', infot, nout, lerr, ok )
162 CALL dgtrfs(
'N', 2, 1, c, d, e, cf, df, ef, f, ip, b, 2, x, 1,
163 $ r1, r2, w, iw, info )
164 CALL chkxer(
'DGTRFS', infot, nout, lerr, ok )
170 CALL dgtcon(
'/', 0, c, d, e, f, ip, anorm, rcond, w, iw,
172 CALL chkxer(
'DGTCON', infot, nout, lerr, ok )
174 CALL dgtcon(
'I', -1, c, d, e, f, ip, anorm, rcond, w, iw,
176 CALL chkxer(
'DGTCON', infot, nout, lerr, ok )
178 CALL dgtcon(
'I', 0, c, d, e, f, ip, -anorm, rcond, w, iw,
180 CALL chkxer(
'DGTCON', infot, nout, lerr, ok )
182 ELSE IF( lsamen( 2, c2,
'PT' ) )
THEN
191 CALL dpttrf( -1, d, e, info )
192 CALL chkxer(
'DPTTRF', infot, nout, lerr, ok )
198 CALL dpttrs( -1, 0, d, e, x, 1, info )
199 CALL chkxer(
'DPTTRS', infot, nout, lerr, ok )
201 CALL dpttrs( 0, -1, d, e, x, 1, info )
202 CALL chkxer(
'DPTTRS', infot, nout, lerr, ok )
204 CALL dpttrs( 2, 1, d, e, x, 1, info )
205 CALL chkxer(
'DPTTRS', infot, nout, lerr, ok )
211 CALL dptrfs( -1, 0, d, e, df, ef, b, 1, x, 1, r1, r2, w, info )
212 CALL chkxer(
'DPTRFS', infot, nout, lerr, ok )
214 CALL dptrfs( 0, -1, d, e, df, ef, b, 1, x, 1, r1, r2, w, info )
215 CALL chkxer(
'DPTRFS', infot, nout, lerr, ok )
217 CALL dptrfs( 2, 1, d, e, df, ef, b, 1, x, 2, r1, r2, w, info )
218 CALL chkxer(
'DPTRFS', infot, nout, lerr, ok )
220 CALL dptrfs( 2, 1, d, e, df, ef, b, 2, x, 1, r1, r2, w, info )
221 CALL chkxer(
'DPTRFS', infot, nout, lerr, ok )
227 CALL dptcon( -1, d, e, anorm, rcond, w, info )
228 CALL chkxer(
'DPTCON', infot, nout, lerr, ok )
230 CALL dptcon( 0, d, e, -anorm, rcond, w, info )
231 CALL chkxer(
'DPTCON', infot, nout, lerr, ok )
236 CALL alaesm( path, ok, nout )