anchor.diff:      integer(kind=kint ), pointer :: INL(:), INU(:)
       integer(kind=kint ), pointer :: IAL(:)
       integer(kind=kint ), pointer :: IAU(:)
  -
  -    real   (kind=kreal):: ALUtmp(3,3)
  -    integer(kind=kint ):: ip
   
       N = hecMAT%N
       NP = hecMAT%NP
  
positive.code:           endif
          enddo
        enddo
      enddo
  
      deallocate (IW1, IW2)
  
negative.code: 
          enddo
        enddo
  
        call ILU1a33 (DkINV, &
             Dlu0(9*i-8), Dlu0(9*i-7), Dlu0(9*i-6), &
             Dlu0(9*i-5), Dlu0(9*i-4), Dlu0(9*i-3), &
             Dlu0(9*i-2), Dlu0(9*i-1), Dlu0(9*i  ))
        Dlu0(9*i-8)= DkINV(1,1)
        Dlu0(9*i-7)= DkINV(1,2)
        Dlu0(9*i-6)= DkINV(1,3)
        Dlu0(9*i-5)= DkINV(2,1)
        Dlu0(9*i-4)= DkINV(2,2)
        Dlu0(9*i-3)= DkINV(2,3)
        Dlu0(9*i-2)= DkINV(3,1)
        Dlu0(9*i-1)= DkINV(3,2)
        Dlu0(9*i  )= DkINV(3,3)
      enddo
  
      deallocate (IW1, IW2)
  