light-curve-dmdt 0.2.6

dm-dt maps generator
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

use conv::*;
use criterion::{black_box, Criterion};
use light_curve_dmdt::{DmDt, ErfFloat, ErrorFunction, Grid};
use ndarray::Array1;

pub fn bench_gausses<T>(c: &mut Criterion)
where
    T: ErfFloat + ValueFrom<f32>,
{
    let dmdt = DmDt {
        lgdt_grid: Grid::new(T::zero(), 2.0_f32.value_as::<T>().unwrap(), 32),
        dm_grid: Grid::new(
            (-1.25_f32).value_as::<T>().unwrap(),
            1.25_f32.value_as::<T>().unwrap(),
            32,
        ),
    };

    let t = Array1::linspace(T::zero(), 100.0_f32.value_as::<T>().unwrap(), 101);
    let m = t.mapv(T::sin);
    // err is ~0.03
    let err2 = Array1::from_elem(101, 0.001_f32.value_as::<T>().unwrap());

    let erfs = [ErrorFunction::Exact, ErrorFunction::Eps1Over1e3];

    for erf in erfs.iter() {
        c.bench_function(
            format!(
                "DmDt<{}>::convert_lc_to_gausses, erf type is {:?}",
                std::any::type_name::<T>(),
                erf
            )
            .as_str(),
            |b| {
                b.iter(|| {
                    black_box(dmdt.gausses(
                        t.as_slice_memory_order().unwrap(),
                        m.as_slice_memory_order().unwrap(),
                        err2.as_slice_memory_order().unwrap(),
                        &erf,
                    ));
                })
            },
        );
    }
}