rust: Hook NNLO QCD ad in

NNPDF · May 28, 2024 · 6cba59b · 6cba59b
1 parent 2891926
commit 6cba59b
Showing 1 changed file with 14 additions and 0 deletions.
diff --git a/crates/ekore/src/anomalous_dimensions/unpolarized/spacelike.rs b/crates/ekore/src/anomalous_dimensions/unpolarized/spacelike.rs
@@ -22,6 +22,16 @@ pub fn gamma_ns_qcd(order_qcd: usize, mode: u16, c: &mut Cache, nf: u8) -> Vec<C
         };
         gamma_ns[1] = gamma_ns_1
     }
+    // NNLO and beyond
+    if order_qcd >= 3 {
+        let gamma_ns_2 = match mode {
+            PID_NSP => as3::gamma_nsp(c, nf),
+            PID_NSM => as3::gamma_nsm(c, nf),
+            PID_NSV => as3::gamma_nsv(c, nf),
+            _ => panic!("Unkown non-singlet sector element"),
+        };
+        gamma_ns[2] = gamma_ns_2
+    }
     gamma_ns
 }
 
@@ -39,5 +49,9 @@ pub fn gamma_singlet_qcd(order_qcd: usize, c: &mut Cache, nf: u8) -> Vec<[[Compl
     if order_qcd >= 2 {
         gamma_S[1] = as2::gamma_singlet(c, nf);
     }
+    // NNLO and beyond
+    if order_qcd >= 3 {
+        gamma_S[2] = as3::gamma_singlet(c, nf);
+    }
     gamma_S
 }