anchor.diff:            end if
          end do
       end do
  +    !$OMP END PARALLEL DO
     end subroutine filter_neighbors
   
   end module modphi
  
positive.code:        DEALLOCATE(VSa)
         DEALLOCATE(mask)
      end do
    end subroutine normal_vector_stroke
  
    double precision function dpair(i,j)
  
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))
  