FCLoopFindTopologyMappings bug

[ygtw1]

The programs automatically quit when running findMappings in FCLoopFindTopologyMappings. My programs are as follows.

<< FeynCalc`;

FeynCalc 10.1.0 (dev version). For help, use the online documentation, visit the forum and have a look at the supplied examples. The PDF-version of the manual can be downloaded here.

If you use FeynCalc in your research, please evaluate FeynCalcHowToCite[] to learn how to cite this software.

Please keep in mind that the proper academic attribution of our work is crucial to ensure the future development of this package!

LaunchKernels[];
$ParallelizeFeynCalc = True;

SPD[P] = 4;
SPD[k1] = 0;
SPD[P, k1] = 2;

topos1 = {FCTopology[
"qtopo1", {SFAD[{{q3, 0}, {0, 1}, 1}], SFAD[{{q2, 0}, {0, 1}, 1}],
SFAD[{{q1, 0}, {0, 1}, 1}], SFAD[{{q2 + q3, 0}, {0, 1}, 1}],
SFAD[{{q1 - q2, 0}, {0, 1}, 1}],
SFAD[{{P/2 + q3, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q2, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q1, 0}, {1, 1}, 1}],
SFAD[{{k1 - P/2 + q3, 0}, {1, 1}, 1}],
SFAD[{{k1 - P/2 - q2, 0}, {1, 1}, 1}]}, {q1, q2, q3}, {k1,
P}, {Hold[Pair][Momentum[P, D], Momentum[P, D]] -> 4,
Hold[Pair][Momentum[k1, D], Momentum[k1, D]] -> 0,
Hold[Pair][Momentum[k1, D], Momentum[P, D]] -> 2}, {}],
FCTopology[
"qtopo10", {SFAD[{{q3, 0}, {1, 1}, 1}],
SFAD[{{q2, 0}, {0, 1}, 1}], SFAD[{{q1, 0}, {0, 1}, 1}],
SFAD[{{q2 + q3, 0}, {1, 1}, 1}], SFAD[{{q1 - q2, 0}, {0, 1}, 1}],
SFAD[{{P/2 + q3, 0}, {0, 1}, 1}],
SFAD[{{P/2 + q2, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q1, 0}, {1, 1}, 1}],
SFAD[{{k1 - P/2 - q2, 0}, {1, 1}, 1}],
SFAD[{{k1 - P - q2 - q3, 0}, {1, 1}, 1}]}, {q1, q2, q3}, {k1,
P}, {Hold[Pair][Momentum[P, D], Momentum[P, D]] -> 4,
Hold[Pair][Momentum[k1, D], Momentum[k1, D]] -> 0,
Hold[Pair][Momentum[k1, D], Momentum[P, D]] -> 2}, {}],
FCTopology[
"qtopo1006", {SFAD[{{q3, 0}, {1, 1}, 1}],
SFAD[{{q2, 0}, {1, 1}, 1}], SFAD[{{q1, 0}, {1, 1}, 1}],
SFAD[{{q1 + q2, 0}, {0, 1}, 1}],
SFAD[{{q1 + q2 + q3, 0}, {1, 1}, 1}],
SFAD[{{k1 + q1 + q2 + q3, 0}, {1, 1}, 1}],
SFAD[{{k1 - P + q2, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q3, 0}, {0, 1}, 1}],
SFAD[{{P/2 + q1, 0}, {0, 1}, 1}],
SFAD[{{k1 - P/2 + q1 + q2, 0}, {1, 1}, 1}]}, {q1, q2, q3}, {k1,
P}, {Hold[Pair][Momentum[P, D], Momentum[P, D]] -> 4,
Hold[Pair][Momentum[k1, D], Momentum[k1, D]] -> 0,
Hold[Pair][Momentum[k1, D], Momentum[P, D]] -> 2}, {}],
FCTopology[
"qtopo1013", {SFAD[{{q3, 0}, {1, 1}, 1}],
SFAD[{{q2, 0}, {1, 1}, 1}], SFAD[{{q1, 0}, {1, 1}, 1}],
SFAD[{{q1 + q2, 0}, {0, 1}, 1}],
SFAD[{{q1 + q2 + q3, 0}, {1, 1}, 1}],
SFAD[{{k1 - P + q3, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q2, 0}, {0, 1}, 1}],
SFAD[{{P/2 + q1, 0}, {0, 1}, 1}],
SFAD[{{k1 - P/2 + q2 + q3, 0}, {1, 1}, 1}],
SFAD[{{P/2 - q2 - q3, 0}, {1, 1}, 1}]}, {q1, q2, q3}, {k1,
P}, {Hold[Pair][Momentum[P, D], Momentum[P, D]] -> 4,
Hold[Pair][Momentum[k1, D], Momentum[k1, D]] -> 0,
Hold[Pair][Momentum[k1, D], Momentum[P, D]] -> 2}, {}],
FCTopology[
"qtopo1021", {SFAD[{{q3, 0}, {0, 1}, 1}],
SFAD[{{q2, 0}, {0, 1}, 1}], SFAD[{{q1, 0}, {1, 1}, 1}],
SFAD[{{q1 + q2, 0}, {1, 1}, 1}], SFAD[{{k1 - q3, 0}, {0, 1}, 1}],
SFAD[{{P/2 + q2, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q1, 0}, {0, 1}, 1}],
SFAD[{{k1 - P/2 + q2, 0}, {1, 1}, 1}],
SFAD[{{k1 - P/2 - q1, 0}, {0, 1}, 1}],
SFAD[{{k1 - P/2 - q1 - q3, 0}, {0, 1}, 1}]}, {q1, q2, q3}, {k1,
P}, {Hold[Pair][Momentum[P, D], Momentum[P, D]] -> 4,
Hold[Pair][Momentum[k1, D], Momentum[k1, D]] -> 0,
Hold[Pair][Momentum[k1, D], Momentum[P, D]] -> 2}, {}],
FCTopology[
"qtopo1030", {SFAD[{{q3, 0}, {1, 1}, 1}],
SFAD[{{q2, 0}, {1, 1}, 1}], SFAD[{{q1, 0}, {1, 1}, 1}],
SFAD[{{q1 + q2, 0}, {0, 1}, 1}],
SFAD[{{q1 + q2 + q3, 0}, {1, 1}, 1}],
SFAD[{{P - q3, 0}, {1, 1}, 1}],
SFAD[{{k1 - P + q3, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q2, 0}, {0, 1}, 1}],
SFAD[{{P/2 + q1, 0}, {0, 1}, 1}],
SFAD[{{P/2 - q2 - q3, 0}, {1, 1}, 1}]}, {q1, q2, q3}, {k1,
P}, {Hold[Pair][Momentum[P, D], Momentum[P, D]] -> 4,
Hold[Pair][Momentum[k1, D], Momentum[k1, D]] -> 0,
Hold[Pair][Momentum[k1, D], Momentum[P, D]] -> 2}, {}],
FCTopology[
"qtopo104", {SFAD[{{q3, 0}, {0, 1}, 1}],
SFAD[{{q2, 0}, {0, 1}, 1}], SFAD[{{q1, 0}, {1, 1}, 1}],
SFAD[{{q2 + q3, 0}, {0, 1}, 1}], SFAD[{{k1 - q1, 0}, {1, 1}, 1}],
SFAD[{{k1 - P + q3, 0}, {0, 1}, 1}],
SFAD[{{k1 - P - q2, 0}, {0, 1}, 1}],
SFAD[{{P/2 - q1, 0}, {0, 1}, 1}],
SFAD[{{P/2 - q1 + q2, 0}, {0, 1}, 1}],
SFAD[{{k1 - P/2 - q2, 0}, {1, 1}, 1}]}, {q1, q2, q3}, {k1,
P}, {Hold[Pair][Momentum[P, D], Momentum[P, D]] -> 4,
Hold[Pair][Momentum[k1, D], Momentum[k1, D]] -> 0,
Hold[Pair][Momentum[k1, D], Momentum[P, D]] -> 2}, {}],
FCTopology[
"qtopo105", {SFAD[{{q3, 0}, {1, 1}, 1}],
SFAD[{{q2, 0}, {0, 1}, 1}], SFAD[{{q1, 0}, {0, 1}, 1}],
SFAD[{{q1 + q2, 0}, {0, 1}, 1}],
SFAD[{{q1 + q2 + q3, 0}, {1, 1}, 1}],
SFAD[{{k1 - q3, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q2, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q1, 0}, {1, 1}, 1}],
SFAD[{{k1 - P/2 - q1 - q3, 0}, {0, 1}, 1}],
SFAD[{{k1 - P - q1 - q2 - q3, 0}, {1, 1}, 1}]}, {q1, q2,
q3}, {k1, P}, {Hold[Pair][Momentum[P, D], Momentum[P, D]] -> 4,
Hold[Pair][Momentum[k1, D], Momentum[k1, D]] -> 0,
Hold[Pair][Momentum[k1, D], Momentum[P, D]] -> 2}, {}],
FCTopology[
"qtopo1066", {SFAD[{{q3, 0}, {1, 1}, 1}],
SFAD[{{q2, 0}, {1, 1}, 1}], SFAD[{{q1, 0}, {1, 1}, 1}],
SFAD[{{q1 + q2, 0}, {0, 1}, 1}],
SFAD[{{q1 + q2 + q3, 0}, {1, 1}, 1}],
SFAD[{{k1 - q3, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q2, 0}, {0, 1}, 1}],
SFAD[{{P/2 + q1, 0}, {0, 1}, 1}],
SFAD[{{k1 - P/2 - q2 - q3, 0}, {1, 1}, 1}],
SFAD[{{k1 - P - q1 - q2 - q3, 0}, {1, 1}, 1}]}, {q1, q2,
q3}, {k1, P}, {Hold[Pair][Momentum[P, D], Momentum[P, D]] -> 4,
Hold[Pair][Momentum[k1, D], Momentum[k1, D]] -> 0,
Hold[Pair][Momentum[k1, D], Momentum[P, D]] -> 2}, {}],
FCTopology[
"qtopo1085", {SFAD[{{q3, 0}, {1, 1}, 1}],
SFAD[{{q2, 0}, {1, 1}, 1}], SFAD[{{q1, 0}, {1, 1}, 1}],
SFAD[{{q1 + q2, 0}, {0, 1}, 1}],
SFAD[{{q1 + q2 + q3, 0}, {1, 1}, 1}],
SFAD[{{k1 - q3, 0}, {1, 1}, 1}],
SFAD[{{k1 - P + q2, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q1, 0}, {0, 1}, 1}],
SFAD[{{k1 - P/2 - q3, 0}, {0, 1}, 1}],
SFAD[{{P/2 - q2 - q3, 0}, {1, 1}, 1}]}, {q1, q2, q3}, {k1,
P}, {Hold[Pair][Momentum[P, D], Momentum[P, D]] -> 4,
Hold[Pair][Momentum[k1, D], Momentum[k1, D]] -> 0,
Hold[Pair][Momentum[k1, D], Momentum[P, D]] -> 2}, {}],
FCTopology[
"qtopo1092", {SFAD[{{q3, 0}, {0, 1}, 1}],
SFAD[{{q2, 0}, {0, 1}, 1}], SFAD[{{q1, 0}, {1, 1}, 1}],
SFAD[{{q1 + q2, 0}, {1, 1}, 1}],
SFAD[{{k1 + q1 + q2, 0}, {1, 1}, 1}],
SFAD[{{P + q1 + q2 + q3, 0}, {1, 1}, 1}],
SFAD[{{k1 + q1 + q2 + q3, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q2, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q1, 0}, {0, 1}, 1}],
SFAD[{{P/2 + q1 + q3, 0}, {0, 1}, 1}]}, {q1, q2, q3}, {k1,
P}, {Hold[Pair][Momentum[P, D], Momentum[P, D]] -> 4,
Hold[Pair][Momentum[k1, D], Momentum[k1, D]] -> 0,
Hold[Pair][Momentum[k1, D], Momentum[P, D]] -> 2}, {}],
FCTopology[
"qtopo1093", {SFAD[{{q3, 0}, {1, 1}, 1}],
SFAD[{{q2, 0}, {1, 1}, 1}], SFAD[{{q1, 0}, {1, 1}, 1}],
SFAD[{{q1 + q2, 0}, {0, 1}, 1}],
SFAD[{{q1 + q2 + q3, 0}, {1, 1}, 1}],
SFAD[{{k1 - q3, 0}, {1, 1}, 1}],
SFAD[{{k1 - P + q2, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q1, 0}, {0, 1}, 1}],
SFAD[{{k1 - P/2 + q1 + q2, 0}, {1, 1}, 1}],
SFAD[{{k1 - P/2 - q3, 0}, {0, 1}, 1}]}, {q1, q2, q3}, {k1,
P}, {Hold[Pair][Momentum[P, D], Momentum[P, D]] -> 4,
Hold[Pair][Momentum[k1, D], Momentum[k1, D]] -> 0,
Hold[Pair][Momentum[k1, D], Momentum[P, D]] -> 2}, {}],
FCTopology[
"qtopo1100", {SFAD[{{q3, 0}, {1, 1}, 1}],
SFAD[{{q2, 0}, {1, 1}, 1}], SFAD[{{q1, 0}, {1, 1}, 1}],
SFAD[{{q1 + q2, 0}, {0, 1}, 1}],
SFAD[{{q1 + q2 + q3, 0}, {1, 1}, 1}],
SFAD[{{k1 - q3, 0}, {1, 1}, 1}],
SFAD[{{k1 - P + q2, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q1, 0}, {0, 1}, 1}],
SFAD[{{k1 - P/2 + q2, 0}, {0, 1}, 1}],
SFAD[{{P/2 - q2 - q3, 0}, {1, 1}, 1}]}, {q1, q2, q3}, {k1,
P}, {Hold[Pair][Momentum[P, D], Momentum[P, D]] -> 4,
Hold[Pair][Momentum[k1, D], Momentum[k1, D]] -> 0,
Hold[Pair][Momentum[k1, D], Momentum[P, D]] -> 2}, {}],
FCTopology[
"qtopo1105", {SFAD[{{q3, 0}, {1, 1}, 1}],
SFAD[{{q2, 0}, {1, 1}, 1}], SFAD[{{q1, 0}, {1, 1}, 1}],
SFAD[{{q1 + q2, 0}, {0, 1}, 1}],
SFAD[{{q1 + q2 + q3, 0}, {1, 1}, 1}],
SFAD[{{k1 - q3, 0}, {1, 1}, 1}],
SFAD[{{k1 - P + q2, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q1, 0}, {0, 1}, 1}],
SFAD[{{k1 - P/2 + q2, 0}, {0, 1}, 1}],
SFAD[{{k1 - P/2 - q1 - q3, 0}, {1, 1}, 1}]}, {q1, q2, q3}, {k1,
P}, {Hold[Pair][Momentum[P, D], Momentum[P, D]] -> 4,
Hold[Pair][Momentum[k1, D], Momentum[k1, D]] -> 0,
Hold[Pair][Momentum[k1, D], Momentum[P, D]] -> 2}, {}],
FCTopology[
"qtopo1111", {SFAD[{{q3, 0}, {1, 1}, 1}],
SFAD[{{q2, 0}, {1, 1}, 1}], SFAD[{{q1, 0}, {1, 1}, 1}],
SFAD[{{q1 + q2, 0}, {0, 1}, 1}],
SFAD[{{q1 + q2 + q3, 0}, {1, 1}, 1}],
SFAD[{{k1 + q1 + q2, 0}, {0, 1}, 1}],
SFAD[{{k1 - q3, 0}, {1, 1}, 1}],
SFAD[{{k1 - P + q2, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q1, 0}, {0, 1}, 1}],
SFAD[{{k1 - P/2 + q1 + q2, 0}, {1, 1}, 1}]}, {q1, q2, q3}, {k1,
P}, {Hold[Pair][Momentum[P, D], Momentum[P, D]] -> 4,
Hold[Pair][Momentum[k1, D], Momentum[k1, D]] -> 0,
Hold[Pair][Momentum[k1, D], Momentum[P, D]] -> 2}, {}],
FCTopology[
"qtopo1119", {SFAD[{{q3, 0}, {1, 1}, 1}],
SFAD[{{q2, 0}, {1, 1}, 1}], SFAD[{{q1, 0}, {1, 1}, 1}],
SFAD[{{q1 + q2, 0}, {0, 1}, 1}],
SFAD[{{q1 + q2 + q3, 0}, {1, 1}, 1}],
SFAD[{{k1 + q1 + q2, 0}, {0, 1}, 1}],
SFAD[{{k1 - q3, 0}, {1, 1}, 1}],
SFAD[{{k1 - P + q2, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q1, 0}, {0, 1}, 1}],
SFAD[{{k1 - P/2 + q2, 0}, {0, 1}, 1}]}, {q1, q2, q3}, {k1,
P}, {Hold[Pair][Momentum[P, D], Momentum[P, D]] -> 4,
Hold[Pair][Momentum[k1, D], Momentum[k1, D]] -> 0,
Hold[Pair][Momentum[k1, D], Momentum[P, D]] -> 2}, {}],
FCTopology[
"qtopo1127", {SFAD[{{q3, 0}, {1, 1}, 1}],
SFAD[{{q2, 0}, {1, 1}, 1}], SFAD[{{q1, 0}, {1, 1}, 1}],
SFAD[{{q1 + q2, 0}, {0, 1}, 1}],
SFAD[{{q1 + q2 + q3, 0}, {1, 1}, 1}],
SFAD[{{k1 - P + q2, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q3, 0}, {0, 1}, 1}],
SFAD[{{P/2 + q1, 0}, {0, 1}, 1}],
SFAD[{{k1 - P/2 + q1 + q2, 0}, {1, 1}, 1}],
SFAD[{{P/2 - q1 - q2, 0}, {1, 1}, 1}]}, {q1, q2, q3}, {k1,
P}, {Hold[Pair][Momentum[P, D], Momentum[P, D]] -> 4,
Hold[Pair][Momentum[k1, D], Momentum[k1, D]] -> 0,
Hold[Pair][Momentum[k1, D], Momentum[P, D]] -> 2}, {}],
FCTopology[
"qtopo1132", {SFAD[{{q3, 0}, {1, 1}, 1}],
SFAD[{{q2, 0}, {1, 1}, 1}], SFAD[{{q1, 0}, {1, 1}, 1}],
SFAD[{{q1 + q2, 0}, {0, 1}, 1}],
SFAD[{{q1 + q2 + q3, 0}, {1, 1}, 1}],
SFAD[{{P - q2, 0}, {1, 1}, 1}],
SFAD[{{k1 - P + q2, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q3, 0}, {0, 1}, 1}],
SFAD[{{P/2 + q1, 0}, {0, 1}, 1}],
SFAD[{{P/2 - q2 - q3, 0}, {1, 1}, 1}]}, {q1, q2, q3}, {k1,
P}, {Hold[Pair][Momentum[P, D], Momentum[P, D]] -> 4,
Hold[Pair][Momentum[k1, D], Momentum[k1, D]] -> 0,
Hold[Pair][Momentum[k1, D], Momentum[P, D]] -> 2}, {}],
FCTopology[
"qtopo1149", {SFAD[{{q3, 0}, {1, 1}, 1}],
SFAD[{{q2, 0}, {1, 1}, 1}], SFAD[{{q1, 0}, {1, 1}, 1}],
SFAD[{{q1 + q2, 0}, {0, 1}, 1}],
SFAD[{{q1 + q2 + q3, 0}, {1, 1}, 1}],
SFAD[{{k1 - q2, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q3, 0}, {0, 1}, 1}],
SFAD[{{P/2 + q1, 0}, {0, 1}, 1}],
SFAD[{{P/2 - q2 - q3, 0}, {1, 1}, 1}],
SFAD[{{k1 - P/2 - q2 - q3, 0}, {1, 1}, 1}]}, {q1, q2, q3}, {k1,
P}, {Hold[Pair][Momentum[P, D], Momentum[P, D]] -> 4,
Hold[Pair][Momentum[k1, D], Momentum[k1, D]] -> 0,
Hold[Pair][Momentum[k1, D], Momentum[P, D]] -> 2}, {}],
FCTopology[
"qtopo1155", {SFAD[{{q3, 0}, {1, 1}, 1}],
SFAD[{{q2, 0}, {1, 1}, 1}], SFAD[{{q1, 0}, {1, 1}, 1}],
SFAD[{{q1 + q2, 0}, {0, 1}, 1}],
SFAD[{{q1 + q2 + q3, 0}, {1, 1}, 1}],
SFAD[{{k1 - q2, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q3, 0}, {0, 1}, 1}],
SFAD[{{P/2 + q1, 0}, {0, 1}, 1}],
SFAD[{{k1 - P/2 - q2 - q3, 0}, {1, 1}, 1}],
SFAD[{{k1 - P - q1 - q2 - q3, 0}, {1, 1}, 1}]}, {q1, q2,
q3}, {k1, P}, {Hold[Pair][Momentum[P, D], Momentum[P, D]] -> 4,
Hold[Pair][Momentum[k1, D], Momentum[k1, D]] -> 0,
Hold[Pair][Momentum[k1, D], Momentum[P, D]] -> 2}, {}],
FCTopology[
"qtopo1158", {SFAD[{{q3, 0}, {0, 1}, 1}],
SFAD[{{q2, 0}, {0, 1}, 1}], SFAD[{{q1, 0}, {1, 1}, 1}],
SFAD[{{q1 + q2, 0}, {1, 1}, 1}],
SFAD[{{k1 + q1 + q2, 0}, {1, 1}, 1}],
SFAD[{{k1 + q1 + q2 + q3, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q2, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q1, 0}, {0, 1}, 1}],
SFAD[{{P/2 + q1 + q3, 0}, {0, 1}, 1}],
SFAD[{{k1 - P/2 + q2, 0}, {1, 1}, 1}]}, {q1, q2, q3}, {k1,
P}, {Hold[Pair][Momentum[P, D], Momentum[P, D]] -> 4,
Hold[Pair][Momentum[k1, D], Momentum[k1, D]] -> 0,
Hold[Pair][Momentum[k1, D], Momentum[P, D]] -> 2}, {}],
FCTopology[
"qtopo116", {SFAD[{{q3, 0}, {1, 1}, 1}],
SFAD[{{q2, 0}, {0, 1}, 1}], SFAD[{{q1, 0}, {0, 1}, 1}],
SFAD[{{q1 + q2, 0}, {0, 1}, 1}],
SFAD[{{q1 + q2 + q3, 0}, {1, 1}, 1}],
SFAD[{{k1 - q3, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q2, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q1, 0}, {1, 1}, 1}],
SFAD[{{k1 - P/2 + q2, 0}, {1, 1}, 1}],
SFAD[{{P/2 - q2 - q3, 0}, {0, 1}, 1}]}, {q1, q2, q3}, {k1,
P}, {Hold[Pair][Momentum[P, D], Momentum[P, D]] -> 4,
Hold[Pair][Momentum[k1, D], Momentum[k1, D]] -> 0,
Hold[Pair][Momentum[k1, D], Momentum[P, D]] -> 2}, {}],
FCTopology[
"qtopo1226", {SFAD[{{q3, 0}, {1, 1}, 1}],
SFAD[{{q2, 0}, {1, 1}, 1}], SFAD[{{q1, 0}, {1, 1}, 1}],
SFAD[{{q1 + q2, 0}, {0, 1}, 1}],
SFAD[{{q1 + q2 - q3, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q2, 0}, {0, 1}, 1}],
SFAD[{{P/2 + q1, 0}, {0, 1}, 1}],
SFAD[{{P/2 - q2 + q3, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q2 - q3, 0}, {1, 1}, 1}],
SFAD[{{k1 - P/2 - q2 + q3, 0}, {1, 1}, 1}]}, {q1, q2, q3}, {k1,
P}, {Hold[Pair][Momentum[P, D], Momentum[P, D]] -> 4,
Hold[Pair][Momentum[k1, D], Momentum[k1, D]] -> 0,
Hold[Pair][Momentum[k1, D], Momentum[P, D]] -> 2}, {}],
FCTopology[
"qtopo1239", {SFAD[{{q3, 0}, {0, 1}, 1}],
SFAD[{{q2, 0}, {0, 1}, 1}], SFAD[{{q1, 0}, {1, 1}, 1}],
SFAD[{{q1 + q2, 0}, {1, 1}, 1}],
SFAD[{{k1 + q1 + q2, 0}, {1, 1}, 1}],
SFAD[{{k1 - P + q3, 0}, {0, 1}, 1}],
SFAD[{{P/2 + q2, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q1, 0}, {0, 1}, 1}],
SFAD[{{k1 - P/2 + q1, 0}, {0, 1}, 1}],
SFAD[{{P/2 + q1 - q3, 0}, {0, 1}, 1}]}, {q1, q2, q3}, {k1,
P}, {Hold[Pair][Momentum[P, D], Momentum[P, D]] -> 4,
Hold[Pair][Momentum[k1, D], Momentum[k1, D]] -> 0,
Hold[Pair][Momentum[k1, D], Momentum[P, D]] -> 2}, {}],
FCTopology[
"qtopo1240", {SFAD[{{q3, 0}, {1, 1}, 1}],
SFAD[{{q2, 0}, {1, 1}, 1}], SFAD[{{q1, 0}, {1, 1}, 1}],
SFAD[{{q1 + q2, 0}, {0, 1}, 1}],
SFAD[{{q1 + q2 - q3, 0}, {1, 1}, 1}],
SFAD[{{k1 - P + q2, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q1, 0}, {0, 1}, 1}],
SFAD[{{k1 - P/2 + q2, 0}, {0, 1}, 1}],
SFAD[{{P/2 - q2 + q3, 0}, {1, 1}, 1}],
SFAD[{{k1 - P/2 + q2 - q3, 0}, {1, 1}, 1}]}, {q1, q2, q3}, {k1,
P}, {Hold[Pair][Momentum[P, D], Momentum[P, D]] -> 4,
Hold[Pair][Momentum[k1, D], Momentum[k1, D]] -> 0,
Hold[Pair][Momentum[k1, D], Momentum[P, D]] -> 2}, {}],
FCTopology[
"qtopo1251", {SFAD[{{q3, 0}, {1, 1}, 1}],
SFAD[{{q2, 0}, {1, 1}, 1}], SFAD[{{q1, 0}, {1, 1}, 1}],
SFAD[{{q1 + q2, 0}, {0, 1}, 1}], SFAD[{{k1 - q3, 0}, {1, 1}, 1}],
SFAD[{{k1 - P + q2, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q1, 0}, {0, 1}, 1}],
SFAD[{{k1 - P/2 + q2, 0}, {0, 1}, 1}],
SFAD[{{k1 - P/2 - q1, 0}, {0, 1}, 1}],
SFAD[{{k1 - P/2 - q1 - q3, 0}, {1, 1}, 1}]}, {q1, q2, q3}, {k1,
P}, {Hold[Pair][Momentum[P, D], Momentum[P, D]] -> 4,
Hold[Pair][Momentum[k1, D], Momentum[k1, D]] -> 0,
Hold[Pair][Momentum[k1, D], Momentum[P, D]] -> 2}, {}],
FCTopology[
"qtopo127", {SFAD[{{q3, 0}, {1, 1}, 1}],
SFAD[{{q2, 0}, {0, 1}, 1}], SFAD[{{q1, 0}, {0, 1}, 1}],
SFAD[{{q1 + q2, 0}, {0, 1}, 1}],
SFAD[{{q1 + q2 + q3, 0}, {1, 1}, 1}],
SFAD[{{k1 - q3, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q2, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q1, 0}, {1, 1}, 1}],
SFAD[{{k1 - P/2 + q2, 0}, {1, 1}, 1}],
SFAD[{{k1 - P/2 - q1 - q3, 0}, {0, 1}, 1}]}, {q1, q2, q3}, {k1,
P}, {Hold[Pair][Momentum[P, D], Momentum[P, D]] -> 4,
Hold[Pair][Momentum[k1, D], Momentum[k1, D]] -> 0,
Hold[Pair][Momentum[k1, D], Momentum[P, D]] -> 2}, {}],
FCTopology[
"qtopo1279", {SFAD[{{q3, 0}, {1, 1}, 1}],
SFAD[{{q2, 0}, {1, 1}, 1}], SFAD[{{q1, 0}, {1, 1}, 1}],
SFAD[{{k1 - q2, 0}, {1, 1}, 1}],
SFAD[{{k1 - P + q1, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q3, 0}, {0, 1}, 1}],
SFAD[{{k1 - P/2 + q1 + q3, 0}, {1, 1}, 1}],
SFAD[{{P/2 - q1, 0}, {0, 1}, 1}],
SFAD[{{P/2 - q1 - q2, 0}, {1, 1}, 1}],
SFAD[{{P/2 - q1 - q2 - q3, 0}, {0, 1}, 1}]}, {q1, q2, q3}, {k1,
P}, {Hold[Pair][Momentum[P, D], Momentum[P, D]] -> 4,
Hold[Pair][Momentum[k1, D], Momentum[k1, D]] -> 0,
Hold[Pair][Momentum[k1, D], Momentum[P, D]] -> 2}, {}],
FCTopology[
"qtopo1288", {SFAD[{{q3, 0}, {1, 1}, 1}],
SFAD[{{q2, 0}, {1, 1}, 1}], SFAD[{{q1, 0}, {1, 1}, 1}],
SFAD[{{k1 - q2, 0}, {1, 1}, 1}],
SFAD[{{k1 - P + q1, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q3, 0}, {0, 1}, 1}],
SFAD[{{P - q1 - q2, 0}, {0, 1}, 1}],
SFAD[{{P/2 - q1, 0}, {0, 1}, 1}],
SFAD[{{P/2 - q1 - q2, 0}, {1, 1}, 1}],
SFAD[{{P/2 - q1 - q2 - q3, 0}, {0, 1}, 1}]}, {q1, q2, q3}, {k1,
P}, {Hold[Pair][Momentum[P, D], Momentum[P, D]] -> 4,
Hold[Pair][Momentum[k1, D], Momentum[k1, D]] -> 0,
Hold[Pair][Momentum[k1, D], Momentum[P, D]] -> 2}, {}],
FCTopology[
"qtopo1295", {SFAD[{{q3, 0}, {1, 1}, 1}],
SFAD[{{q2, 0}, {1, 1}, 1}], SFAD[{{q1, 0}, {1, 1}, 1}],
SFAD[{{k1 - q2, 0}, {1, 1}, 1}],
SFAD[{{k1 - P + q1, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q3, 0}, {0, 1}, 1}],
SFAD[{{P/2 - q1, 0}, {0, 1}, 1}],
SFAD[{{P/2 - q1 - q2, 0}, {1, 1}, 1}],
SFAD[{{k1 - P/2 - q2 - q3, 0}, {1, 1}, 1}],
SFAD[{{P/2 - q1 - q2 - q3, 0}, {0, 1}, 1}]}, {q1, q2, q3}, {k1,
P}, {Hold[Pair][Momentum[P, D], Momentum[P, D]] -> 4,
Hold[Pair][Momentum[k1, D], Momentum[k1, D]] -> 0,
Hold[Pair][Momentum[k1, D], Momentum[P, D]] -> 2}, {}],
FCTopology[
"qtopo1303", {SFAD[{{q3, 0}, {1, 1}, 1}],
SFAD[{{q2, 0}, {1, 1}, 1}], SFAD[{{q1, 0}, {1, 1}, 1}],
SFAD[{{k1 - q2, 0}, {1, 1}, 1}],
SFAD[{{k1 - P + q1, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q3, 0}, {0, 1}, 1}],
SFAD[{{P/2 - q1, 0}, {0, 1}, 1}],
SFAD[{{k1 - P/2 - q2, 0}, {0, 1}, 1}],
SFAD[{{P/2 - q1 - q2, 0}, {1, 1}, 1}],
SFAD[{{k1 - P/2 - q2 - q3, 0}, {1, 1}, 1}]}, {q1, q2, q3}, {k1,
P}, {Hold[Pair][Momentum[P, D], Momentum[P, D]] -> 4,
Hold[Pair][Momentum[k1, D], Momentum[k1, D]] -> 0,
Hold[Pair][Momentum[k1, D], Momentum[P, D]] -> 2}, {}],
FCTopology[
"qtopo1309", {SFAD[{{q3, 0}, {1, 1}, 1}],
SFAD[{{q2, 0}, {1, 1}, 1}], SFAD[{{q1, 0}, {1, 1}, 1}],
SFAD[{{k1 - q2, 0}, {1, 1}, 1}],
SFAD[{{k1 - P + q1, 0}, {1, 1}, 1}],
SFAD[{{P/2 - q1, 0}, {0, 1}, 1}],
SFAD[{{k1 - P/2 - q2, 0}, {0, 1}, 1}],
SFAD[{{k1 - P/2 - q2 + q3, 0}, {1, 1}, 1}],
SFAD[{{P/2 - q1 - q2, 0}, {1, 1}, 1}],
SFAD[{{P/2 - q1 - q2 + q3, 0}, {0, 1}, 1}]}, {q1, q2, q3}, {k1,
P}, {Hold[Pair][Momentum[P, D], Momentum[P, D]] -> 4,
Hold[Pair][Momentum[k1, D], Momentum[k1, D]] -> 0,
Hold[Pair][Momentum[k1, D], Momentum[P, D]] -> 2}, {}],
FCTopology[
"qtopo1323", {SFAD[{{q3, 0}, {1, 1}, 1}],
SFAD[{{q2, 0}, {1, 1}, 1}], SFAD[{{q1, 0}, {1, 1}, 1}],
SFAD[{{k1 - P + q1, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q2, 0}, {0, 1}, 1}],
SFAD[{{P/2 + q2 + q3, 0}, {1, 1}, 1}],
SFAD[{{k1 - P/2 + q1 + q2, 0}, {1, 1}, 1}],
SFAD[{{P/2 - q1, 0}, {0, 1}, 1}],
SFAD[{{k1 - P/2 + q1 + q2 + q3, 0}, {0, 1}, 1}],
SFAD[{{P/2 - q1 - q2, 0}, {1, 1}, 1}]}, {q1, q2, q3}, {k1,
P}, {Hold[Pair][Momentum[P, D], Momentum[P, D]] -> 4,
Hold[Pair][Momentum[k1, D], Momentum[k1, D]] -> 0,
Hold[Pair][Momentum[k1, D], Momentum[P, D]] -> 2}, {}],
FCTopology[
"qtopo1329", {SFAD[{{q3, 0}, {1, 1}, 1}],
SFAD[{{q2, 0}, {1, 1}, 1}], SFAD[{{q1, 0}, {1, 1}, 1}],
SFAD[{{k1 - P + q1, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q2, 0}, {0, 1}, 1}],
SFAD[{{P/2 - q1, 0}, {0, 1}, 1}],
SFAD[{{P/2 + q2 - q3, 0}, {1, 1}, 1}],
SFAD[{{k1 - P/2 - q2 + q3, 0}, {1, 1}, 1}],
SFAD[{{P/2 - q1 - q2, 0}, {1, 1}, 1}],
SFAD[{{P/2 - q1 - q2 + q3, 0}, {0, 1}, 1}]}, {q1, q2, q3}, {k1,
P}, {Hold[Pair][Momentum[P, D], Momentum[P, D]] -> 4,
Hold[Pair][Momentum[k1, D], Momentum[k1, D]] -> 0,
Hold[Pair][Momentum[k1, D], Momentum[P, D]] -> 2}, {}],
FCTopology[
"qtopo1334", {SFAD[{{q3, 0}, {1, 1}, 1}],
SFAD[{{q2, 0}, {1, 1}, 1}], SFAD[{{q1, 0}, {1, 1}, 1}],
SFAD[{{k1 - q1, 0}, {1, 1}, 1}],
SFAD[{{k1 - P + q3, 0}, {1, 1}, 1}],
SFAD[{{k1 - q1 + q2, 0}, {0, 1}, 1}],
SFAD[{{k1 - q1 + q2 + q3, 0}, {1, 1}, 1}],
SFAD[{{k1 - P/2 + q3, 0}, {0, 1}, 1}],
SFAD[{{P/2 - q1, 0}, {0, 1}, 1}],
SFAD[{{P/2 - q1 + q2, 0}, {1, 1}, 1}]}, {q1, q2, q3}, {k1,
P}, {Hold[Pair][Momentum[P, D], Momentum[P, D]] -> 4,
Hold[Pair][Momentum[k1, D], Momentum[k1, D]] -> 0,
Hold[Pair][Momentum[k1, D], Momentum[P, D]] -> 2}, {}],
FCTopology[
"qtopo1339", {SFAD[{{q3, 0}, {1, 1}, 1}],
SFAD[{{q2, 0}, {1, 1}, 1}], SFAD[{{q1, 0}, {1, 1}, 1}],
SFAD[{{k1 - q1, 0}, {1, 1}, 1}],
SFAD[{{k1 - P + q3, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q2, 0}, {0, 1}, 1}],
SFAD[{{k1 - P/2 + q2 + q3, 0}, {1, 1}, 1}],
SFAD[{{P/2 - q1, 0}, {0, 1}, 1}],
SFAD[{{P/2 - q1 - q2, 0}, {1, 1}, 1}],
SFAD[{{P/2 - q1 - q2 - q3, 0}, {0, 1}, 1}]}, {q1, q2, q3}, {k1,
P}, {Hold[Pair][Momentum[P, D], Momentum[P, D]] -> 4,
Hold[Pair][Momentum[k1, D], Momentum[k1, D]] -> 0,
Hold[Pair][Momentum[k1, D], Momentum[P, D]] -> 2}, {}],
FCTopology[
"qtopo1340", {SFAD[{{q3, 0}, {1, 1}, 1}],
SFAD[{{q2, 0}, {1, 1}, 1}], SFAD[{{q1, 0}, {1, 1}, 1}],
SFAD[{{k1 - q1, 0}, {1, 1}, 1}],
SFAD[{{k1 - P + q3, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q2, 0}, {0, 1}, 1}],
SFAD[{{P - q1 - q3, 0}, {0, 1}, 1}],
SFAD[{{P/2 - q1, 0}, {0, 1}, 1}],
SFAD[{{P/2 - q1 - q2, 0}, {1, 1}, 1}],
SFAD[{{P/2 - q1 - q2 - q3, 0}, {0, 1}, 1}]}, {q1, q2, q3}, {k1,
P}, {Hold[Pair][Momentum[P, D], Momentum[P, D]] -> 4,
Hold[Pair][Momentum[k1, D], Momentum[k1, D]] -> 0,
Hold[Pair][Momentum[k1, D], Momentum[P, D]] -> 2}, {}],
FCTopology[
"qtopo1341", {SFAD[{{q3, 0}, {1, 1}, 1}],
SFAD[{{q2, 0}, {1, 1}, 1}], SFAD[{{q1, 0}, {1, 1}, 1}],
SFAD[{{k1 - q1, 0}, {1, 1}, 1}],
SFAD[{{k1 - P + q3, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q2, 0}, {0, 1}, 1}],
SFAD[{{k1 - P/2 + q2, 0}, {0, 1}, 1}],
SFAD[{{k1 - P/2 + q2 + q3, 0}, {1, 1}, 1}],
SFAD[{{P/2 - q1, 0}, {0, 1}, 1}],
SFAD[{{P/2 - q1 - q2, 0}, {1, 1}, 1}]}, {q1, q2, q3}, {k1,
P}, {Hold[Pair][Momentum[P, D], Momentum[P, D]] -> 4,
Hold[Pair][Momentum[k1, D], Momentum[k1, D]] -> 0,
Hold[Pair][Momentum[k1, D], Momentum[P, D]] -> 2}, {}],
FCTopology[
"qtopo1342", {SFAD[{{q3, 0}, {1, 1}, 1}],
SFAD[{{q2, 0}, {1, 1}, 1}], SFAD[{{q1, 0}, {1, 1}, 1}],
SFAD[{{k1 - q1, 0}, {1, 1}, 1}],
SFAD[{{k1 - P + q3, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q2, 0}, {0, 1}, 1}],
SFAD[{{P/2 - q1, 0}, {0, 1}, 1}],
SFAD[{{P/2 - q1 - q2, 0}, {1, 1}, 1}],
SFAD[{{k1 - P/2 - q1 - q2, 0}, {1, 1}, 1}],
SFAD[{{P/2 - q1 - q2 - q3, 0}, {0, 1}, 1}]}, {q1, q2, q3}, {k1,
P}, {Hold[Pair][Momentum[P, D], Momentum[P, D]] -> 4,
Hold[Pair][Momentum[k1, D], Momentum[k1, D]] -> 0,
Hold[Pair][Momentum[k1, D], Momentum[P, D]] -> 2}, {}],
FCTopology[
"qtopo1345", {SFAD[{{q3, 0}, {0, 1}, 1}],
SFAD[{{q2, 0}, {0, 1}, 1}], SFAD[{{q1, 0}, {1, 1}, 1}],
SFAD[{{k1 - q2, 0}, {0, 1}, 1}],
SFAD[{{k1 - q2 + q3, 0}, {0, 1}, 1}],
SFAD[{{k1 - P + q1, 0}, {1, 1}, 1}],
SFAD[{{P/2 - q1, 0}, {0, 1}, 1}],
SFAD[{{k1 - P/2 - q2 + q3, 0}, {1, 1}, 1}],
SFAD[{{P/2 - q1 - q2, 0}, {0, 1}, 1}],
SFAD[{{P/2 - q1 - q2 + q3, 0}, {0, 1}, 1}]}, {q1, q2, q3}, {k1,
P}, {Hold[Pair][Momentum[P, D], Momentum[P, D]] -> 4,
Hold[Pair][Momentum[k1, D], Momentum[k1, D]] -> 0,
Hold[Pair][Momentum[k1, D], Momentum[P, D]] -> 2}, {}],
FCTopology[
"qtopo1356", {SFAD[{{q3, 0}, {0, 1}, 1}],
SFAD[{{q2, 0}, {0, 1}, 1}], SFAD[{{q1, 0}, {1, 1}, 1}],
SFAD[{{k1 - q2, 0}, {0, 1}, 1}],
SFAD[{{k1 - P + q1, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q3, 0}, {1, 1}, 1}],
SFAD[{{k1 - P/2 + q1 + q3, 0}, {0, 1}, 1}],
SFAD[{{P/2 - q1, 0}, {0, 1}, 1}],
SFAD[{{P/2 - q1 - q2, 0}, {0, 1}, 1}],
SFAD[{{P/2 - q1 - q2 - q3, 0}, {0, 1}, 1}]}, {q1, q2, q3}, {k1,
P}, {Hold[Pair][Momentum[P, D], Momentum[P, D]] -> 4,
Hold[Pair][Momentum[k1, D], Momentum[k1, D]] -> 0,
Hold[Pair][Momentum[k1, D], Momentum[P, D]] -> 2}, {}],
FCTopology[
"qtopo1367", {SFAD[{{q3, 0}, {0, 1}, 1}],
SFAD[{{q2, 0}, {0, 1}, 1}], SFAD[{{q1, 0}, {1, 1}, 1}],
SFAD[{{k1 - P + q1, 0}, {1, 1}, 1}],
SFAD[{{P - q1 + q2 + q3, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q3, 0}, {1, 1}, 1}],
SFAD[{{k1 - P + q1 - q2, 0}, {1, 1}, 1}],
SFAD[{{P/2 - q1, 0}, {0, 1}, 1}],
SFAD[{{P/2 - q1 + q2, 0}, {0, 1}, 1}],
SFAD[{{k1 - P + q1 - q2 - q3, 0}, {1, 1}, 1}]}, {q1, q2,
q3}, {k1, P}, {Hold[Pair][Momentum[P, D], Momentum[P, D]] -> 4,
Hold[Pair][Momentum[k1, D], Momentum[k1, D]] -> 0,
Hold[Pair][Momentum[k1, D], Momentum[P, D]] -> 2}, {}],
FCTopology[
"qtopo1378", {SFAD[{{q3, 0}, {0, 1}, 1}],
SFAD[{{q2, 0}, {0, 1}, 1}], SFAD[{{q1, 0}, {0, 1}, 1}],
SFAD[{{q1 + q2, 0}, {0, 1}, 1}],
SFAD[{{q1 + q2 - q3, 0}, {0, 1}, 1}],
SFAD[{{P/2 + q2, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q1, 0}, {1, 1}, 1}],
SFAD[{{P/2 - q2 + q3, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q2 - q3, 0}, {1, 1}, 1}],
SFAD[{{k1 - P/2 + q2 - q3, 0}, {1, 1}, 1}]}, {q1, q2, q3}, {k1,
P}, {Hold[Pair][Momentum[P, D], Momentum[P, D]] -> 4,
Hold[Pair][Momentum[k1, D], Momentum[k1, D]] -> 0,
Hold[Pair][Momentum[k1, D], Momentum[P, D]] -> 2}, {}],
FCTopology[
"qtopo1379", {SFAD[{{q3, 0}, {0, 1}, 1}],
SFAD[{{q2, 0}, {0, 1}, 1}], SFAD[{{q1, 0}, {1, 1}, 1}],
SFAD[{{k1 - P + q1, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q2, 0}, {1, 1}, 1}],
SFAD[{{P/2 - q1, 0}, {0, 1}, 1}],
SFAD[{{P/2 + q2 - q3, 0}, {1, 1}, 1}],
SFAD[{{k1 - P/2 - q2 + q3, 0}, {1, 1}, 1}],
SFAD[{{P/2 - q1 - q2, 0}, {0, 1}, 1}],
SFAD[{{P/2 - q1 - q2 + q3, 0}, {0, 1}, 1}]}, {q1, q2, q3}, {k1,
P}, {Hold[Pair][Momentum[P, D], Momentum[P, D]] -> 4,
Hold[Pair][Momentum[k1, D], Momentum[k1, D]] -> 0,
Hold[Pair][Momentum[k1, D], Momentum[P, D]] -> 2}, {}],
FCTopology[
"qtopo138", {SFAD[{{q3, 0}, {1, 1}, 1}],
SFAD[{{q2, 0}, {0, 1}, 1}], SFAD[{{q1, 0}, {0, 1}, 1}],
SFAD[{{q1 + q2, 0}, {0, 1}, 1}],
SFAD[{{k1 - P + q3, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q2, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q1, 0}, {1, 1}, 1}],
SFAD[{{k1 - P/2 + q1, 0}, {1, 1}, 1}],
SFAD[{{k1 - P/2 + q1 + q3, 0}, {0, 1}, 1}],
SFAD[{{k1 - P/2 - q2, 0}, {1, 1}, 1}]}, {q1, q2, q3}, {k1,
P}, {Hold[Pair][Momentum[P, D], Momentum[P, D]] -> 4,
Hold[Pair][Momentum[k1, D], Momentum[k1, D]] -> 0,
Hold[Pair][Momentum[k1, D], Momentum[P, D]] -> 2}, {}],
FCTopology[
"qtopo1390", {SFAD[{{q3, 0}, {0, 1}, 1}],
SFAD[{{q2, 0}, {0, 1}, 1}], SFAD[{{q1, 0}, {1, 1}, 1}],
SFAD[{{k1 - P + q1, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q2, 0}, {1, 1}, 1}],
SFAD[{{P/2 - q1, 0}, {0, 1}, 1}],
SFAD[{{k1 - P/2 - q2, 0}, {1, 1}, 1}],
SFAD[{{k1 - P/2 - q2 + q3, 0}, {1, 1}, 1}],
SFAD[{{P/2 - q1 - q2, 0}, {0, 1}, 1}],
SFAD[{{P/2 - q1 - q2 + q3, 0}, {0, 1}, 1}]}, {q1, q2, q3}, {k1,
P}, {Hold[Pair][Momentum[P, D], Momentum[P, D]] -> 4,
Hold[Pair][Momentum[k1, D], Momentum[k1, D]] -> 0,
Hold[Pair][Momentum[k1, D], Momentum[P, D]] -> 2}, {}],
FCTopology[
"qtopo1401", {SFAD[{{q3, 0}, {0, 1}, 1}],
SFAD[{{q2, 0}, {0, 1}, 1}], SFAD[{{q1, 0}, {1, 1}, 1}],
SFAD[{{k1 - q1, 0}, {1, 1}, 1}],
SFAD[{{k1 - P + q3, 0}, {0, 1}, 1}],
SFAD[{{P/2 + q2, 0}, {1, 1}, 1}],
SFAD[{{P/2 - q1, 0}, {0, 1}, 1}],
SFAD[{{P/2 - q1 - q2, 0}, {0, 1}, 1}],
SFAD[{{k1 - P/2 - q1 - q2, 0}, {0, 1}, 1}],
SFAD[{{P/2 - q1 - q2 - q3, 0}, {0, 1}, 1}]}, {q1, q2, q3}, {k1,
P}, {Hold[Pair][Momentum[P, D], Momentum[P, D]] -> 4,
Hold[Pair][Momentum[k1, D], Momentum[k1, D]] -> 0,
Hold[Pair][Momentum[k1, D], Momentum[P, D]] -> 2}, {}],
FCTopology[
"qtopo1432", {SFAD[{{q3, 0}, {0, 1}, 1}],
SFAD[{{q2, 0}, {0, 1}, 1}], SFAD[{{q1, 0}, {1, 1}, 1}],
SFAD[{{k1 - q1, 0}, {1, 1}, 1}],
SFAD[{{k1 - q1 + q2, 0}, {1, 1}, 1}],
SFAD[{{k1 - q1 + q2 + q3, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q3, 0}, {1, 1}, 1}],
SFAD[{{k1 - P/2 + q3, 0}, {1, 1}, 1}],
SFAD[{{P/2 - q1, 0}, {0, 1}, 1}],
SFAD[{{P/2 - q1 + q2, 0}, {0, 1}, 1}]}, {q1, q2, q3}, {k1,
P}, {Hold[Pair][Momentum[P, D], Momentum[P, D]] -> 4,
Hold[Pair][Momentum[k1, D], Momentum[k1, D]] -> 0,
Hold[Pair][Momentum[k1, D], Momentum[P, D]] -> 2}, {}],
FCTopology[
"qtopo1481", {SFAD[{{q3, 0}, {0, 1}, 1}],
SFAD[{{q2, 0}, {1, 1}, 1}], SFAD[{{q1, 0}, {0, 1}, 1}],
SFAD[{{q1 + q2, 0}, {1, 1}, 1}],
SFAD[{{P + q1 + q2 + q3, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q2, 0}, {0, 1}, 1}],
SFAD[{{P/2 + q1, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q1 + q3, 0}, {1, 1}, 1}],
SFAD[{{k1 - P - q1 - q2, 0}, {1, 1}, 1}],
SFAD[{{k1 - P - q1 - q2 - q3, 0}, {1, 1}, 1}]}, {q1, q2,
q3}, {k1, P}, {Hold[Pair][Momentum[P, D], Momentum[P, D]] -> 4,
Hold[Pair][Momentum[k1, D], Momentum[k1, D]] -> 0,
Hold[Pair][Momentum[k1, D], Momentum[P, D]] -> 2}, {}],
FCTopology[
"qtopo1491", {SFAD[{{q3, 0}, {0, 1}, 1}],
SFAD[{{q2, 0}, {1, 1}, 1}], SFAD[{{q1, 0}, {0, 1}, 1}],
SFAD[{{q1 + q2, 0}, {1, 1}, 1}],
SFAD[{{q1 + q2 + q3, 0}, {1, 1}, 1}],
SFAD[{{k1 - P + q2, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q3, 0}, {1, 1}, 1}],
SFAD[{{P/2 + q1, 0}, {1, 1}, 1}],
SFAD[{{k1 - P/2 - q3, 0}, {1, 1}, 1}],
SFAD[{{P/2 - q2 - q3, 0}, {0, 1}, 1}]}, {q1, q2, q3}, {k1,
P}, {Hold[Pair][Momentum[P, D], Momentum[P, D]] -> 4,
Hold[Pair][Momentum[k1, D], Momentum[k1, D]] -> 0,
Hold[Pair][Momentum[k1, D], Momentum[P, D]] -> 2}, {}]};

subtopos = FCLoopFindSubtopologies[topos1, FCParallelize -> True];

FCLoopFindTopologyMappings[topos1, PreferredTopologies -> subtopos,
Momentum -> All, FCParallelize -> False, FCVerbose -> 1];

FCLoopFindTopologyMappings: Entering.

FCLoopFindTopologyMappings: Calling FCLoopFindIntegralMappings.

FCLoopFindTopologyMappings: FCLoopFindIntegralMappings done, timing: 456.9

FCLoopFindTopologyMappings: Calling findMappings.

[ygtw1]

The version of Mathematica is 13.1.0.0.

[ygtw1]

I have found the problem comes from the bug of SelectNotFree.
SelectNotFree[expr, x] will quit when Lengh[x]>15000.

[ygtw1]

I have two ways to solve this problem.

One is replacing the code target= SelectNotFree[input, preferred]; in line 226 of FCLoopFindTopologyMappings.m with target= SelectNotFree[input, Intersection[preferred,Flatten[{#[[1]],Values[#[[-1]]]}&/@Flatten[input]]]];

The other is replacing it with target =Block[{myFreeQ}, myFreeQ[expr_, xlist_] := AllTrue[xlist, FreeQ[expr, #] &]; Select[input, ! myFreeQ[#, preferred] &]];

[vsht]

I've also encountered such issues in the past, mostly related to Select(Not)Free(2) and FreeQ2. To have an arbitrary number of arguments after the first slot these functions use recursion. This mostly works well in the daily life but was never intended for lists with thousands or tens of thousands entries as they appear when doing topology mappings.

In many cases one can rewrite the relevant code pieces using Alternatives, since the list of objects to check for consists of atomic symbols. If the objects are more complicated, one can easily run into situations where Alternatives won't work correctly, though.

In general, it's a lot of case by case considerations, though. Sometimes a solution that looks good will kill the performance, sometimes it will work fine but fail when the number of list elements is doubled or tripled.

Apart from that, Mathematica sub kernels sometimes happen to die unexpected when the load gets too high. This is something one can't do much about as a user.

[vsht]

You can try it with the new FreeQ2. In general, I don't recommend using FCLoopFindSubtopologies, as it rarely gives you more than a few more mappings but takes ages to run.

[vsht]

I think that Commit 2fbce53 should have taken care of it, because there's no more recursion when using FreeQ2 and other routines relying on it. Does the issue still occur?

[ygtw1]

In[1]:= << FeynCalc`

