Enum splashsurf_lib::ParticleDensityComputationStrategy

pub enum ParticleDensityComputationStrategy {
    Global,
    SynchronizeSubdomains,
    IndependentSubdomains,
}

Available strategies for the computation of the particle densities

Variants

Global

Compute the particle densities globally before performing domain decomposition.

With this approach the particle densities are computed globally on all particles before any domain decomposition is performed.

This approach is guaranteed to lead to consistent results and does not depend on the following decomposition. However, it is also by far the slowest method as global operations (especially the global neighborhood search) are much slower.

SynchronizeSubdomains

Compute particle densities for all particles locally followed by a synchronization step.

This is the recommended approach. The particle densities will be evaluated for all particles per subdomain, possibly in parallel. Afterwards, the values for all non-ghost particles are written to a global array. This happens in a separate step before performing any reconstructions For the following reconstruction procedure, each subdomain will update the densities of its ghost particles from this global array. This ensures that all ghost-particles receive correct density values without requiring to double the width of the ghost-particle margin just to ensure correct values for the actual inner ghost-particles (i.e. in contrast to the completely local approach).

The actual synchronization overhead is relatively low and this approach is often the fastest method.

This approach should always lead consistent results. Only in very rare cases when a particle is not uniquely assigned during domain decomposition this might lead to problems. If you encounter such problems with this approach please report it as a bug.

IndependentSubdomains

Compute densities locally per subdomain without global synchronization.

The particle densities will be evaluated per subdomain on-the-fly just before the reconstruction of the subdomain happens. In order to compute correct densities for the ghost particles of each subdomain it is required that the ghost-particle margin is at least two times the kernel compact support radius. This may add a lot of additional ghost-particles to each subdomain.

If the ghost-particle margin is not set wide enough, this may lead to density differences on subdomain boundaries. Otherwise this approach robust with respect to the classification of particles into the subdomains.

Trait Implementations

impl Clone for ParticleDensityComputationStrategy

fn clone(&self) -> ParticleDensityComputationStrategy

Returns a copy of the value.

fn clone_from(&mut self, source: &Self)

Performs copy-assignment from source.

impl Debug for ParticleDensityComputationStrategy

fn fmt(&self, f: &mut Formatter<'_>) -> Result

Formats the value using the given formatter.

impl Copy for ParticleDensityComputationStrategy