Struct StrangSplitting 

pub struct StrangSplitting {
    pub g: f64,
    /* private fields */
}

Second-order symplectic time integrator via Strang splitting.

Executes drift(dt/2) -> Poisson solve -> kick(dt) -> drift(dt/2) each step. After the kick, applies hypercollision damping when the representation is SpectralV.

Fields

g: f64

Gravitational constant G used when solving the Poisson equation.

Implementations

impl StrangSplitting

pub fn new(g: f64) -> Self

Create a new Strang splitting integrator with the given gravitational constant.

Trait Implementations

impl TimeIntegrator for StrangSplitting

fn advance( &mut self, repr: &mut dyn PhaseSpaceRepr, solver: &dyn PoissonSolver, advector: &dyn Advector, dt: f64, ) -> Result<StepProducts, CausticError>

Advance the phase-space representation by one time step dt using Strang splitting.

Sequence: drift(dt/2) -> Poisson solve -> kick(dt) -> [hypercollision] -> drift(dt/2). Returns the end-of-step density, potential, and acceleration for reuse by diagnostics.

fn max_dt(&self, repr: &dyn PhaseSpaceRepr, cfl_factor: f64) -> f64

Estimate the maximum stable time step from the dynamical time t_dyn = 1/sqrt(G*rho_max).

fn last_step_timings(&self) -> Option<&StepTimings>

Return the timing breakdown from the most recent step.

fn set_progress(&mut self, progress: Arc<StepProgress>)

Attach a shared progress reporter for intra-step phase tracking.

fn suggested_dt(&self) -> Option<f64>

Return a suggested Δt from the adaptive controller (if any). Default returns None; adaptive integrators override this.