pub fn kernel_poly6(r: f64, h: f64) -> f64
Poly6 SPH kernel: W_poly6(r, h).
Returns 315/(64·π·h⁹) · (h²−r²)³ for r < h, else 0.
315/(64·π·h⁹) · (h²−r²)³
r < h
0