diff --git a/extras/n3lo_bench/plot_msht.py b/extras/n3lo_bench/plot_msht.py
index c922c1c03..8ddc92d2a 100644
--- a/extras/n3lo_bench/plot_msht.py
+++ b/extras/n3lo_bench/plot_msht.py
@@ -74,7 +74,7 @@ def msht_splitting(entry, x, nf, posterior):
         )
         qqns = np.array(
             [
-                msht.pqqps3a(x, ans0),
+                msht.p3nsa(x, ans0, 1, nf),
                 msht.p3nsa(x, ans1, 1, nf),
                 msht.p3nsa(x, ans2, 1, nf),
             ]
diff --git a/src/eko/evolution_operator/__init__.py b/src/eko/evolution_operator/__init__.py
index 3c04c79ef..0584ff87f 100644
--- a/src/eko/evolution_operator/__init__.py
+++ b/src/eko/evolution_operator/__init__.py
@@ -256,7 +256,7 @@ def quad_ker(
     is_threshold : boolean
         is this an intermediate threshold operator?
     n3lo_ad_variation : tuple
-        |N3LO| anomalous dimension variation ``(gg_var, gq_var, qg_var, qq_var, nsp_var, nsm_var, nsv_var)``
+        |N3LO| anomalous dimension variation ``(gg, gq, qg, qq, nsp, nsm, nsv)``
     is_polarized : boolean
         is polarized evolution ?
     is_time_like : boolean
@@ -369,7 +369,7 @@ def quad_ker_qcd(
     is_threshold : boolean
         is this an itermediate threshold operator?
     n3lo_ad_variation : tuple
-        |N3LO| anomalous dimension variation ``(gg_var, gq_var, qg_var, qq_var, nsp_var, nsm_var, nsv_var)``
+        |N3LO| anomalous dimension variation ``(gg, gq, qg, qq, nsp, nsm, nsv)``
     use_fhmruvv : bool
         if True use the Falcioni, Herzog, Moch, Ruijl, Ueda, Vermaseren, Vogt N3LO anomalous dimensions
 
@@ -956,11 +956,12 @@ def compute(self):
                 self.a_em[1],
             )
         logger.info(
-            "%s: order: (%d, %d), solution strategy: %s",
+            "%s: order: (%d, %d), solution strategy: %s, use fhmruvv: %s",
             self.log_label,
             self.order[0],
             self.order[1],
             self.config["method"],
+            self.config["use_fhmruvv"],
         )
 
         self.integrate()
diff --git a/src/eko/io/runcards.py b/src/eko/io/runcards.py
index faa57bf67..4087afbac 100644
--- a/src/eko/io/runcards.py
+++ b/src/eko/io/runcards.py
@@ -46,7 +46,7 @@ class TheoryCard(DictLike):
     xif: float
     """Ratio between factorization scale and process scale."""
     n3lo_ad_variation: N3LOAdVariation
-    """|N3LO| anomalous dimension variation: ``(gg_var, gq_var, qg_var, qq_var, nsp_var, nsm_var, nsv_var)``."""
+    """|N3LO| anomalous dimension variation: ``(gg, gq, qg, qq, nsp, nsm, nsv)``."""
     use_fhmruvv: Optional[bool] = False
     """If True use the Falcioni, Herzog, Moch, Ruijl, Ueda, Vermaseren, Vogt N3LO anomalous dimensions"""
     matching_order: Optional[Order] = None
diff --git a/src/eko/runner/legacy.py b/src/eko/runner/legacy.py
index e41b774a3..39659f121 100644
--- a/src/eko/runner/legacy.py
+++ b/src/eko/runner/legacy.py
@@ -78,9 +78,7 @@ def __init__(
             interpol_dispatcher=bfd,
             n3lo_ad_variation=new_theory.n3lo_ad_variation,
             use_fhmruvv=new_theory.use_fhmruvv,
-            matching_order=new_theory.matching_order
-            if new_theory.matching_order is not None
-            else new_theory.order,
+            matching_order=new_theory.matching_order,
         )
 
         with EKO.create(path) as builder:
diff --git a/src/ekore/anomalous_dimensions/unpolarized/space_like/__init__.py b/src/ekore/anomalous_dimensions/unpolarized/space_like/__init__.py
index 3d3bc7723..932b7e985 100644
--- a/src/ekore/anomalous_dimensions/unpolarized/space_like/__init__.py
+++ b/src/ekore/anomalous_dimensions/unpolarized/space_like/__init__.py
@@ -39,7 +39,7 @@ def gamma_ns(order, mode, n, nf, n3lo_ad_variation, use_fhmruvv=False):
     nf : int
         Number of active flavors
     n3lo_ad_variation : tuple
-        |N3LO| anomalous dimension variation ``(gg_var, gq_var, qg_var, qq_var, nsp_var, nsm_var, nsv_var)``
+        |N3LO| anomalous dimension variation ``(gg, gq, qg, qq, nsp, nsm, nsv)``
     use_fhmruvv: bool
         if True use the Falcioni Herzog Moch Vogt N3LO anomalous dimensions
 
@@ -111,7 +111,7 @@ def gamma_singlet(order, n, nf, n3lo_ad_variation, use_fhmruvv=False):
     nf : int
         Number of active flavors
     n3lo_ad_variation : tuple
-        |N3LO| anomalous dimension variation ``(gg_var, gq_var, qg_var, qq_var, nsp_var, nsm_var, nsv_var)``
+        |N3LO| anomalous dimension variation ``(gg, gq, qg, qq, nsp, nsm, nsv)``
     use_fhmruvv: bool
         if True use the Falcioni Herzog Moch Vogt N3LO anomalous dimensions