From 6cba59b5b5b7ee40394160387356cd906218c36a Mon Sep 17 00:00:00 2001
From: Felix Hekhorn <felix.hekhorn@uni-tuebingen.de>
Date: Tue, 28 May 2024 11:07:00 +0300
Subject: [PATCH] rust: Hook NNLO QCD ad in

---
 .../anomalous_dimensions/unpolarized/spacelike.rs  | 14 ++++++++++++++
 1 file changed, 14 insertions(+)

diff --git a/crates/ekore/src/anomalous_dimensions/unpolarized/spacelike.rs b/crates/ekore/src/anomalous_dimensions/unpolarized/spacelike.rs
index e1bfafb58..dec73bf04 100644
--- 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
 }