109 SUBROUTINE dpptrs( UPLO, N, NRHS, AP, B, LDB, INFO )
118 INTEGER info, ldb, n, nrhs
121 DOUBLE PRECISION ap( * ),
b( ldb, * )
145 upper =
lsame( uplo,
'U' )
146 IF( .NOT.upper .AND. .NOT.
lsame( uplo,
'L' ) )
THEN
148 ELSE IF( n.LT.0 )
THEN
150 ELSE IF( nrhs.LT.0 )
THEN
152 ELSE IF( ldb.LT.max( 1, n ) )
THEN
156 CALL
xerbla(
'DPPTRS', -info )
162 IF( n.EQ.0 .OR. nrhs.EQ.0 )
173 CALL
dtpsv(
'Upper',
'Transpose',
'Non-unit', n, ap,
178 CALL
dtpsv(
'Upper',
'No transpose',
'Non-unit', n, ap,
189 CALL
dtpsv(
'Lower',
'No transpose',
'Non-unit', n, ap,
194 CALL
dtpsv(
'Lower',
'Transpose',
'Non-unit', n, ap,
subroutine xerbla(SRNAME, INFO)
XERBLA
subroutine dtpsv(UPLO, TRANS, DIAG, N, AP, X, INCX)
DTPSV
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real b(3) integer i
subroutine dpptrs(UPLO, N, NRHS, AP, B, LDB, INFO)
DPPTRS
logical function lsame(CA, CB)
LSAME