76 REAL R1( NMAX ), R2( NMAX ), RW( NMAX )
77 COMPLEX A( NMAX, NMAX ), B( NMAX ), W( NMAX ),
95 COMMON / infoc / infot, nout, ok, lerr
96 COMMON / srnamc / srnamt
101 WRITE( nout, fmt = * )
111 IF( lsamen( 2, c2,
'TR' ) )
THEN
117 CALL ctrtri(
'/',
'N', 0, a, 1, info )
118 CALL chkxer(
'CTRTRI', infot, nout, lerr, ok )
120 CALL ctrtri(
'U',
'/', 0, a, 1, info )
121 CALL chkxer(
'CTRTRI', infot, nout, lerr, ok )
123 CALL ctrtri(
'U',
'N', -1, a, 1, info )
124 CALL chkxer(
'CTRTRI', infot, nout, lerr, ok )
126 CALL ctrtri(
'U',
'N', 2, a, 1, info )
127 CALL chkxer(
'CTRTRI', infot, nout, lerr, ok )
133 CALL ctrti2(
'/',
'N', 0, a, 1, info )
134 CALL chkxer(
'CTRTI2', infot, nout, lerr, ok )
136 CALL ctrti2(
'U',
'/', 0, a, 1, info )
137 CALL chkxer(
'CTRTI2', infot, nout, lerr, ok )
139 CALL ctrti2(
'U',
'N', -1, a, 1, info )
140 CALL chkxer(
'CTRTI2', infot, nout, lerr, ok )
142 CALL ctrti2(
'U',
'N', 2, a, 1, info )
143 CALL chkxer(
'CTRTI2', infot, nout, lerr, ok )
150 CALL ctrtrs(
'/',
'N',
'N', 0, 0, a, 1, x, 1, info )
151 CALL chkxer(
'CTRTRS', infot, nout, lerr, ok )
153 CALL ctrtrs(
'U',
'/',
'N', 0, 0, a, 1, x, 1, info )
154 CALL chkxer(
'CTRTRS', infot, nout, lerr, ok )
156 CALL ctrtrs(
'U',
'N',
'/', 0, 0, a, 1, x, 1, info )
157 CALL chkxer(
'CTRTRS', infot, nout, lerr, ok )
159 CALL ctrtrs(
'U',
'N',
'N', -1, 0, a, 1, x, 1, info )
160 CALL chkxer(
'CTRTRS', infot, nout, lerr, ok )
162 CALL ctrtrs(
'U',
'N',
'N', 0, -1, a, 1, x, 1, info )
163 CALL chkxer(
'CTRTRS', infot, nout, lerr, ok )
170 CALL ctrrfs(
'/',
'N',
'N', 0, 0, a, 1, b, 1, x, 1, r1, r2, w,
172 CALL chkxer(
'CTRRFS', infot, nout, lerr, ok )
174 CALL ctrrfs(
'U',
'/',
'N', 0, 0, a, 1, b, 1, x, 1, r1, r2, w,
176 CALL chkxer(
'CTRRFS', infot, nout, lerr, ok )
178 CALL ctrrfs(
'U',
'N',
'/', 0, 0, a, 1, b, 1, x, 1, r1, r2, w,
180 CALL chkxer(
'CTRRFS', infot, nout, lerr, ok )
182 CALL ctrrfs(
'U',
'N',
'N', -1, 0, a, 1, b, 1, x, 1, r1, r2, w,
184 CALL chkxer(
'CTRRFS', infot, nout, lerr, ok )
186 CALL ctrrfs(
'U',
'N',
'N', 0, -1, a, 1, b, 1, x, 1, r1, r2, w,
188 CALL chkxer(
'CTRRFS', infot, nout, lerr, ok )
190 CALL ctrrfs(
'U',
'N',
'N', 2, 1, a, 1, b, 2, x, 2, r1, r2, w,
192 CALL chkxer(
'CTRRFS', infot, nout, lerr, ok )
194 CALL ctrrfs(
'U',
'N',
'N', 2, 1, a, 2, b, 1, x, 2, r1, r2, w,
196 CALL chkxer(
'CTRRFS', infot, nout, lerr, ok )
198 CALL ctrrfs(
'U',
'N',
'N', 2, 1, a, 2, b, 2, x, 1, r1, r2, w,
200 CALL chkxer(
'CTRRFS', infot, nout, lerr, ok )
206 CALL ctrcon(
'/',
'U',
'N', 0, a, 1, rcond, w, rw, info )
207 CALL chkxer(
'CTRCON', infot, nout, lerr, ok )
209 CALL ctrcon(
'1',
'/',
'N', 0, a, 1, rcond, w, rw, info )
210 CALL chkxer(
'CTRCON', infot, nout, lerr, ok )
212 CALL ctrcon(
'1',
'U',
'/', 0, a, 1, rcond, w, rw, info )
213 CALL chkxer(
'CTRCON', infot, nout, lerr, ok )
215 CALL ctrcon(
'1',
'U',
'N', -1, a, 1, rcond, w, rw, info )
216 CALL chkxer(
'CTRCON', infot, nout, lerr, ok )
218 CALL ctrcon(
'1',
'U',
'N', 2, a, 1, rcond, w, rw, info )
219 CALL chkxer(
'CTRCON', infot, nout, lerr, ok )
225 CALL clatrs(
'/',
'N',
'N',
'N', 0, a, 1, x, scale, rw, info )
226 CALL chkxer(
'CLATRS', infot, nout, lerr, ok )
228 CALL clatrs(
'U',
'/',
'N',
'N', 0, a, 1, x, scale, rw, info )
229 CALL chkxer(
'CLATRS', infot, nout, lerr, ok )
231 CALL clatrs(
'U',
'N',
'/',
'N', 0, a, 1, x, scale, rw, info )
232 CALL chkxer(
'CLATRS', infot, nout, lerr, ok )
234 CALL clatrs(
'U',
'N',
'N',
'/', 0, a, 1, x, scale, rw, info )
235 CALL chkxer(
'CLATRS', infot, nout, lerr, ok )
237 CALL clatrs(
'U',
'N',
'N',
'N', -1, a, 1, x, scale, rw, info )
238 CALL chkxer(
'CLATRS', infot, nout, lerr, ok )
240 CALL clatrs(
'U',
'N',
'N',
'N', 2, a, 1, x, scale, rw, info )
241 CALL chkxer(
'CLATRS', infot, nout, lerr, ok )
245 ELSE IF( lsamen( 2, c2,
'TP' ) )
THEN
251 CALL ctptri(
'/',
'N', 0, a, info )
252 CALL chkxer(
'CTPTRI', infot, nout, lerr, ok )
254 CALL ctptri(
'U',
'/', 0, a, info )
255 CALL chkxer(
'CTPTRI', infot, nout, lerr, ok )
257 CALL ctptri(
'U',
'N', -1, a, info )
258 CALL chkxer(
'CTPTRI', infot, nout, lerr, ok )
264 CALL ctptrs(
'/',
'N',
'N', 0, 0, a, x, 1, info )
265 CALL chkxer(
'CTPTRS', infot, nout, lerr, ok )
267 CALL ctptrs(
'U',
'/',
'N', 0, 0, a, x, 1, info )
268 CALL chkxer(
'CTPTRS', infot, nout, lerr, ok )
270 CALL ctptrs(
'U',
'N',
'/', 0, 0, a, x, 1, info )
271 CALL chkxer(
'CTPTRS', infot, nout, lerr, ok )
273 CALL ctptrs(
'U',
'N',
'N', -1, 0, a, x, 1, info )
274 CALL chkxer(
'CTPTRS', infot, nout, lerr, ok )
276 CALL ctptrs(
'U',
'N',
'N', 0, -1, a, x, 1, info )
277 CALL chkxer(
'CTPTRS', infot, nout, lerr, ok )
279 CALL ctptrs(
'U',
'N',
'N', 2, 1, a, x, 1, info )
280 CALL chkxer(
'CTPTRS', infot, nout, lerr, ok )
286 CALL ctprfs(
'/',
'N',
'N', 0, 0, a, b, 1, x, 1, r1, r2, w, rw,
288 CALL chkxer(
'CTPRFS', infot, nout, lerr, ok )
290 CALL ctprfs(
'U',
'/',
'N', 0, 0, a, b, 1, x, 1, r1, r2, w, rw,
292 CALL chkxer(
'CTPRFS', infot, nout, lerr, ok )
294 CALL ctprfs(
'U',
'N',
'/', 0, 0, a, b, 1, x, 1, r1, r2, w, rw,
296 CALL chkxer(
'CTPRFS', infot, nout, lerr, ok )
298 CALL ctprfs(
'U',
'N',
'N', -1, 0, a, b, 1, x, 1, r1, r2, w,
300 CALL chkxer(
'CTPRFS', infot, nout, lerr, ok )
302 CALL ctprfs(
'U',
'N',
'N', 0, -1, a, b, 1, x, 1, r1, r2, w,
304 CALL chkxer(
'CTPRFS', infot, nout, lerr, ok )
306 CALL ctprfs(
'U',
'N',
'N', 2, 1, a, b, 1, x, 2, r1, r2, w, rw,
308 CALL chkxer(
'CTPRFS', infot, nout, lerr, ok )
310 CALL ctprfs(
'U',
'N',
'N', 2, 1, a, b, 2, x, 1, r1, r2, w, rw,
312 CALL chkxer(
'CTPRFS', infot, nout, lerr, ok )
318 CALL ctpcon(
'/',
'U',
'N', 0, a, rcond, w, rw, info )
319 CALL chkxer(
'CTPCON', infot, nout, lerr, ok )
321 CALL ctpcon(
'1',
'/',
'N', 0, a, rcond, w, rw, info )
322 CALL chkxer(
'CTPCON', infot, nout, lerr, ok )
324 CALL ctpcon(
'1',
'U',
'/', 0, a, rcond, w, rw, info )
325 CALL chkxer(
'CTPCON', infot, nout, lerr, ok )
327 CALL ctpcon(
'1',
'U',
'N', -1, a, rcond, w, rw, info )
328 CALL chkxer(
'CTPCON', infot, nout, lerr, ok )
334 CALL clatps(
'/',
'N',
'N',
'N', 0, a, x, scale, rw, info )
335 CALL chkxer(
'CLATPS', infot, nout, lerr, ok )
337 CALL clatps(
'U',
'/',
'N',
'N', 0, a, x, scale, rw, info )
338 CALL chkxer(
'CLATPS', infot, nout, lerr, ok )
340 CALL clatps(
'U',
'N',
'/',
'N', 0, a, x, scale, rw, info )
341 CALL chkxer(
'CLATPS', infot, nout, lerr, ok )
343 CALL clatps(
'U',
'N',
'N',
'/', 0, a, x, scale, rw, info )
344 CALL chkxer(
'CLATPS', infot, nout, lerr, ok )
346 CALL clatps(
'U',
'N',
'N',
'N', -1, a, x, scale, rw, info )
347 CALL chkxer(
'CLATPS', infot, nout, lerr, ok )
351 ELSE IF( lsamen( 2, c2,
'TB' ) )
THEN
357 CALL ctbtrs(
'/',
'N',
'N', 0, 0, 0, a, 1, x, 1, info )
358 CALL chkxer(
'CTBTRS', infot, nout, lerr, ok )
360 CALL ctbtrs(
'U',
'/',
'N', 0, 0, 0, a, 1, x, 1, info )
361 CALL chkxer(
'CTBTRS', infot, nout, lerr, ok )
363 CALL ctbtrs(
'U',
'N',
'/', 0, 0, 0, a, 1, x, 1, info )
364 CALL chkxer(
'CTBTRS', infot, nout, lerr, ok )
366 CALL ctbtrs(
'U',
'N', 'n
', -1, 0, 0, A, 1, X, 1, INFO )
367 CALL CHKXER( 'ctbtrs', INFOT, NOUT, LERR, OK )
369 CALL CTBTRS( 'u
', 'n
', 'n
', 0, -1, 0, A, 1, X, 1, INFO )
370 CALL CHKXER( 'ctbtrs', INFOT, NOUT, LERR, OK )
372 CALL CTBTRS( 'u
', 'n
', 'n
', 0, 0, -1, A, 1, X, 1, INFO )
373 CALL CHKXER( 'ctbtrs', INFOT, NOUT, LERR, OK )
375 CALL CTBTRS( 'u
', 'n
', 'n
', 2, 1, 1, A, 1, X, 2, INFO )
376 CALL CHKXER( 'ctbtrs', INFOT, NOUT, LERR, OK )
378 CALL CTBTRS( 'u
', 'n
', 'n
', 2, 0, 1, A, 1, X, 1, INFO )
379 CALL CHKXER( 'ctbtrs', INFOT, NOUT, LERR, OK )
385 CALL CTBRFS( '/
', 'n
', 'n
', 0, 0, 0, A, 1, B, 1, X, 1, R1, R2,
387 CALL CHKXER( 'ctbrfs', INFOT, NOUT, LERR, OK )
389 CALL CTBRFS( 'u
', '/
', 'n
', 0, 0, 0, A, 1, B, 1, X, 1, R1, R2,
391 CALL CHKXER( 'ctbrfs', INFOT, NOUT, LERR, OK )
393 CALL CTBRFS( 'u
', 'n
', '/
', 0, 0, 0, A, 1, B, 1, X, 1, R1, R2,
395 CALL CHKXER( 'ctbrfs', INFOT, NOUT, LERR, OK )
397 CALL CTBRFS( 'u
', 'n
', 'n
', -1, 0, 0, A, 1, B, 1, X, 1, R1, R2,
399 CALL CHKXER( 'ctbrfs', INFOT, NOUT, LERR, OK )
401 CALL CTBRFS( 'u
', 'n
', 'n
', 0, -1, 0, A, 1, B, 1, X, 1, R1, R2,
403 CALL CHKXER( 'ctbrfs', INFOT, NOUT, LERR, OK )
405 CALL CTBRFS( 'u
', 'n
', 'n
', 0, 0, -1, A, 1, B, 1, X, 1, R1, R2,
407 CALL CHKXER( 'ctbrfs', INFOT, NOUT, LERR, OK )
409 CALL CTBRFS( 'u
', 'n
', 'n
', 2, 1, 1, A, 1, B, 2, X, 2, R1, R2,
411 CALL CHKXER( 'ctbrfs', INFOT, NOUT, LERR, OK )
413 CALL CTBRFS( 'u
', 'n
', 'n', 2, 1, 1, a, 2, b, 1, x, 2, r1, r2,
415 CALL chkxer(
'CTBRFS', infot, nout, lerr, ok )
417 CALL ctbrfs(
'U',
'N',
'N', 2, 1, 1, a, 2, b, 2, x, 1, r1, r2,
419 CALL chkxer(
'CTBRFS', infot, nout, lerr, ok )
425 CALL ctbcon(
'/',
'U',
'N', 0, 0, a, 1, rcond, w, rw, info )
426 CALL chkxer(
'CTBCON', infot, nout, lerr, ok )
428 CALL ctbcon( '1
', '/
', 'n
', 0, 0, A, 1, RCOND, W, RW, INFO )
429 CALL CHKXER( 'ctbcon', INFOT, NOUT, LERR, OK )
431 CALL CTBCON( '1
', 'u
', '/
', 0, 0, A, 1, RCOND, W, RW, INFO )
432 CALL CHKXER( 'ctbcon', INFOT, NOUT, LERR, OK )
434 CALL CTBCON( '1
', 'u
', 'n
', -1, 0, A, 1, RCOND, W, RW, INFO )
435 CALL CHKXER( 'ctbcon', INFOT, NOUT, LERR, OK )
437 CALL CTBCON( '1
', 'u
', 'n
', 0, -1, A, 1, RCOND, W, RW, INFO )
438 CALL CHKXER( 'ctbcon', INFOT, NOUT, LERR, OK )
440 CALL CTBCON( '1
', 'u
', 'n
', 2, 1, A, 1, RCOND, W, RW, INFO )
441 CALL CHKXER( 'ctbcon', INFOT, NOUT, LERR, OK )
447 CALL CLATBS( '/
', 'n
', 'n',
'N', 0, 0, a, 1, x, scale, rw,
449 CALL chkxer(
'CLATBS', infot, nout, lerr, ok )
451 CALL clatbs(
'U',
'/',
'N',
'N', 0, 0, a, 1, x, scale, rw,
455 CALL CLATBS( 'u
', 'n
', '/
', 'n
', 0, 0, A, 1, X, SCALE, RW,
457 CALL CHKXER( 'clatbs', INFOT, NOUT, LERR, OK )
459 CALL CLATBS( 'u
', 'n
', 'n
', '/
', 0, 0, A, 1, X, SCALE, RW,
461 CALL CHKXER( 'clatbs', INFOT, NOUT, LERR, OK )
463 CALL CLATBS( 'u
', 'n
', 'n
', 'n
', -1, 0, A, 1, X, SCALE, RW,
465 CALL CHKXER( 'clatbs', INFOT, NOUT, LERR, OK )
467 CALL CLATBS( 'u
', 'n
', 'n
', 'n
', 1, -1, A, 1, X, SCALE, RW,
469 CALL CHKXER( 'clatbs', INFOT, NOUT, LERR, OK )
471 CALL CLATBS( 'u
', 'n
', 'n
', 'n
', 2, 1, A, 1, X, SCALE, RW,
473 CALL CHKXER( 'clatbs', INFOT, NOUT, LERR, OK )
478 CALL ALAESM( PATH, OK, NOUT )