Struct salva2d::solver::DFSPHSolver

pub struct DFSPHSolver<N: RealField, KernelDensity: Kernel = CubicSplineKernel, KernelGradient: Kernel = CubicSplineKernel> {
    pub min_pressure_iter: usize,
    pub max_pressure_iter: usize,
    pub max_density_error: N,
    pub min_divergence_iter: usize,
    pub max_divergence_iter: usize,
    pub max_divergence_error: N,
    // some fields omitted
}

A DFSPH (Divergence Free Smoothed Particle Hydrodynamics) pressure solver.

Fields

min_pressure_iter: usize

Minimum number of iterations that must be executed for pressure resolution.

max_pressure_iter: usize

Maximum number of iterations that must be executed for pressure resolution.

max_density_error: N

Maximum acceptable density error (in percents).

The pressure solver will continue iterating until the density error drops bellow this threshold, or until the maximum number of pressure iterations is reached.

min_divergence_iter: usize

Minimum number of iterations that must be executed for divergence resolution.

max_divergence_iter: usize

Maximum number of iterations that must be executed for divergence resolution.

max_divergence_error: N

Maximum acceptable divergence error (in percents).

The pressure solver will continue iterating until the divergence error drops bellow this threshold, or until the maximum number of pressure iterations is reached.

Methods

impl<N, KernelDensity, KernelGradient> DFSPHSolver<N, KernelDensity, KernelGradient> where
    N: RealField,
    KernelDensity: Kernel,
    KernelGradient: Kernel, 

pub fn new() -> Self

Initialize a new DFSPH pressure solver.

Trait Implementations

impl<N, KernelDensity, KernelGradient> PressureSolver<N> for DFSPHSolver<N, KernelDensity, KernelGradient> where
    N: RealField,
    KernelDensity: Kernel,
    KernelGradient: Kernel, 

fn init_with_fluids(&mut self, fluids: &[Fluid<N>])

Initialize this solver with the given fluids.

fn init_with_boundaries(&mut self, _boundaries: &[Boundary<N>])

Initialize this solver with the given boundaries.

fn predict_advection(
    &mut self,
    timestep: &TimestepManager<N>,
    kernel_radius: N,
    contact_manager: &ContactManager<N>,
    gravity: &Vector<N>,
    fluids: &mut [Fluid<N>],
    boundaries: &[Boundary<N>]
)

Predicts advection with the given gravity.

fn evaluate_kernels(
    &mut self,
    kernel_radius: N,
    contact_manager: &mut ContactManager<N>,
    fluids: &[Fluid<N>],
    boundaries: &[Boundary<N>]
)

Evaluate the SPH kernels for all the contacts in contact_manager.

fn compute_densities(
    &mut self,
    contact_manager: &ContactManager<N>,
    fluids: &[Fluid<N>],
    boundaries: &mut [Boundary<N>]
)

Compute the densities of all the boundary and fluid particles.

fn step(
    &mut self,
    counters: &mut Counters,
    timestep: &mut TimestepManager<N>,
    gravity: &Vector<N>,
    contact_manager: &mut ContactManager<N>,
    kernel_radius: N,
    fluids: &mut [Fluid<N>],
    boundaries: &[Boundary<N>]
)

Solves pressure and non-pressure force for the given fluids and boundaries.