andam 0.1.2

A comprehensive Rust library for cosmological calculations, visualizations and simulations
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

//! High-resolution BBN evolution with equations

use andam::early_universe::*;
use andam::visualization::equation_plots::*;

fn main() -> Result<(), Box<dyn std::error::Error>> {
    println!("Computing BBN evolution...");
    println!("Using analytical freeze-out approach\n");

    // BBN parameters
    let params = BBNParameters::default();

    // Run BBN simulation
    let result = run_bbn(&params);

    // Print results
    println!("=== BBN Results ===");
    println!("Initial temperature: {:.2e} K", params.temp_initial);
    println!("Baryon-to-photon ratio η: {:.2e}", params.eta);
    println!("\nFinal Abundances:");
    println!("  Y_p (⁴He): {:.6}", result.yp());
    println!("  D/H: {:.6e}", result.dh_ratio());
    println!("  ³He/H: {:.6e}", result.he3h_ratio());
    println!("  ⁷Li/H: {:.6e}", result.li7h_ratio());

    // Create high-resolution plot
    let config = PublicationConfig {
        width: 2400,
        height: 1600,
        dpi: 300,
        ..Default::default()
    };

    plot_abundance_evolution(
        "bbn_evolution_hires.png",
        &result.evolution,
        &config,
    )?;

    println!("\nCreated: bbn_evolution_hires.png");
    Ok(())
}