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 )
subroutine zpptri(UPLO, N, AP, INFO)
ZPPTRI
subroutine zdscal(N, DA, ZX, INCX)
ZDSCAL
subroutine ztpmv(UPLO, TRANS, DIAG, N, AP, X, INCX)
ZTPMV
complex *16 function zdotc(N, ZX, INCX, ZY, INCY)
ZDOTC
subroutine xerbla(SRNAME, INFO)
XERBLA
subroutine ztptri(UPLO, DIAG, N, AP, INFO)
ZTPTRI
logical function lsame(CA, CB)
LSAME
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j
subroutine zhpr(UPLO, N, ALPHA, X, INCX, AP)
ZHPR