Struct MergerIC 

pub struct MergerIC {
    pub body1: Box<dyn IsolatedEquilibrium>,
    pub mass1: Decimal,
    pub body2: Box<dyn IsolatedEquilibrium>,
    pub mass2: Decimal,
    pub separation: [f64; 3],
    pub relative_velocity: [f64; 3],
    pub impact_parameter: Decimal,
    /* private fields */
}

Binary merger initial conditions: two displaced and boosted equilibria.

The combined distribution is f(x,v) = f1(x-x1, v-v1) + f2(x-x2, v-v2), exact for collisionless systems before the halos begin to interact.

Fields

body1: Box<dyn IsolatedEquilibrium>

First body’s equilibrium distribution (e.g. Plummer, King, NFW).

mass1: Decimal

Total mass of the first body.

body2: Box<dyn IsolatedEquilibrium>

Second body’s equilibrium distribution.

mass2: Decimal

Total mass of the second body.

separation: [f64; 3]

Initial separation vector between the two centres (body2 - body1).

relative_velocity: [f64; 3]

Relative velocity of body 2 with respect to body 1.

impact_parameter: Decimal

Transverse offset for off-axis encounters; zero gives a head-on merger.

Implementations

impl MergerIC

pub fn new( body1: Box<dyn IsolatedEquilibrium>, mass1: f64, body2: Box<dyn IsolatedEquilibrium>, mass2: f64, separation: [f64; 3], relative_velocity: [f64; 3], impact_parameter: f64, ) -> Self

Create a MergerIC from f64 parameters (backward-compatible).

pub fn new_decimal( body1: Box<dyn IsolatedEquilibrium>, mass1: Decimal, body2: Box<dyn IsolatedEquilibrium>, mass2: Decimal, separation: [f64; 3], relative_velocity: [f64; 3], impact_parameter: Decimal, ) -> Self

Create a MergerIC from Decimal parameters (exact config).

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

Sample both components on the grid and sum. Body 1 is centred at (-sep/2, 0, 0) with velocity (-v_rel/2, 0, 0). Body 2 is centred at (+sep/2, 0, 0) with velocity (+v_rel/2, 0, 0).

pub fn validate(&self, domain: &Domain) -> Result<()>

Check that both components fit within the domain.