pub fn sph_kernel_spiky_grad(r: Vec3, r_len: f32, h: f32) -> Vec3
Spiky kernel gradient (for pressure) — W_spiky = (15/πh^6)*(h-|r|)^3