94 SUBROUTINE spptri( UPLO, N, AP, INFO )
113 parameter( one = 1.0e+0 )
117 INTEGER j, jc, jj, jjn
133 upper =
lsame( uplo,
'U' )
134 IF( .NOT.upper .AND. .NOT.
lsame( uplo,
'L' ) )
THEN
136 ELSE IF( n.LT.0 )
THEN
140 CALL
xerbla(
'SPPTRI', -info )
151 CALL
stptri( uplo,
'Non-unit', n, ap, info )
164 $ CALL
sspr(
'Upper',
j-1, one, ap( jc ), 1, ap )
166 CALL
sscal(
j, ajj, ap( jc ), 1 )
176 ap( jj ) =
sdot( n-
j+1, ap( jj ), 1, ap( jj ), 1 )
178 $ CALL
stpmv(
'Lower',
'Transpose',
'Non-unit', n-
j,
179 $ ap( jjn ), ap( jj+1 ), 1 )
subroutine stptri(UPLO, DIAG, N, AP, INFO)
STPTRI
real function sdot(N, SX, INCX, SY, INCY)
SDOT
subroutine xerbla(SRNAME, INFO)
XERBLA
logical function lsame(CA, CB)
LSAME
subroutine stpmv(UPLO, TRANS, DIAG, N, AP, X, INCX)
STPMV
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j
subroutine sspr(UPLO, N, ALPHA, X, INCX, AP)
SSPR
subroutine spptri(UPLO, N, AP, INFO)
SPPTRI
subroutine sscal(N, SA, SX, INCX)
SSCAL