ekore 0.15.4

EKO expressions
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

//! The unpolarized, space-like |OME| at various couplings power.
use crate::harmonics::cache::Cache;
use num::complex::Complex;
use num::Zero;
pub mod as1;
pub mod as2;

/// Compute the tower of the singlet |OME|.
pub fn A_singlet(
    matching_order_qcd: usize,
    c: &mut Cache,
    nf: u8,
    L: f64,
) -> Vec<[[Complex<f64>; 3]; 3]> {
    let mut A_s = vec![
        [
            [
                Complex::<f64>::zero(),
                Complex::<f64>::zero(),
                Complex::<f64>::zero()
            ],
            [
                Complex::<f64>::zero(),
                Complex::<f64>::zero(),
                Complex::<f64>::zero()
            ],
            [
                Complex::<f64>::zero(),
                Complex::<f64>::zero(),
                Complex::<f64>::zero()
            ]
        ];
        matching_order_qcd
    ];
    if matching_order_qcd >= 1 {
        A_s[0] = as1::A_singlet(c, nf, L);
    }
    if matching_order_qcd >= 2 {
        // TODO recover MSbar mass
        A_s[1] = as2::A_singlet(c, nf, L, false);
    }
    A_s
}

/// Compute the tower of the non-singlet |OME|.
pub fn A_non_singlet(
    matching_order_qcd: usize,
    c: &mut Cache,
    nf: u8,
    L: f64,
) -> Vec<[[Complex<f64>; 2]; 2]> {
    let mut A_ns = vec![
        [
            [Complex::<f64>::zero(), Complex::<f64>::zero()],
            [Complex::<f64>::zero(), Complex::<f64>::zero()]
        ];
        matching_order_qcd
    ];
    if matching_order_qcd >= 1 {
        A_ns[0] = as1::A_ns(c, nf, L);
    }
    if matching_order_qcd >= 2 {
        A_ns[1] = as2::A_ns(c, nf, L);
    }
    A_ns
}