78 REAL A( NMAX, NMAX ), AF( NMAX, NMAX ), B( NMAX ),
79 $ R1( NMAX ), R2( NMAX ), W( 3*NMAX ), X( NMAX )
97 COMMON / infoc / infot, nout, ok, lerr
98 COMMON / srnamc / srnamt
106 WRITE( nout, fmt = * )
113 a( i, j ) = 1. / real( i+j )
114 af( i, j ) = 1. / real( i+j )
125 IF( lsamen( 2, c2,
'PO' ) )
THEN
134 CALL spotrf(
'/', 0, a, 1, info )
135 CALL chkxer(
'SPOTRF', infot, nout, lerr, ok )
137 CALL spotrf(
'U', -1, a, 1, info )
138 CALL chkxer(
'SPOTRF', infot, nout, lerr, ok )
140 CALL spotrf(
'U', 2, a, 1, info )
141 CALL chkxer(
'SPOTRF', infot, nout, lerr, ok )
147 CALL spotf2(
'/', 0, a, 1, info )
148 CALL chkxer(
'SPOTF2', infot, nout, lerr, ok )
150 CALL spotf2(
'U', -1, a, 1, info )
151 CALL chkxer(
'SPOTF2', infot, nout, lerr, ok )
153 CALL spotf2(
'U', 2, a, 1, info )
154 CALL chkxer(
'SPOTF2', infot, nout, lerr, ok )
160 CALL spotri(
'/', 0, a, 1, info )
161 CALL chkxer(
'SPOTRI', infot, nout, lerr, ok )
163 CALL spotri(
'U', -1, a, 1, info )
164 CALL chkxer(
'SPOTRI', infot, nout, lerr, ok )
166 CALL spotri(
'U', 2, a, 1, info )
167 CALL chkxer(
'SPOTRI', infot, nout, lerr, ok )
173 CALL spotrs(
'/', 0, 0, a, 1, b, 1, info )
174 CALL chkxer(
'SPOTRS', infot, nout, lerr, ok )
176 CALL spotrs(
'U', -1, 0, a, 1, b, 1, info )
177 CALL chkxer(
'SPOTRS', infot, nout, lerr, ok )
179 CALL spotrs(
'U', 0, -1, a, 1, b, 1, info )
180 CALL chkxer(
'SPOTRS', infot, nout, lerr, ok )
182 CALL spotrs(
'U', 2, 1, a, 1, b, 2, info )
183 CALL chkxer(
'SPOTRS', infot, nout, lerr, ok )
185 CALL spotrs(
'U', 2, 1, a, 2, b, 1, info )
186 CALL chkxer(
'SPOTRS', infot, nout, lerr, ok )
192 CALL sporfs(
'/', 0, 0, a, 1, af, 1, b, 1, x, 1, r1, r2, w, iw,
194 CALL chkxer(
'SPORFS', infot, nout, lerr, ok )
196 CALL sporfs(
'U', -1, 0, a, 1, af, 1, b, 1, x, 1, r1, r2, w,
198 CALL chkxer(
'SPORFS', infot, nout, lerr, ok )
200 CALL sporfs(
'U', 0, -1, a, 1, af, 1, b, 1, x, 1, r1, r2, w,
202 CALL chkxer(
'SPORFS', infot, nout, lerr, ok )
204 CALL sporfs(
'U', 2, 1, a, 1, af, 2, b, 2, x, 2, r1, r2, w, iw,
206 CALL chkxer(
'SPORFS', infot, nout, lerr, ok )
208 CALL sporfs(
'U', 2, 1, a, 2, af, 1, b, 2, x, 2, r1, r2, w, iw,
210 CALL chkxer(
'SPORFS', infot, nout, lerr, ok )
212 CALL sporfs(
'U', 2, 1, a, 2, af, 2, b, 1, x, 2, r1, r2, w, iw,
214 CALL chkxer(
'SPORFS', infot, nout, lerr, ok )
216 CALL sporfs(
'U', 2, 1, a, 2, af, 2, b, 2, x, 1, r1, r2, w, iw,
218 CALL chkxer(
'SPORFS', infot, nout, lerr, ok )
224 CALL spocon(
'/', 0, a, 1, anrm, rcond, w, iw, info )
225 CALL chkxer(
'SPOCON', infot, nout, lerr, ok )
227 CALL spocon(
'U', -1, a, 1, anrm, rcond, w, iw, info )
228 CALL chkxer(
'SPOCON', infot, nout, lerr, ok )
230 CALL spocon(
'U', 2, a, 1, anrm, rcond, w, iw, info )
231 CALL chkxer(
'SPOCON', infot, nout, lerr, ok )
237 CALL spoequ( -1, a, 1, r1, rcond, anrm, info )
238 CALL chkxer(
'SPOEQU', infot, nout, lerr, ok )
240 CALL spoequ( 2, a, 1, r1, rcond, anrm, info )
241 CALL chkxer(
'SPOEQU', infot, nout, lerr, ok )
243 ELSE IF( lsamen( 2, c2,
'PP' ) )
THEN
252 CALL spptrf(
'/', 0, a, info )
253 CALL chkxer(
'SPPTRF', infot, nout, lerr, ok )
255 CALL spptrf(
'U', -1, a, info )
256 CALL chkxer(
'SPPTRF', infot, nout, lerr, ok )
262 CALL spptri(
'/', 0, a, info )
263 CALL chkxer(
'SPPTRI', infot, nout, lerr, ok )
265 CALL spptri(
'U', -1, a, info )
266 CALL chkxer(
'SPPTRI', infot, nout, lerr, ok )
272 CALL spptrs(
'/', 0, 0, a, b, 1, info )
273 CALL chkxer(
'SPPTRS', infot, nout, lerr, ok )
275 CALL spptrs(
'U', -1, 0, a, b, 1, info )
276 CALL chkxer(
'SPPTRS', infot, nout, lerr, ok )
278 CALL spptrs(
'U', 0, -1, a, b, 1, info )
279 CALL chkxer(
'SPPTRS', infot, nout, lerr, ok )
281 CALL spptrs(
'U', 2, 1, a, b, 1, info )
282 CALL chkxer(
'SPPTRS', infot, nout, lerr, ok )
288 CALL spprfs(
'/', 0, 0, a, af, b, 1, x, 1, r1, r2, w, iw,
290 CALL chkxer(
'SPPRFS', infot, nout, lerr, ok )
292 CALL spprfs(
'U', -1, 0, a, af, b, 1, x, 1, r1, r2, w, iw,
294 CALL chkxer(
'SPPRFS', infot, nout, lerr, ok )
296 CALL spprfs(
'U', 0, -1, a, af, b, 1, x, 1, r1, r2, w, iw,
298 CALL chkxer(
'SPPRFS', infot, nout, lerr, ok )
300 CALL spprfs(
'U', 2, 1, a, af, b, 1, x, 2, r1, r2, w, iw,
302 CALL chkxer(
'SPPRFS', infot, nout, lerr, ok )
304 CALL spprfs(
'U', 2, 1, a, af, b, 2, x, 1, r1, r2, w, iw,
306 CALL chkxer(
'SPPRFS', infot, nout, lerr, ok )
312 CALL sppcon(
'/', 0, a, anrm, rcond, w, iw, info )
313 CALL chkxer(
'SPPCON', infot, nout, lerr, ok )
315 CALL sppcon(
'U', -1, a, anrm, rcond, w, iw, info )
316 CALL chkxer(
'SPPCON', infot, nout, lerr, ok )
322 CALL sppequ(
'/', 0, a, r1, rcond, anrm, info )
323 CALL chkxer(
'SPPEQU', infot, nout, lerr, ok )
325 CALL sppequ(
'U', -1, a, r1, rcond, anrm, info )
326 CALL chkxer(
'SPPEQU', infot, nout, lerr, ok )
328 ELSE IF( lsamen( 2, c2,
'PB' ) )
THEN
337 CALL spbtrf(
'/', 0, 0, a, 1, info )
338 CALL chkxer(
'SPBTRF', infot, nout, lerr, ok )
340 CALL spbtrf(
'U', -1, 0, a, 1, info )
341 CALL chkxer(
'SPBTRF', infot, nout, lerr, ok )
343 CALL spbtrf(
'U', 1, -1, a, 1, info )
344 CALL chkxer(
'SPBTRF', infot, nout, lerr, ok )
346 CALL spbtrf(
'U', 2, 1, a, 1, info )
347 CALL chkxer(
'SPBTRF', infot, nout, lerr, ok )
353 CALL spbtf2(
'/', 0, 0, a, 1, info )
354 CALL chkxer(
'SPBTF2', infot, nout, lerr, ok )
356 CALL spbtf2(
'U', -1, 0, a, 1, info )
357 CALL chkxer(
'SPBTF2', infot, nout, lerr, ok )
359 CALL spbtf2(
'U', 1, -1, a, 1, info )
360 CALL chkxer(
'SPBTF2', infot, nout, lerr, ok )
362 CALL spbtf2(
'U', 2, 1, a, 1, info )
363 CALL chkxer(
'SPBTF2', infot, nout, lerr, ok )
369 CALL spbtrs(
'/', 0, 0, 0, a, 1, b, 1, info )
370 CALL chkxer(
'SPBTRS', infot, nout, lerr, ok )
372 CALL spbtrs(
'U', -1, 0, 0, a, 1, b, 1, info )
373 CALL chkxer(
'SPBTRS', infot, nout, lerr, ok )
375 CALL spbtrs(
'U', 1, -1, 0, a, 1, b, 1, info )
376 CALL chkxer(
'SPBTRS', infot, nout, lerr, ok )
378 CALL spbtrs(
'U', 0, 0, -1, a, 1, b, 1, info )
379 CALL chkxer(
'SPBTRS', infot, nout, lerr, ok )
381 CALL spbtrs(
'U', 2, 1, 1, a, 1, b, 1, info )
382 CALL chkxer(
'SPBTRS', infot, nout, lerr, ok )
384 CALL spbtrs(
'U', 2, 0, 1, a, 1, b, 1, info )
385 CALL chkxer(
'SPBTRS', infot, nout, lerr, ok )
391 CALL spbrfs(
'/', 0, 0, 0, a, 1, af, 1, b, 1, x, 1, r1, r2, w,
393 CALL chkxer(
'SPBRFS', infot, nout, lerr, ok )
395 CALL spbrfs(
'U', -1, 0, 0, a, 1, af, 1, b, 1, x, 1, r1, r2, w,
397 CALL chkxer(
'SPBRFS', infot, nout, lerr, ok )
399 CALL spbrfs(
'U', 1, -1, 0, a, 1, af, 1, b, 1, x, 1, r1, r2, w,
401 CALL chkxer(
'SPBRFS', infot, nout, lerr, ok )
403 CALL spbrfs(
'U', 0, 0, -1, a, 1, af, 1, b, 1, x, 1, r1, r2, w,
405 CALL chkxer(
'SPBRFS', infot, nout, lerr, ok )
407 CALL spbrfs(
'U', 2, 1, 1, a, 1, af, 2, b, 2, x, 2, r1, r2, w,
409 CALL chkxer(
'SPBRFS', infot, nout, lerr, ok )
411 CALL spbrfs(
'U', 2, 1, 1, a, 2, af, 1, b, 2, x, 2, r1, r2, w,
413 CALL chkxer(
'SPBRFS', infot, nout, lerr, ok )
415 CALL spbrfs(
'U', 2, 0, 1, a, 1, af, 1, b, 1, x, 2, r1, r2, w,
417 CALL chkxer(
'SPBRFS', infot, nout, lerr, ok )
419 CALL spbrfs(
'U', 2, 0, 1, a, 1, af, 1, b, 2, x, 1, r1, r2, w,
421 CALL chkxer(
'SPBRFS', infot, nout, lerr, ok )
427 CALL spbcon(
'/', 0, 0, a, 1, anrm, rcond, w, iw, info )
428 CALL chkxer(
'SPBCON', infot, nout, lerr, ok )
430 CALL spbcon(
'U', -1, 0, a, 1, anrm, rcond, w, iw, info )
431 CALL chkxer(
'SPBCON', infot, nout, lerr, ok )
433 CALL spbcon(
'U', 1, -1, a, 1, anrm, rcond, w, iw, info )
434 CALL chkxer(
'SPBCON', infot, nout, lerr, ok )
436 CALL spbcon(
'U', 2, 1, a, 1, anrm, rcond, w, iw, info )
437 CALL chkxer(
'SPBCON', infot, nout, lerr, ok )
443 CALL spbequ(
'/', 0, 0, a, 1, r1, rcond, anrm, info )
444 CALL chkxer(
'SPBEQU', infot, nout, lerr, ok )
446 CALL spbequ(
'U', -1, 0, a, 1, r1, rcond, anrm, info )
447 CALL chkxer(
'SPBEQU', infot, nout, lerr, ok )
449 CALL spbequ(
'U', 1, -1, a, 1, r1, rcond, anrm, info )
450 CALL chkxer(
'SPBEQU', infot, nout, lerr, ok )
452 CALL spbequ(
'U', 2, 1, a, 1, r1, rcond, anrm, info )
453 CALL chkxer(
'SPBEQU', infot, nout, lerr, ok )
458 CALL alaesm( path, ok, nout )
subroutine chkxer(SRNAMT, INFOT, NOUT, LERR, OK)
subroutine alaesm(PATH, OK, NOUT)
ALAESM
subroutine spptrf(UPLO, N, AP, INFO)
SPPTRF
subroutine spbrfs(UPLO, N, KD, NRHS, AB, LDAB, AFB, LDAFB, B, LDB, X, LDX, FERR, BERR, WORK, IWORK, INFO)
SPBRFS
subroutine sppequ(UPLO, N, AP, S, SCOND, AMAX, INFO)
SPPEQU
subroutine spptrs(UPLO, N, NRHS, AP, B, LDB, INFO)
SPPTRS
subroutine spprfs(UPLO, N, NRHS, AP, AFP, B, LDB, X, LDX, FERR, BERR, WORK, IWORK, INFO)
SPPRFS
subroutine spbtf2(UPLO, N, KD, AB, LDAB, INFO)
SPBTF2 computes the Cholesky factorization of a symmetric/Hermitian positive definite band matrix (un...
subroutine spbtrs(UPLO, N, KD, NRHS, AB, LDAB, B, LDB, INFO)
SPBTRS
subroutine sppcon(UPLO, N, AP, ANORM, RCOND, WORK, IWORK, INFO)
SPPCON
subroutine spptri(UPLO, N, AP, INFO)
SPPTRI
subroutine spbequ(UPLO, N, KD, AB, LDAB, S, SCOND, AMAX, INFO)
SPBEQU
subroutine spbcon(UPLO, N, KD, AB, LDAB, ANORM, RCOND, WORK, IWORK, INFO)
SPBCON
subroutine spbtrf(UPLO, N, KD, AB, LDAB, INFO)
SPBTRF
subroutine spotri(UPLO, N, A, LDA, INFO)
SPOTRI
subroutine spotf2(UPLO, N, A, LDA, INFO)
SPOTF2 computes the Cholesky factorization of a symmetric/Hermitian positive definite matrix (unblock...
subroutine spotrf(UPLO, N, A, LDA, INFO)
SPOTRF
subroutine spocon(UPLO, N, A, LDA, ANORM, RCOND, WORK, IWORK, INFO)
SPOCON
subroutine sporfs(UPLO, N, NRHS, A, LDA, AF, LDAF, B, LDB, X, LDX, FERR, BERR, WORK, IWORK, INFO)
SPORFS
subroutine spoequ(N, A, LDA, S, SCOND, AMAX, INFO)
SPOEQU
subroutine spotrs(UPLO, N, NRHS, A, LDA, B, LDB, INFO)
SPOTRS
subroutine serrpo(PATH, NUNIT)
SERRPO