reimplemented 1/(r+d_2) smoothing, NlogN still without comments