anchor.diff:  
       neighbors(:,:) = 0
   
  +    call OMP_SET_NUM_THREADS(omp_threads)
  +
  +    !$OMP PARALLEL DO &
  +    !$OMP DEFAULT(SHARED) PRIVATE(ki)
       do i=1,numnodesi
          ki = 0
          do j=1,numnodesi
  
positive.code:           end if
         end do
      end do
    end subroutine filter_neighbors
  
  end module modphi
  
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))
  