Function cubic_spline_kernel_grad 

pub fn cubic_spline_kernel_grad(dx: f64, dy: f64, dz: f64, h: f64) -> [f64; 3]

Gradient of the cubic-spline kernel ∇W(r, h) along the displacement vector.

Returns [dW/dx, dW/dy, dW/dz].