anchor.diff:         DEALLOCATE(VSa)
          DEALLOCATE(mask)
       end do
  +    !$OMP END PARALLEL DO
     end subroutine normal_vector_stroke
   
     double precision function dpair(i,j)
  
positive.code: 
      neighbors(:,:) = 0
  
      do i=1,numnodesi
         ki = 0
         do j=1,numnodesi
  
negative.code: 
      tmpVector(:) = 0.
  
      call OMP_SET_NUM_THREADS(omp_threads)
  
      !$OMP PARALLEL DO &
      !$OMP DEFAULT(SHARED) PRIVATE(tmpVector,ind,mask,VSa)
      do i=1, numnodesi
         ind = minloc(neighbors(i,:),dim=1) - 1
         allocate (mask(ind))
  