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 zpptri(UPLO, N, AP, INFO)
ZPPTRI
subroutine zdscal(N, DA, ZX, INCX)
ZDSCAL
subroutine ztpmv(UPLO, TRANS, DIAG, N, AP, X, INCX)
ZTPMV
subroutine xerbla(SRNAME, INFO)
XERBLA
subroutine ztptri(UPLO, DIAG, N, AP, INFO)
ZTPTRI
subroutine zhpr(UPLO, N, ALPHA, X, INCX, AP)
ZHPR