75 INTEGER i, info, j, rank
78 REAL a( nmax, nmax ), work( 2*nmax )
90 COMMON / infoc / infot, nout, ok, lerr
91 COMMON / srnamc / srnamt
99 WRITE( nout, fmt = * )
105 a( i, j ) = 1.0 /
REAL( i+j )
123 CALL spstrf(
'/', 0, a, 1, piv, rank, -1.0, work, info )
124 CALL chkxer(
'SPSTRF', infot, nout, lerr, ok )
126 CALL spstrf(
'U', -1, a, 1, piv, rank, -1.0, work, info )
127 CALL chkxer(
'SPSTRF', infot, nout, lerr, ok )
129 CALL spstrf(
'U', 2, a, 1, piv, rank, -1.0, work, info )
130 CALL chkxer(
'SPSTRF', infot, nout, lerr, ok )
136 CALL spstf2(
'/', 0, a, 1, piv, rank, -1.0, work, info )
137 CALL chkxer(
'SPSTF2', infot, nout, lerr, ok )
139 CALL spstf2(
'U', -1, a, 1, piv, rank, -1.0, work, info )
140 CALL chkxer(
'SPSTF2', infot, nout, lerr, ok )
142 CALL spstf2(
'U', 2, a, 1, piv, rank, -1.0, work, info )
143 CALL chkxer(
'SPSTF2', infot, nout, lerr, ok )
148 CALL alaesm( path, ok, nout )
subroutine spstf2(UPLO, N, A, LDA, PIV, RANK, TOL, WORK, INFO)
SPSTF2 computes the Cholesky factorization with complete pivoting of a real symmetric or complex Herm...
subroutine spstrf(UPLO, N, A, LDA, PIV, RANK, TOL, WORK, INFO)
SPSTRF
subroutine chkxer(SRNAMT, INFOT, NOUT, LERR, OK)
subroutine alaesm(PATH, OK, NOUT)
ALAESM