Function sph_compute_forces 

pub fn sph_compute_forces(particles: &mut Vec<SphParticle>, config: &SphConfig)

Compute SPH forces: pressure gradient + viscosity + gravity + surface tension.