Skip to content

Commit 4b55d9d

Browse files
AMFlow Refactorclaude
AMFlow Refactor
and
claude
committed
bench: triangle_1L_3mass_eps01 oracle #88
Co-Authored-By: Claude Opus 4.7 <noreply@anthropic.com>
1 parent 2a62f31 commit 4b55d9d

3 files changed

Lines changed: 226 additions & 0 deletions

‎tools/bench/triangle_1L_3mass_eps01_black_box_amflow_cpp.json‎

Lines changed: 66 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,66 @@
1+
{
2+
"mode": "black_box_amflow",
3+
"options": {
4+
"working_pre": 120,
5+
"chop_pre": 20,
6+
"rationalize_pre": 100,
7+
"silent_mode": true,
8+
"x_order": 200,
9+
"extra_x_order": 240,
10+
"learn_x_order": -1,
11+
"test_x_order": 5
12+
},
13+
"family": {
14+
"name": "tri3me1",
15+
"loops": [
16+
"l"
17+
],
18+
"legs": [
19+
"p1",
20+
"p2"
21+
],
22+
"conservation": {},
23+
"replacement": {
24+
"p1^2": "0",
25+
"p2^2": "0",
26+
"(p1 + p2)^2": "s"
27+
},
28+
"propagators": [
29+
"l^2 - mAsq",
30+
"(l + p1)^2 - mBsq",
31+
"(l + p1 + p2)^2 - mCsq"
32+
]
33+
},
34+
"integrals": [
35+
{
36+
"indices": [
37+
1,
38+
1,
39+
1
40+
]
41+
}
42+
],
43+
"eps_samples": [
44+
"1/100"
45+
],
46+
"work_dir": "/tmp/amflow_bench_triangle_1L_3mass_eps01_cpp",
47+
"amf_options": {
48+
"ending_schemes": [
49+
"Tradition",
50+
"SingleMass"
51+
],
52+
"max_recursion_depth": 16,
53+
"blackbox": {
54+
"ibp_rank": 0,
55+
"ibp_dot": 3,
56+
"integral_order": 5,
57+
"n_thread": 4,
58+
"numeric_values": {
59+
"s": "-3",
60+
"mAsq": "1",
61+
"mBsq": "4",
62+
"mCsq": "9"
63+
}
64+
}
65+
}
66+
}

‎tools/bench/triangle_1L_3mass_eps01_black_box_amflow_mma.wl‎

