94 SUBROUTINE zpptri( UPLO, N, AP, INFO )
113 parameter ( one = 1.0d+0 )
117 INTEGER J, JC, JJ, JJN
123 EXTERNAL lsame, zdotc
136 upper = lsame( uplo,
'U' )
137 IF( .NOT.upper .AND. .NOT.lsame( uplo,
'L' ) )
THEN
139 ELSE IF( n.LT.0 )
THEN
143 CALL xerbla(
'ZPPTRI', -info )
154 CALL ztptri( uplo,
'Non-unit', n, ap, info )
166 $
CALL zhpr(
'Upper', j-1, one, ap( jc ), 1, ap )
168 CALL zdscal( j, ajj, ap( jc ), 1 )
178 ap( jj ) = dble( zdotc( n-j+1, ap( jj ), 1, ap( jj ), 1 ) )
180 $
CALL ztpmv(
'Lower',
'Conjugate transpose',
'Non-unit',
181 $ n-j, ap( jjn ), ap( jj+1 ), 1 )
subroutine ztptri(UPLO, DIAG, N, AP, INFO)
ZTPTRI
subroutine xerbla(SRNAME, INFO)
XERBLA
subroutine ztpmv(UPLO, TRANS, DIAG, N, AP, X, INCX)
ZTPMV
subroutine zhpr(UPLO, N, ALPHA, X, INCX, AP)
ZHPR
subroutine zpptri(UPLO, N, AP, INFO)
ZPPTRI
subroutine zdscal(N, DA, ZX, INCX)
ZDSCAL