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.

Uses Xorshift64 RNG with Box-Muller transform for Gaussian random field generation in Fourier space, then IFFT to real space.

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

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

f(a) ≈ Ω_m(a)^{0.55} is the logarithmic growth rate.

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

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

f(x,v) = ρ̄ / ((2π)^{3/2} σ_v³) · exp(-|v - v₀(x)|² / (2σ_v²)) where σ_v = 0.3 · max(Δv) keeps the distribution cold but resolved.