darkmatter 0.0.2

Dark matter simulation engine — gravitational fields, particle dynamics, halo stability, and cosmological constants
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
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137

use darkmatter::candidates::fuzzy_dm::{
    FuzzyDarkMatter, condensation_fraction, core_halo_mass_relation_solar, dynamical_heating_rate,
    granule_mass_solar, gravitational_pressure, interference_fringe_scale_kpc, occupation_number,
    quantum_jeans_wavenumber, quantum_pressure, soliton_central_density, soliton_density_profile,
    soliton_mass_solar, superradiance_rate, suppression_transfer_function,
};

fn approx_eq(a: f64, b: f64, rel_tol: f64) -> bool {
    if a == 0.0 && b == 0.0 {
        return true;
    }
    let denom = a.abs().max(b.abs());
    (a - b).abs() / denom < rel_tol
}

#[test]
fn struct_basic() {
    let fdm = FuzzyDarkMatter { mass_ev: 1e-22 };
    assert!(fdm.mass_kg() > 0.0);
    assert!(approx_eq(fdm.m22(), 1.0, 1e-6));
    assert!(fdm.compton_wavelength_m() > 0.0);
}

#[test]
fn de_broglie() {
    let fdm = FuzzyDarkMatter { mass_ev: 1e-22 };
    let lambda = fdm.de_broglie_wavelength_kpc(200.0);
    assert!(lambda > 0.0);
    assert!(lambda < 100.0);
}

#[test]
fn jeans_scales() {
    let fdm = FuzzyDarkMatter { mass_ev: 1e-22 };
    let rho = 1e-25;
    let l_j = fdm.jeans_length_kpc(rho);
    let m_j = fdm.jeans_mass_solar(rho);
    assert!(l_j > 0.0);
    assert!(m_j > 0.0);
}

#[test]
fn soliton() {
    let fdm = FuzzyDarkMatter { mass_ev: 1e-22 };
    let r_c = fdm.soliton_core_radius_kpc(1e12);
    assert!(r_c > 0.0);
    assert!(r_c < 10.0);

    let rho_c = soliton_central_density(1e-22, r_c);
    assert!(rho_c > 0.0);

    let rho_r = soliton_density_profile(0.5, rho_c, r_c);
    assert!(rho_r > 0.0);
    assert!(rho_r <= rho_c);

    let m_sol = soliton_mass_solar(rho_c, r_c);
    assert!(m_sol > 0.0);
}

#[test]
fn core_halo_relation() {
    let m_core = core_halo_mass_relation_solar(1e-22, 1e12);
    assert!(m_core > 0.0);
    assert!(m_core < 1e12);
}

#[test]
fn quantum_vs_gravitational_pressure() {
    let mass_ev = 1e-22;
    let rho = 1e-25;
    let r = 3.0857e19;
    let p_q = quantum_pressure(mass_ev, rho, r);
    let p_g = gravitational_pressure(rho, r);
    assert!(p_q > 0.0);
    assert!(p_g > 0.0);
}

#[test]
fn jeans_wavenumber() {
    let k_j = quantum_jeans_wavenumber(1e-22, 1e-25);
    assert!(k_j > 0.0);
}

#[test]
fn suppression() {
    let t1 = suppression_transfer_function(0.01, 1e-22);
    let t2 = suppression_transfer_function(100.0, 1e-22);
    assert!(t1 < t2);
    assert!(t2 > 0.0);
}

#[test]
fn interference_fringe() {
    let scale = interference_fringe_scale_kpc(1e-22, 200.0);
    assert!(scale > 0.0);
}

#[test]
fn granule_mass() {
    let m_g = granule_mass_solar(1e-22, 1e-25, 200.0);
    assert!(m_g > 0.0);
}

#[test]
fn dynamical_heating() {
    let rate = dynamical_heating_rate(1e-22, 1e-25, 200.0);
    assert!(rate > 0.0);
}

#[test]
fn superradiance() {
    let rate = superradiance_rate(1e-12, 10.0, 0.9);
    assert!(rate >= 0.0);
}

#[test]
fn occupation_number_positive() {
    let n_occ = occupation_number(1e-25, 1e-22, 200.0);
    assert!(n_occ > 0.0);
}

#[test]
fn condensation() {
    let f = condensation_fraction(1.0, 10.0);
    assert!(f > 0.0);
    assert!(f < 1.0);

    let f_above = condensation_fraction(20.0, 10.0);
    assert_eq!(f_above, 0.0);
}

#[test]
fn half_mode_mass() {
    let fdm = FuzzyDarkMatter { mass_ev: 1e-22 };
    let m_hm = fdm.half_mode_mass_solar();
    assert!(m_hm > 0.0);
}