starting ome

crates/ekore/src/lib.rs

@@ -7,4 +7,5 @@ pub mod anomalous_dimensions;
 pub mod bib;
 mod constants;
 pub mod harmonics;
+pub mod operator_matrix_elements;
 pub mod util;

crates/ekore/src/operator_matrix_elements.rs

@@ -0,0 +1,3 @@
+//! |OME|.
+
+pub mod unpolarized;

crates/ekore/src/operator_matrix_elements/unpolarized.rs

@@ -0,0 +1,3 @@
+//! The unpolarized |OME| for space-like and time-like kinematics.
+
+pub mod spacelike;

crates/ekore/src/operator_matrix_elements/unpolarized/spacelike.rs

@@ -0,0 +1,12 @@
+//! The unpolarized, space-like |OME| at various couplings power.
+pub mod as1;
+
+// Compute the tower of the singlet |OME|.
+// pub fn A_singlet() -> Vec<Complex<f64>>{
+
+// }
+
+// Compute the tower of the non-singlet |OME|.
+// pub fn A_non_singlet() -> Vec<Complex<f64>>{
+
+// }

crates/ekore/src/operator_matrix_elements/unpolarized/spacelike/as1.rs

@@ -0,0 +1,25 @@
+//! NLO QCD
+
+use num::complex::Complex;
+
+use crate::constants::{CA, CF, TR};
+use crate::harmonics::cache::{Cache, K};
+
+/// Compute heavy-heavy |OME| :math:`A_{HH}^{(1)}`
+///
+/// Defined in Eq. () of .
+pub fn A_hh(c: &mut Cache, _nf: u8, L: f64) -> Complex<f64> {
+    let N = c.n;
+    let S1m = c.get(K::S1) - 1. / N;
+    let S2m = c.get(K::S2) - 1. / N.powu(2);
+    let ahh_l = (2. + N - 3. * N.powu(2)) / (N * (1. + N)) + 4. * S1m;
+    let ahh = 2.
+        * (2. + 5. * N + N.powu(2)
+            - 6. * N.powu(3)
+            - 2. * N.powu(4)
+            - 2. * N * (-1. - 2. * N + N.powu(3)) * S1m)
+        / (N * (1. + N)).powu(2)
+        + 4. * (S1m.powu(2) + S2m);
+
+    -CF * (ahh_l * L + ahh)
+}