131 SUBROUTINE dtptrs( UPLO, TRANS, DIAG, N, NRHS, AP, B, LDB, INFO )
139 CHARACTER diag, trans, uplo
140 INTEGER info, ldb, n, nrhs
143 DOUBLE PRECISION ap( * ),
b( ldb, * )
149 DOUBLE PRECISION zero
150 parameter( zero = 0.0d+0 )
153 LOGICAL nounit, upper
171 upper =
lsame( uplo,
'U' )
172 nounit =
lsame( diag,
'N' )
173 IF( .NOT.upper .AND. .NOT.
lsame( uplo,
'L' ) )
THEN
175 ELSE IF( .NOT.
lsame( trans,
'N' ) .AND. .NOT.
176 $
lsame( trans,
'T' ) .AND. .NOT.
lsame( trans,
'C' ) )
THEN
178 ELSE IF( .NOT.nounit .AND. .NOT.
lsame( diag,
'U' ) )
THEN
180 ELSE IF( n.LT.0 )
THEN
182 ELSE IF( nrhs.LT.0 )
THEN
184 ELSE IF( ldb.LT.max( 1, n ) )
THEN
188 CALL
xerbla(
'DTPTRS', -info )
203 IF( ap( jc+info-1 ).EQ.zero )
210 IF( ap( jc ).EQ.zero )
212 jc = jc + n - info + 1
221 CALL
dtpsv( uplo, trans, diag, n, ap,
b( 1,
j ), 1 )
LOGICAL function lsame(CA, CB)
LSAME
subroutine xerbla(SRNAME, INFO)
XERBLA
subroutine dtptrs(UPLO, TRANS, DIAG, N, NRHS, AP, B, LDB, INFO)
DTPTRS
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
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j