anchor.diff:      integer info
       if(self%ludecomposed) then
          if(.not.ispacked(self))then
  -          call ZGETRS( 'N', self%rank, 1, get2dpointer(self), self%rank, self%IPIV, &
  +          call ZGETRS( 'N', self%rank, 1, self%za, self%rank, self%IPIV, &
                  b , self%rank, INFO)
          else
  -          call ZHPTRS( 'U', self%rank, 1, getpackedpointer(self), self%IPIV, b, self%rank, INFO )
  +          call ZHPTRS( 'U', self%rank, 1, self%zap, self%IPIV, b, self%rank, INFO )
          endif
          if (info.ne.0)then
             write(*,*)"error in iterativearpacksecequn Hermiteanmatrixlinsolve ",info
  
positive.code:     allocate(self%ipiv(self%rank))
      if(.not.self%ludecomposed)then
         if(.not.ispacked(self))then
            call ZGETRF( self%rank,self%rank, get2dpointer(self), self%rank, self%IPIV, INFO )
         else
            call ZHPTRF('U',self%rank, getpackedpointer(self), self%IPIV, INFO )
         endif
         if (info.ne.0)then
            write(*,*)"error in iterativearpacksecequn  HermiteanmatrixLU ",info
  
negative.code:     complex(8),intent(inout)::vout(:)
      
      if(self%packed.eqv..true.)then
         call zhpmv("U",self%rank,alpha,self%zap,vin, 1,beta,vout, 1)
      else
         call zhemv("U",self%rank,alpha,self%za,self%rank,vin, 1,beta,vout, 1)
      endif
    end subroutine Hermiteanmatrixvector
    
  