Lines changed: 86 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,86 @@
1+
(* ::Package:: *)
2+
3+
(*
4+
1-loop triangle 3mass at eps=1/100 (eps variation).
5+
*)
6+
7+
current = If[$FrontEnd === Null, $InputFileName, NotebookFileName[]] // DirectoryName;
8+
repo = FileNameJoin[{current, "..", ".."}];
9+
Get[FileNameJoin[{repo, "reference", "amflow-master", "AMFlow.m"}]];
10+
11+
SetReductionOptions[
12+
"IBPReducer" -> "Kira",
13+
"BlackBoxRank" -> 0,
14+
"BlackBoxDot" -> 3,
15+
"DeleteBlackBoxDirectory" -> False
16+
];
17+
Kira`$KiraExecutable = "/usr/local/bin/kira";
18+
Kira`$FermatExecutable = "/usr/share/Ferl7/fer64";
19+
20+
SetAMFOptions[
21+
"AMFMode" -> {"Prescription", "Mass", "Propagator"},
22+
"EndingScheme" -> {"Tradition", "SingleMass"},
23+
"WorkingPre" -> 120, "ChopPre" -> 20, "RationalizePre" -> 100,
24+
"XOrder" -> 200, "ExtraXOrder" -> 240,
25+
"LearnXOrder" -> -1, "TestXOrder" -> 5
26+
];
27+
28+
AMFlowInfo["Family"] = tri3me1;
29+
AMFlowInfo["Loop"] = {l};
30+
AMFlowInfo["Leg"] = {p1, p2};
31+
AMFlowInfo["Conservation"] = {};
32+
AMFlowInfo["Replacement"] = {
33+
p1^2 -> 0,
34+
p2^2 -> 0,
35+
(p1 + p2)^2 -> s
36+
};
37+
AMFlowInfo["Propagator"] = {
38+
l^2 - mAsq,
39+
(l + p1)^2 - mBsq,
40+
(l + p1 + p2)^2 - mCsq
41+
};
42+
AMFlowInfo["Numeric"] = {
43+
s -> -3,
44+
mAsq -> 1,
45+
mBsq -> 4,
46+
mCsq -> 9
47+
};
48+
AMFlowInfo["NThread"] = 4;
49+
50+
targets = { j[tri3me1, 1, 1, 1] };
51+
epslist = {1/100};
52+
cache = Environment["AMF_BENCH_MMA_CACHE"];
53+
If[cache === $Failed || cache === None,
54+
cache = FileNameJoin[{repo, "tools", "bench", "mma_refs",
55+
"triangle_1L_3mass_eps01_mma_cache"}]];
56+
57+
Print["==BENCH== mma 1-loop triangle 3mass at eps=1/100 (eps variation)."];
58+
Print["==CACHE== ", cache];
59+
Print["==EPS== ", epslist // InputForm];
60+
61+
elapsed = AbsoluteTiming[
62+
sol = BlackBoxAMFlow[targets, epslist, cache];
63+
][[1]];
64+
65+
Print["==TIME== ", elapsed];
66+
Print["==RESULT_INPUTFORM== ", sol // InputForm];
67+
68+
Do[
69+
key = Keys[sol][[i]];
70+
val = Values[sol][[i, 1]];
71+
Print["==VALUE== ", key // InputForm, " :: ",
72+
N[Re[val], 50] // InputForm, " :: ", N[Im[val], 50] // InputForm],
73+
{i, Length[Keys[sol]]}
74+
];
75+
76+
sysSol = Get[FileNameJoin[{cache, "1", "solution"}]];
77+
If[Head[sysSol] === List,
78+
Do[
79+
key = sysSol[[i, 1]];
80+
val = sysSol[[i, 2, 1]];
81+
Print["==SAMPLED== ", key // InputForm, " :: ",
82+
N[Re[val], 50] // InputForm, " :: ", N[Im[val], 50] // InputForm],
83+
{i, Length[sysSol]}
84+
]];
85+
86+
Quit[];

‎tools/bench/triangle_1L_3mass_eps01_mma_reference.json‎

Lines changed: 74 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,74 @@
1+
{
2+
"family": "tri3me1",
3+
"source": "1-loop triangle 3-mass at eps=1/100.",
4+
"eps": "1/100",
5+
"numeric": {
6+
"s": "-3",
7+
"mAsq": "1",
8+
"mBsq": "4",
9+
"mCsq": "9"
10+
},
11+
"options": {
12+
"WorkingPre": 120,
13+
"ChopPre": 20,
14+
"RationalizePre": 100,
15+
"XOrder": 240,
16+
"ExtraXOrder": 280,
17+
"LearnXOrder": -1,
18+
"TestXOrder": 5,
19+
"BlackBoxRank": 0,
20+
"BlackBoxDot": 3,
21+
"EndingScheme": [
22+
"Tradition",
23+
"SingleMass"
24+
]
25+
},
26+
"mma_black_box_amflow_time_seconds": 42.294965,
27+
"mma_cache": "tools/bench/mma_refs/triangle_1L_3mass_eps01_mma_cache",
28+
"requested_targets": [
29+
{
30+
"integral": "j[tri3me1, 1, 1, 1]",
31+
"value": {
32+
"re": "-0.11536467115323117059547042679405492696795709796478950741113050691898634035966",
33+
"im": "1.17198768532832343377722629997973751080006806117516603980439831203068825e-64"
34+
}
35+
}
36+
],
37+
"sampled_values_from_same_mma_run": [
38+
{
39+
"integral": "j[tri3me1, 1, 1, 1]",
40+
"value": {
41+
"re": "-0.11536467115323117059547042679405492696795709796478950741113050691898634035966",
42+
"im": "1.17198768532832343377722629997973751080006806117516603980439831203068825e-64"
43+
}
44+
},
45+
{
46+
"integral": "j[tri3me1, 1, 0, 0]",
47+
"value": {
48+
"re": "100.43695466580869061973029178839469025710053075296053705860416896003287937199187",
49+
"im": "-1.3158348773542797770445893981204167514902189087763771776381982506383205e-65"
50+
}
51+
},
52+
{
53+
"integral": "j[tri3me1, 0, 1, 0]",
54+
"value": {
55+
"re": "396.2168377245496567194219509235809964156100178755608433512402",
56+
"im": "3.2013486382109397472483845742544427655728317416044216329011008413431202e-65"
57+
}
58+
},
59+
{
60+
"integral": "j[tri3me1, 0, 0, 1]",
61+
"value": {
62+
"re": "884.2877736436795351054755603379866803372462062732095027274345",
63+
"im": "-3.4990662739835032252525881517093035631082236660117353775136316252894371e-65"
64+
}
65+
},
66+
{
67+
"integral": "j[tri3me1, 1, 0, 1]",
68+
"value": {
69+
"re": "97.87287032848707738890046432766073357357911114843254030597236948189243356317957",
70+
"im": "3.59344296803787307505388492926136041093609462335900691655927089168933582e-63"
71+
}
72+
}
73+
]
74+
}

0 commit comments

Comments
 (0)