Function density_summation_kernel 

pub fn density_summation_kernel(
    positions: &[[f64; 3]],
    masses: &[f64],
    h: f64,
    kernel: SphKernel,
) -> Vec<f64>

Compute per-particle densities with a specified kernel type.