Fix apfel benchmarks

src/ekomark/benchmark/external/apfel_utils.py

@@ -72,8 +72,8 @@ def compute_apfel_data(
 
     # Run
     apf_tabs = {}
-    for mu2 in operators["Q2grid"]:
-        apfel.EvolveAPFEL(theory["Q0"], np.sqrt(mu2))
+    for mu in operators["mugrid"]:
+        apfel.EvolveAPFEL(theory["Q0"], np.sqrt(mu**2))
         print(f"Executing APFEL took {(time.perf_counter() - apf_start)} s")
 
         tab = {}
@@ -111,7 +111,7 @@ def compute_apfel_data(
             evol_pdf = rotate_flavor_to_evolution @ pdfs
             tab = dict(zip(evol_basis, evol_pdf))
 
-        apf_tabs[mu2] = tab
+        apf_tabs[mu**2] = tab
 
     ref = {
         "target_xgrid": target_xgrid,

src/ekomark/benchmark/external/pegasus_utils.py

@@ -82,12 +82,12 @@ def compute_pegasus_data(theory, operators, skip_pdfs, rotate_to_evolution_basis
 
     # run pegaus
     out_tabs = {}
-    for mu2 in operators["Q2grid"]:
+    for mu in operators["mugrid"]:
         tab = {}
         for x in target_xgrid:
             # last two numbers are the min and max pid to calculate,
             # keep everthing for simplicity.
-            xf, _ = pegasus.xparton(x, mu2, -6, 6)
+            xf, _ = pegasus.xparton(x, mu**2, -6, 6)
             temp = dict(zip(labels, xf))
             for pid in labels:
                 if pid in skip_pdfs:
@@ -114,7 +114,7 @@ def compute_pegasus_data(theory, operators, skip_pdfs, rotate_to_evolution_basis
             evol_pdf = rotate_flavor_to_evolution @ pdfs
             tab = dict(zip(evol_basis, evol_pdf))
 
-        out_tabs[mu2] = tab
+        out_tabs[mu**2] = tab
 
     ref = {"target_xgrid": target_xgrid, "values": out_tabs}