Struct TidalIC 

pub struct TidalIC {
    pub host_potential: Box<dyn Fn([f64; 3]) -> f64 + Send + Sync>,
    pub progenitor: Box<dyn IsolatedEquilibrium>,
    pub progenitor_position: [f64; 3],
    pub progenitor_velocity: [f64; 3],
}

Tidal stream initial conditions: a satellite on an orbit in a fixed host potential.

The progenitor’s distribution function is sampled in its own rest frame and placed at the specified position and velocity within the host.

Fields

host_potential: Box<dyn Fn([f64; 3]) -> f64 + Send + Sync>

Fixed external host potential Phi_host(x). Does not evolve self-consistently.

progenitor: Box<dyn IsolatedEquilibrium>

Equilibrium model for the satellite (e.g. Plummer, King).

progenitor_position: [f64; 3]

Initial centre-of-mass position of the satellite in the host frame.

progenitor_velocity: [f64; 3]

Initial centre-of-mass velocity of the satellite in the host frame.

Implementations

impl TidalIC

pub fn new( host_potential: Box<dyn Fn([f64; 3]) -> f64 + Send + Sync>, progenitor: Box<dyn IsolatedEquilibrium>, progenitor_position: [f64; 3], progenitor_velocity: [f64; 3], ) -> Self

Create a tidal-stream IC from a host potential, satellite model, and orbital state.

pub fn sample_on_grid( &self, domain: &Domain, progress: Option<&StepProgress>, ) -> PhaseSpaceSnapshot

Sample progenitor f centred on (progenitor_position, progenitor_velocity).

At each (x, v) grid point: compute position and velocity relative to the progenitor centre, evaluate energy in progenitor’s rest frame, then use the progenitor’s distribution function.

pub fn escape_velocity(&self, r: f64) -> f64

Compute escape velocity from host potential at given galactocentric radius r. Particles with velocity > v_esc at radius r are tidal debris.