astroimsim-data 0.1.0

Data package for astroimsim ecosystem
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

use rand::distr::Uniform;
use rand::prelude::IteratorRandom;
use rand_distr::{Binomial, Distribution, Poisson};


pub fn test_binom(){
    let num = 4000*4000*9;
    let mut vals = Vec::with_capacity(num);

    let max_flux = 2000*400*2000*60*60*24*7;
    let min_flux = max_flux/2;
    println!("Binomial Distribution");
    println!("Assuming photon count is random and between {min_flux} and {max_flux} photons...");
    let now = std::time::Instant::now();
    let photons = Uniform::new(min_flux,max_flux).unwrap();
    let mut rng = rand::rng();
    let mut rng2 = rand::rng();
    for i in 0..num{
        let flux =photons.sample(&mut rng);
        let bin = Binomial::new(flux, 0.3).unwrap();
        let v = bin.sample(&mut rng2);
        vals.push(v)
    }

    let sum:u64 = vals.iter().sum();
    let duration = now.elapsed().as_secs();
    println!("Calculated values for 4k*4k*9 = {:?} in {:?}s, sum is {:?}",num,duration, sum);



    let num = 4000*4000*9;
    let mut vals = Vec::with_capacity(num);

    let max_flux:u64 = 2000*400*2000*60*60*24*7;
    let min_flux = max_flux/2;

    println!("Poisson Distribution");
    println!("Assuming photon count is random and between {min_flux} and {max_flux} photons...");
    let now = std::time::Instant::now();
    let photons = Uniform::new(min_flux,max_flux).unwrap();
    let mut rng = rand::rng();
    let mut rng2 = rand::rng();
    for i in 0..num{
        let flux =photons.sample(&mut rng);
        let bin = Poisson::new(flux as f64).unwrap();
        let v = bin.sample(&mut rng2);
        vals.push(v)
    }

    let sum:f64 = vals.iter().sum();
    let duration = now.elapsed().as_secs();
    println!("Calculated values for 4k*4k*9 = {:?} in {:?}s, sum is {:?}",num,duration, sum);




}