anchor.diff:  
     NPFIU = hecMAT%NPU+NPUf1
     NPFIL = hecMAT%NPL+NPLf1
  -  
  +
       FI1L= 0
       FI1U= 0
   
  
positive.code:     enddo
      deallocate(vv)
      deallocate(zz)
      
      do i =1,N
        SAINVD(9*i-8) = 1.0d0/SAINVD(9*i-8)
        SAINVD(9*i-4) = 1.0d0/SAINVD(9*i-4)
  
negative.code:         vv(3*in  )= vv(3*in  ) + AU(9*j-6)*X1 + AU(9*j-3)*X2 + AU(9*j  )*X3
        enddo
      enddo
  
      !{d}={v^t}{z_j}
  
      !dtemp(1) = SAINVD(9*itr-8)
      !dtemp(2) = SAINVD(9*itr-4)
  
  !$OMP PARALLEL DEFAULT(NONE) &
  !$OMP&PRIVATE(i,j,jS,jE,in,X1,X2,X3) &
  !$OMP&FIRSTPRIVATE(vv) &
  !$OMP&SHARED(N,itr,SAINVD,SAINVL,inumFI1L,FI1L)
  !$OMP DO
      do i=itr,N
        SAINVD(9*i-8) = vv(3*i-2)
        SAINVD(9*i-4) = vv(3*i-2)*SAINVD(9*i-7)   + vv(3*i-1)
  