94 SUBROUTINE zpptri( UPLO, N, AP, INFO )
113 parameter( one = 1.0d+0 )
117 INTEGER j, jc, jj, jjn
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 )
LOGICAL function lsame(CA, CB)
LSAME
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
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j
COMPLEX *16 function zdotc(N, ZX, INCX, ZY, INCY)
ZDOTC
subroutine zhpr(UPLO, N, ALPHA, X, INCX, AP)
ZHPR