pub fn sph_kernel_poly6(r_sq: f32, h: f32) -> f32
Poly6 kernel W(r,h) = (315/64πh^9)*(h^2-r^2)^3 for |r|<=h