anchor.diff:            if (o == i)then
               do idof = 1, NDOF
                 do jdof = 1, NDOF
  -!              SAINVU(NDOF2*(n-1)+NDOF*(idof-1)+jdof)=SAINVL(NDOF2*(n-1)+NDOF*(jdof-1)+idof)
  +                SAINVU(NDOF2*(n-1)+NDOF*(jdof-1)+idof)=SAINVL(NDOF2*(j-1)+NDOF*(idof-1)+jdof)
                 end do 
               end do 
  -            SAINVU(9*n-8)=SAINVL(9*j-8)
  -            SAINVU(9*n-7)=SAINVL(9*j-5)
  -            SAINVU(9*n-6)=SAINVL(9*j-2)
  -            SAINVU(9*n-5)=SAINVL(9*j-7)
  -            SAINVU(9*n-4)=SAINVL(9*j-4)
  -            SAINVU(9*n-3)=SAINVL(9*j-1)
  -            SAINVU(9*n-2)=SAINVL(9*j-6)
  -            SAINVU(9*n-1)=SAINVL(9*j-3)
  -            SAINVU(9*n  )=SAINVL(9*j  )
               n = n + 1
               cycle flag1
             endif
  
positive.code:     FILTER= hecMAT%Rarray(5)
  
      Write(*,"(a,F15.8)")"### SAINV FILTER   :",FILTER
      Write(*,"(a)")"This precoditioner is incomplete implementation. Please use other preconditioner"
  
      call hecmw_sainv_nn(hecMAT)
  
  
negative.code:   use m_hecmw_comm_f
    use hecmw_matrix_contact
    use hecmw_matrix_misc
  
    private
  
  