118 SUBROUTINE stptri( UPLO, DIAG, N, AP, INFO )
137 parameter ( one = 1.0e+0, zero = 0.0e+0 )
140 LOGICAL NOUNIT, UPPER
141 INTEGER J, JC, JCLAST, JJ
156 upper = lsame( uplo,
'U' )
157 nounit = lsame( diag,
'N' )
158 IF( .NOT.upper .AND. .NOT.lsame( uplo,
'L' ) )
THEN
160 ELSE IF( .NOT.nounit .AND. .NOT.lsame( diag,
'U' ) )
THEN
162 ELSE IF( n.LT.0 )
THEN
166 CALL xerbla(
'STPTRI', -info )
177 IF( ap( jj ).EQ.zero )
183 IF( ap( jj ).EQ.zero )
185 jj = jj + n - info + 1
198 ap( jc+j-1 ) = one / ap( jc+j-1 )
206 CALL stpmv(
'Upper',
'No transpose', diag, j-1, ap,
208 CALL sscal( j-1, ajj, ap( jc ), 1 )
219 ap( jc ) = one / ap( jc )
228 CALL stpmv(
'Lower',
'No transpose', diag, n-j,
229 $ ap( jclast ), ap( jc+1 ), 1 )
230 CALL sscal( n-j, ajj, ap( jc+1 ), 1 )
subroutine stptri(UPLO, DIAG, N, AP, INFO)
STPTRI
subroutine stpmv(UPLO, TRANS, DIAG, N, AP, X, INCX)
STPMV
subroutine xerbla(SRNAME, INFO)
XERBLA
subroutine sscal(N, SA, SX, INCX)
SSCAL