Struct salva2d::solver::IISPHSolver

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

A IISPH (Implicit Incompressible 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.

Methods

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

pub fn new() -> Self

Initialize a new IISPH pressure solver.

Trait Implementations

impl<N, KernelDensity, KernelGradient> PressureSolver<N> for IISPHSolver<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.