72 INTEGER I, INFO, J, RANK
75 REAL A( NMAX, NMAX ), WORK( 2*NMAX )
87 COMMON / infoc / infot, nout, ok, lerr
88 COMMON / srnamc / srnamt
96 WRITE( nout, fmt = * )
102 a( i, j ) = 1.0 / real( i+j )
120 CALL spstrf(
'/', 0, a, 1, piv, rank, -1.0, work, info )
121 CALL chkxer(
'SPSTRF', infot, nout, lerr, ok )
123 CALL spstrf(
'U', -1, a, 1, piv, rank, -1.0, work, info )
124 CALL chkxer(
'SPSTRF', infot, nout, lerr, ok )
126 CALL spstrf(
'U', 2, a, 1, piv, rank, -1.0, work, info )
127 CALL chkxer(
'SPSTRF', infot, nout, lerr, ok )
133 CALL spstf2(
'/', 0, a, 1, piv, rank, -1.0, work, info )
134 CALL chkxer(
'SPSTF2', infot, nout, lerr, ok )
136 CALL spstf2(
'U', -1, a, 1, piv, rank, -1.0, work, info )
137 CALL chkxer(
'SPSTF2', infot, nout, lerr, ok )
139 CALL spstf2(
'U', 2, a, 1, piv, rank, -1.0, work, info )
140 CALL chkxer(
'SPSTF2', infot, nout, lerr, ok )
145 CALL alaesm( path, ok, nout )
subroutine chkxer(SRNAMT, INFOT, NOUT, LERR, OK)
subroutine alaesm(PATH, OK, NOUT)
ALAESM
subroutine spstrf(UPLO, N, A, LDA, PIV, RANK, TOL, WORK, INFO)
SPSTRF computes the Cholesky factorization with complete pivoting of a real symmetric positive semide...
subroutine spstf2(UPLO, N, A, LDA, PIV, RANK, TOL, WORK, INFO)
SPSTF2 computes the Cholesky factorization with complete pivoting of a real symmetric positive semide...
subroutine serrps(PATH, NUNIT)
SERRPS