pub fn pressure_force(
positions: &[[f64; 3]],
_velocities: &[[f64; 3]],
densities: &[f64],
pressures: &[f64],
masses: &[f64],
h: f64,
) -> Vec<[f64; 3]>Expand description
Compute pressure forces using the SPH symmetric pressure gradient formulation.
F_i^pressure = -sum_j m_j (p_i/rho_i^2 + p_j/rho_j^2) nabla W(r_ij, h)
Returns a Vec<[f64;3]> of forces, one per particle.