Rustify QED

benchmarks/apfel_bench.py

@@ -37,11 +37,7 @@ class BenchmarkVFNS(ApfelBenchmark):
 
     vfns_theory = {
         "FNS": "ZM-VFNS",
-        "ModEv": [
-            "EXA",
-            "EXP",
-            "TRN",
-        ],
+        "ModEv": "EXA",
         "kcThr": 1.0,
         "kbThr": 1.0,
         "ktThr": 1.0,
@@ -132,12 +128,10 @@ class BenchmarkFFNS(ApfelBenchmark):
 
     ffns_theory = {
         "FNS": "FFNS",
-        "ModEv": [
-            "EXA",
-            "EXP",
-            "TRN",
-        ],
+        "ModEv": "EXA",
         "NfFF": 4,
+        "nf0": 4,
+        "nfref": 4,
         "kcThr": 0.0,
         "kbThr": np.inf,
         "ktThr": np.inf,
@@ -159,7 +153,6 @@ def benchmark_plain_pol(self, pto):
 
         th = self.ffns_theory.copy()
         th.update({"PTO": [pto]})
-        th["ModEv"] = ["EXA"]  # TODO for the time one is sufficient
         op = operators.apfel_config.copy()
         op["polarized"] = [True]
         self.run(cartesian_product(th), operators.build(op), ["ToyLH"])
@@ -209,6 +202,8 @@ class BenchmarkFFNS_qed(ApfelBenchmark):
             # "TRN",
         ],
         "NfFF": 5,
+        "nf0": 5,
+        "nfref": 5,
         "kcThr": 0.0,
         "kbThr": 0.0,
         "ktThr": np.inf,
@@ -291,6 +286,8 @@ class BenchmarkVFNS_qed(ApfelBenchmark):
         "kbThr": 1.0,
         "ktThr": 1.0,
         "Q0": 1.25,
+        "nf0": 3,
+        "nfref": 5,
         "alphas": 0.118000,
         "alphaqed": 0.007496,
     }
@@ -319,7 +316,7 @@ def benchmark_plain(self, pto, qed):
     # obj = BenchmarkFFNS()
 
     # obj.benchmark_plain_pol(2)
-    # obj.benchmark_plain(2)
+    # obj.benchmark_plain(0)
     # obj.benchmark_sv(2, "exponentiated")
     # obj.benchmark_kthr(2)
     # obj.benchmark_msbar(2)

benchmarks/sandbox.py

@@ -44,20 +44,20 @@ class Sandbox(Runner):
     sandbox = True
 
     # select here the external program between LHA, LHAPDF, apfel, pegasus
-    external = "apfel"
+    external = "void"
     # external = "pegasus"
 
     # select to plot operators
-    plot_operator = True
+    # plot_operator = True
 
-    rotate_to_evolution_basis = True
+    # rotate_to_evolution_basis = True
 
     @staticmethod
     def generate_operators():
         ops = {
             "ev_op_iterations": [1],
             # "ev_op_max_order": [20],
-            "mugrid": [[10]],
+            "mugrid": [[10.0]],
             # "debug_skip_singlet": [True],
         }
         return ops
@@ -118,12 +118,12 @@ def lha(self):
             "FNS": "FFNS",
             "NfFF": 4,
             "ModEv": "EXA",
-            "XIF": np.sqrt(2),
-            "Q0": np.sqrt(2),
+            "XIF": float(np.sqrt(2)),
+            "Q0": float(np.sqrt(2)),
             "kcThr": 0.0,
             "kbThr": np.inf,
             "ktThr": np.inf,
-            "Qref": np.sqrt(2.0),
+            "Qref": float(np.sqrt(2.0)),
             "alphas": 0.35,
         }
         self.skip_pdfs = lambda _theory: [
@@ -139,11 +139,11 @@ def lha(self):
         ]
         self.run(
             [theory_updates],
-            [{"mugrid": [100], "debug_skip_singlet": True}],
+            [{"mugrid": [100.0], "debug_skip_singlet": True}],
             ["ToyLH"],
         )
 
 
 if __name__ == "__main__":
     sand = Sandbox()
-    sand.doit()
+    sand.lha()