During evaluation of In[1]:= FeynCalc 10.1.0 (dev version). For help, use the online documentation, visit the forum and have a look at the supplied examples. The PDF-version of the manual can be downloaded here.

During evaluation of In[1]:= If you use FeynCalc in your research, please evaluate FeynCalcHowToCite[] to learn how to cite this software.

During evaluation of In[1]:= Please keep in mind that the proper academic attribution of our work is crucial to ensure the future development of this package!

In[2]:= input =
FCTopology[
"x" <> ToString[#], {}, {}, {}, {}, {FCGV[SubtopologyOf] ->
"x"}] & /@ Range[1000] // RandomSample;
preferred = "x" <> ToString[#] & /@ Range[10^6] // RandomSample;
Select[input, !
FreeQ2[preferred,
Flatten[{#[[1]], #[[-1]] /. Rule -> List}]] &] // Timing;
SelectNotFree[input, preferred] // Timing;
{%%[[1]], %[[1]], %%[[2]] === %[[2]]}

Out[6]= {13.4572, 54.3971, True}

I believe there may still be room for optimizing the program to improve its speed, though the impact is likely to be minimal. This issue can be considered resolved.

[vsht]

I think with Intersection one could make it even faster

input = FCTopology[
      "x" <> ToString[#], {}, {}, {}, {}, {FCGV[SubtopologyOf] -> 
        "x"}] & /@ Range[1000] // RandomSample;
preferred = "x" <> ToString[#] & /@ Range[10^6] // RandomSample;
AbsoluteTiming[
 res1 = SelectNotFree[input, Intersection[First /@ input, preferred]];]
(*{0.895677,Null}*)

[ygtw1]

In[1]:= << FeynCalc`

During evaluation of In[1]:= FeynCalc 10.1.0 (dev version). For help, use the online documentation, visit the forum and have a look at the supplied examples. The PDF-version of the manual can be downloaded here.

During evaluation of In[1]:= If you use FeynCalc in your research, please evaluate FeynCalcHowToCite[] to learn how to cite this software.

During evaluation of In[1]:= Please keep in mind that the proper academic attribution of our work is crucial to ensure the future development of this package!

In[2]:= input =FCTopology["x" <> ToString[#], {}, {}, {}, {}, {FCGV[SubtopologyOf] ->"x"}] & /@ Range[-10^3,10^4] // RandomSample;
preferred = "x" <> ToString[#] & /@ Range[10^6] // RandomSample;
way1=With[{temp=Intersection[First /@ input, preferred]},Select[input, MemberQ[temp,#[[1]]]&]]// AbsoluteTiming;
%[[1]]
way2=SelectNotFree[input, Intersection[First /@ input, preferred]]// AbsoluteTiming;
%[[1]]
way3=Select[input, !FreeQ2[preferred,Flatten[{#[[1]], #[[-1]] /. Rule -> List}]] &] // AbsoluteTiming;
%[[1]]
way4=SelectNotFree[input, preferred] // AbsoluteTiming;
%[[1]]
way1[[2]]===way2[[2]]===way3[[2]]===way4[[2]]

Out[5]= 0.936864

Out[7]= 2.48037

Out[9]= 153.987

Out[11]= 760.005

Out[12]= True

[vsht]

Looks good. Actually, in this case Intersection[First /@ input, preferred] is what takes most of the time. You can speed it up even further by switching to associations and using key look ups.

input = FCTopology[
      "x" <> ToString[#], {}, {}, {}, {}, {FCGV[SubtopologyOf] -> 
        "x"}] & /@ Range[-10^3, 10^5] // RandomSample;
preferred = "x" <> ToString[#] & /@ Range[10^6] // RandomSample;

AbsoluteTiming[
 assoc = AssociationThread[preferred -> True];
 res1 = Select[input, KeyExistsQ[assoc, #[[1]]] &];]
AbsoluteTiming[
 res2 = With[{temp = Intersection[First /@ input, preferred]}, 
    Select[input, MemberQ[temp, #[[1]]] &]];]

(*Out[74]= {1.04678,Null}

Out[75]= {83.1486,Null}*)

Sort[res1] === Sort[res2]
(*True*)

It probably would make sense to write a dedicated routine for selecting topologies and integrals. The topology business is something I would still like to do in Mathematica for a while. For most other manipulations I'm planning to switch to FORM
via LoopScalla which will supersede FeynCalc.