anchor.diff:          enddo
         endif
       enddo
  -    
  +
       !------------------------------ iitr = 1 ----------------------------------------
  -    
  +
       zz(:) = 0.0d0
       vv(:) = 0.0d0
  -    
  +
       !{v}=[A]{zi}
   
       zz(3*itr-2)= SAINVD(9*itr-7)
       zz(3*itr-1)= SAINVD(9*itr-4)
       zz(3*itr  )= SAINVD(9*itr-1)
  -    
  +
       zz(3*itr-1)= 1.0d0
  -    
  +
       jS= inumFI1L(itr-1) + 1
       jE= inumFI1L(itr  )
       do j= jS, jE
  
positive.code: 
    NPFIU = hecMAT%NPU+NPUf1
    NPFIL = hecMAT%NPL+NPLf1
    
      FI1L= 0
      FI1U= 0
  
  
negative.code:       enddo flag1
      enddo
    end subroutine hecmw_sainv_make_u_33
  
  !C***
  !C*** FORM_ILU1_33
  !C*** form ILU(1) matrix
    subroutine FORM_ILU0_SAINV_33(hecMAT)
      implicit none
      type(hecmwST_matrix) :: hecMAT
  
      allocate (inumFI1L(0:hecMAT%NP), inumFI1U(0:hecMAT%NP))
      allocate (FI1L (hecMAT%NPL), FI1U (hecMAT%NPU))
  
      inumFI1L = 0
      inumFI1U = 0
      FI1L = 0
      FI1U = 0
  
      inumFI1L = hecMAT%indexL
      inumFI1U = hecMAT%indexU
      FI1L = hecMAT%itemL
      FI1U = hecMAT%itemU
  
      NPFIU = hecMAT%NPU
      NPFIL = hecMAT%NPL
  
    end subroutine FORM_ILU0_SAINV_33
  
  !C***
  !C*** FORM_ILU1_33
  !C*** form ILU(1) matrix
  