Struct ZeldovichIC 

pub struct ZeldovichIC {
    pub mean_density: f64,
    pub h0: f64,
    pub omega_m: f64,
    pub omega_lambda: f64,
    pub scale_factor_init: f64,
    pub random_seed: u64,
}

Zel’dovich pancake IC: cold dark matter sheet. f(x,v,t₀) = ρ̄·δ³(v − v₀(x)) where v₀ is the Zel’dovich velocity field.

Fields

mean_density: f64

h0: f64

omega_m: f64

omega_lambda: f64

scale_factor_init: f64

random_seed: u64

Implementations

impl ZeldovichIC

pub fn new(mean_density: f64, cosmology: CosmologyParams, seed: u64) -> Self

pub fn new_single_mode( mean_density: f64, amplitude: f64, wavenumber: f64, sigma_v: f64, ) -> ZeldovichSingleMode

Create a simple single-mode Zel’dovich pancake for testing.

v₀(x) = -A·sin(k·x₁) f(x,v) = (ρ̄/S) · exp(-(v - v₀(x))² / (2σ_v²))

where S normalises the Gaussian to unity over velocity space.

pub fn displacement_field(&self, _domain: &Domain) -> [Vec<f64>; 3]

Generate the Zel’dovich displacement field s(q) from P(k) via FFT.

Requires rand crate for Gaussian random field generation. Use ZeldovichSingleMode for single-mode validation tests.

pub fn velocity_field(&self, _domain: &Domain) -> [Vec<f64>; 3]

Zel’dovich velocity v₀(q) = H·f·s (growing mode).

Requires displacement_field which needs rand crate.

pub fn sample_on_grid(&self, _domain: &Domain) -> PhaseSpaceSnapshot

Sample onto 6D grid as a thin Gaussian in velocity centred on v₀(x).

Requires velocity_field which needs rand crate.