When two hydro particles are within a small fraction of their smoothing length...
When two hydro particles are within a small fraction of their smoothing length (default 10^-4), add a small delta vector to their distance vector to avoid division by zero