use-nuclear 0.0.1

Radioactive decay and nuclear scalar helpers for RustUse
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

#![allow(clippy::float_cmp)]

use use_nuclear::{
    ATOMIC_MASS_UNIT_MEV_C2, DecayLaw, NuclideNumbers, activity,
    binding_energy_mev_from_mass_defect_u,
};

fn assert_approx_eq(left: f64, right: f64) {
    let tolerance = 1.0e-12 * left.abs().max(right.abs()).max(1.0);
    assert!((left - right).abs() <= tolerance);
}

#[test]
fn basic_usage_covers_decay_energy_and_counts() -> Result<(), &'static str> {
    let decay_law = DecayLaw::from_half_life(10.0).ok_or("expected valid half-life")?;
    let remaining = decay_law
        .remaining_quantity(100.0, 10.0)
        .ok_or("expected valid remaining quantity")?;
    let helium = NuclideNumbers::new(4, 2).ok_or("expected valid nuclide numbers")?;

    assert_approx_eq(remaining, 50.0);
    assert_eq!(activity(2.0, 10.0), Some(20.0));
    assert_eq!(helium.proton_count(), 2);
    assert_eq!(helium.neutron_count(), 2);
    assert_eq!(
        binding_energy_mev_from_mass_defect_u(1.0),
        Some(ATOMIC_MASS_UNIT_MEV_C2),
    );

    Ok(())
}