Struct YoshidaSplitting 

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

Fourth-order symplectic time integrator via Yoshida splitting.

Composes 7 sub-steps (4 drifts, 3 kicks with Poisson solves) using the Yoshida (1990) triple-jump coefficients to achieve O(dt^4) accuracy while preserving symplecticity.

Fields

g: f64

Gravitational constant G used when solving the Poisson equation.

Implementations

impl YoshidaSplitting

pub fn new(g: f64) -> Self

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

Trait Implementations

impl TimeIntegrator for YoshidaSplitting

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 Yoshida splitting.

Executes the 7-sub-step sequence (4 drifts + 3 kick-with-Poisson-solve) and applies SpectralV hypercollision damping after each kick. Returns end-of-step products.

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.