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 serrps(PATH, NUNIT)
SERRPS
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j
subroutine alaesm(PATH, OK, NOUT)
ALAESM