pub fn viscosity_force(
positions: &[[f64; 3]],
velocities: &[[f64; 3]],
densities: &[f64],
masses: &[f64],
h: f64,
mu: f64,
) -> Vec<[f64; 3]>Expand description
Compute viscosity forces using the SPH viscosity formulation.
F_i^visc = mu * sum_j m_j (v_j - v_i) / rho_j * laplacian_W(r_ij, h)