110 SUBROUTINE dpttrs( N, NRHS, D, E, B, LDB, INFO )
118 INTEGER info, ldb, n, nrhs
121 DOUBLE PRECISION b( ldb, * ), d( * ), e( * )
146 ELSE IF( nrhs.LT.0 )
THEN
148 ELSE IF( ldb.LT.max( 1, n ) )
THEN
152 CALL
xerbla(
'DPTTRS', -info )
158 IF( n.EQ.0 .OR. nrhs.EQ.0 )
166 nb = max( 1,
ilaenv( 1,
'DPTTRS',
' ', n, nrhs, -1, -1 ) )
169 IF( nb.GE.nrhs )
THEN
170 CALL
dptts2( n, nrhs, d, e,
b, ldb )
172 DO 10
j = 1, nrhs, nb
173 jb = min( nrhs-
j+1, nb )
174 CALL
dptts2( n, jb, d, e,
b( 1,
j ), ldb )
subroutine xerbla(SRNAME, INFO)
XERBLA
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real b(3) integer i
subroutine dptts2(N, NRHS, D, E, B, LDB)
DPTTS2 solves a tridiagonal system of the form AX=B using the L D LH factorization computed by spttrf...
subroutine dpttrs(N, NRHS, D, E, B, LDB, INFO)
DPTTRS
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j
integer function ilaenv(ISPEC, NAME, OPTS, N1, N2, N3, N4)