Function distortion_from_heating

pub fn distortion_from_heating<F>(
    x_grid: &[f64],
    dq_dz: F,
    z_min: f64,
    z_max: f64,
    n_z: usize,
) -> Vec<f64>
where
    F: Fn(f64) -> f64,

Compute the spectral distortion from an arbitrary energy release history.

ΔI(x) = ∫ G_th(x, z’) · d(Q/ρ_γ)/dz’ dz’

Arguments

• x_grid - frequency grid
• dq_dz - function giving d(Δρ/ρ_γ)/dz at each redshift
• z_min - minimum integration redshift
• z_max - maximum integration redshift
• n_z - number of redshift integration points

Returns Δn(x) at each frequency grid point.