Function density_summation 

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

Compute per-particle densities using SPH summation.

rho_i = sum_j m_j * W(|r_i - r_j|, h)