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 IC: progenitor cluster on an orbit in an external host potential.

Fields

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

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

progenitor: Box<dyn IsolatedEquilibrium>

progenitor_position: [f64; 3]

progenitor_velocity: [f64; 3]

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

pub fn sample_on_grid(&self, domain: &Domain) -> 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.