Fix apfel benchmarks

src/eko/io/runcards.py

@@ -309,7 +309,7 @@ def flavored_mugrid(mugrid: list, masses: list, matching_ratios: list):
 
     """
     atlas = default_atlas(masses, matching_ratios)
-    return [(mu, nf_default(mu**2, atlas)) for mu in mugrid]
+    return [(float(mu), nf_default(mu**2, atlas)) for mu in mugrid]
 
 
 # TODO: move to a more suitable place

src/ekomark/benchmark/external/apfel_utils.py

@@ -72,8 +72,8 @@ def compute_apfel_data(
 
     # Run
     apf_tabs = {}
-    for mu in operators["mugrid"]:
-        apfel.EvolveAPFEL(theory["Q0"], mu)
+    for mu2 in operators["Q2grid"]:
+        apfel.EvolveAPFEL(theory["Q0"], np.sqrt(mu2))
         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[mu**2] = tab
+        apf_tabs[mu2] = tab
 
     ref = {
         "target_xgrid": target_xgrid,

src/ekomark/benchmark/external/pegasus_utils.py

@@ -82,7 +82,7 @@ def compute_pegasus_data(theory, operators, skip_pdfs, rotate_to_evolution_basis
 
     # run pegaus
     out_tabs = {}
-    for mu2 in np.array(operators["mugrid"]) ** 2:
+    for mu2 in operators["Q2grid"]:
         tab = {}
         for x in target_xgrid:
             # last two numbers are the min and max pid to calculate,

src/ekomark/benchmark/runner.py

@@ -4,6 +4,7 @@
 import os
 import sys
 
+import numpy as np
 import pandas as pd
 from banana import cfg as banana_cfg
 from banana.benchmark.runner import BenchmarkRunner
@@ -245,10 +246,11 @@ def log(self, theory, _, pdf, me, ext):
                 flavor_rotation=rotate_to_evolution,
                 qed=qed,
             )
-        for q2 in q2s:
+        for (q2, ref_pdfs), (q2a, res) in zip(ext["values"].items(), pdf_grid.items()):
+            np.testing.assert_allclose(q2, q2a)
             log_tab = dfdict.DFdict()
-            ref_pdfs = ext["values"][q2]
-            res = pdf_grid[q2]
+            # ref_pdfs = ext["values"][q2]
+            # res = pdf_grid[q2]
             my_pdfs = res["pdfs"]
             my_pdf_errs = res["errors"]