unit tests

crates/ekore/src/anomalous_dimensions/unpolarized/spacelike/aem1.rs

@@ -118,4 +118,21 @@ mod tests {
         let me = gamma_ns(&mut c, NF) + gamma_phq(&mut c, NF);
         assert_approx_eq_cmplx!(f64, me, Complex::zero(), epsilon = 1e-12);
     }
+
+    #[test]
+    fn photon_momentum_conservation() {
+        const N: Complex<f64> = cmplx!(2., 0.);
+        let mut c = Cache::new(N);
+
+        for nf in 2..7 {
+            let nu = uplike_flavors(nf);
+            let nd = nf - nu;
+            assert_approx_eq_cmplx!(
+                f64,
+                eu2 * gamma_qph(&mut c, nu) + ed2 * gamma_qph(&mut c, nd) + gamma_phph(&mut c, nf),
+                cmplx!(0., 0.),
+                epsilon = 2e-6
+            );
+        }
+    }
 }