anchor.diff:          SAINVL(9*j-3) = SAINVL(9*j-3)-dd2*zz(3*in  )
         enddo
       endif
  -    
  +
       dd3=SAINVD(9*itr  )
       if(dabs(dd3) > FILTER)then
         SAINVD(9*itr-6)= SAINVD(9*itr-6) - dd3*zz(3*itr-2)
  
positive.code:       zz(3*in-1)= SAINVL(9*j-1)
        zz(3*in  )= SAINVL(9*J  )
      enddo
      
      do i= 1, itr
        X1= zz(3*i-2)
        X2= zz(3*i-1)
  
negative.code:   use m_hecmw_comm_f
    use hecmw_matrix_contact
    use hecmw_matrix_misc
    !$ use omp_lib
  
    private
  
    public:: hecmw_precond_33_SAINV_setup
    public:: hecmw_precond_33_SAINV_apply
    public:: hecmw_precond_33_SAINV_clear
  
    integer(4),parameter :: krealp = 8
  
    integer(kind=kint) :: NPFIU, NPFIL
    integer(kind=kint) :: N
    integer(kind=kint), pointer :: inumFI1L(:) => null()
    integer(kind=kint), pointer :: inumFI1U(:) => null()
    integer(kind=kint), pointer :: FI1L(:) => null()
    integer(kind=kint), pointer :: FI1U(:) => null()
  
    integer(kind=kint), pointer :: indexL(:) => null()
    integer(kind=kint), pointer :: indexU(:) => null()
    integer(kind=kint), pointer :: itemL(:)  => null()
  