gravitational 0.2.6

Gravitational waves from nanohertz to terahertz: strain conventions, detector noise curves (LIGO, Einstein Telescope, LISA), inspirals and binary pulsars, Kerr orbits, ringdown and merger remnants, neutron-star tides, radiated energy, memory, propagation and the graviton mass, graviton counts, stochastic backgrounds and inflation, high-frequency limits, gravity sensors as detectors, laboratory sources, gravitational lensing and ΛCDM distances.
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45

#![doc = include_str!("../README.md")]

pub mod cosmology;
pub mod curves;
pub mod detectors;
pub mod energy;
pub mod hf;
pub mod inspiral;
pub mod laboratory;
pub mod lensing;
pub mod memory;
pub mod propagation;
pub mod quanta;
pub mod remnant;
pub mod ringdown;
pub mod sgwb;
pub mod strain;
pub mod tides;
pub mod ultralight;

/// Hubble constant for h = 1 (100 km s⁻¹ Mpc⁻¹), s⁻¹. Multiply by the dimensionless h.
pub const H100: f64 = 1.0e5 / metrology::MPC;

/// `∫ g(f) df` from `lo` to `hi` (both > 0) by Simpson's rule on a logarithmic grid of `n`
/// intervals (`n` is made even). Accurate for the smooth power-law integrands used here.
pub fn integrate_log(lo: f64, hi: f64, n: usize, g: impl Fn(f64) -> f64) -> f64 {
    let n = (n.max(2) + 1) & !1;
    let (a, b) = (lo.ln(), hi.ln());
    let h = (b - a) / n as f64;
    let term = |k: usize| { let f = (a + k as f64 * h).exp(); g(f) * f };
    let inner: f64 = (1..n).map(|k| if k % 2 == 1 { 4.0 } else { 2.0 } * term(k)).sum();
    (term(0) + inner + term(n)) * h / 3.0
}

#[cfg(test)]
mod tests {
    use super::*;

    #[test]
    fn log_integral_is_exact_for_power_laws() {
        // ∫₁¹⁰ f⁻² df = 0.9
        assert!((integrate_log(1.0, 10.0, 200, |f| f.powi(-2)) - 0.9).abs() < 1e-9);
        assert!((H100 - 3.2408e-18).abs() < 1e-22);
    }
}