Skip to content

A bug when FCParallelizing FCLoopCreatePartialFractioningRules with assigning numbers to symbols after ScalarProducts=symbols #348

@ygtw1

Description

@ygtw1

"test1:"
Quit[]
LaunchKernels[8];
$LoadAddOns={"FeynArts"};
<<FeynCalc`;
$FAVerbose=0;
$ParallelizeFeynCalc=True;
FCClearScalarProducts[];
SPD[P]=4 x;
SPD[k1]=0;
SPD[P,k1]=2 x;
x=1;
glis={GLI["fctopology",{1,1,1,1,1,2,1,1,1}]};
topos={FCTopology["fctopology",
{SFAD[{{q3,0},{1,1},1}],SFAD[{{q2,0},{0,1},1}],SFAD[{{q1,0},{1,1},1}],SFAD[{{q1+q2,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[{{P/2-q2,0},{1,1},1}],SFAD[{{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},{}]};
FCLoopCreatePartialFractioningRules[glis,topos,FCParallelize->True,"KeepApartHead"->True]//FCE
FCLoopCreatePartialFractioningRules[glis,topos,FCParallelize->False,"KeepApartHead"->True]//FCE
%%===%

During evaluation of In[1]:= FeynCalc 10.2.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!

During evaluation of In[1]:= FeynArts 3.12 (24 May 2024) patched for use with FeynCalc, for documentation see the manual or visit www.feynarts.de.

During evaluation of In[1]:= If you use FeynArts in your research, please cite

During evaluation of In[1]:= [Bullet] T. Hahn, Comput. Phys. Commun., 140, 418-431, 2001, arXiv:hep-ph/0012260

During evaluation of In[1]:= Power::infy: Infinite expression 1/0 encountered.

During evaluation of In[1]:= Power::infy: Infinite expression 1/0 encountered.

During evaluation of In[1]:= Power::infy: Infinite expression 1/0 encountered.

During evaluation of In[1]:= General::stop: Further output of Power::infy will be suppressed during this calculation.

Out[13]= {{GLI["fctopology", {1, 1, 1, 1, 1, 2, 1, 1, 1}] ->
FCGV["ApartHead"][ComplexInfinity,
GLI["fctopologyPFR2", {1, 1, 1, 1, 2, 1, 1, 1}]] +
FCGV["ApartHead"][ComplexInfinity,
GLI["fctopologyPFR5", {1, 1, 1, 1, 2, 1, 1, 1}]] +
FCGV["ApartHead"][ComplexInfinity,
GLI["fctopologyPFR8", {1, 1, 1, 1, 1, 2, 1, 1}]]}, {FCTopology[
"fctopologyPFR2", {SFAD[{{q3, 0}, {1, 1}, 1}],
SFAD[{{q1, 0}, {1, 1}, 1}], SFAD[{{q1 + q2, 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[{{P/2 - q2, 0}, {1, 1}, 1}],
SFAD[{{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[
"fctopologyPFR5", {SFAD[{{q3, 0}, {1, 1}, 1}],
SFAD[{{q2, 0}, {0, 1}, 1}], SFAD[{{q1, 0}, {1, 1}, 1}],
SFAD[{{q1 + q2, 0}, {1, 1}, 1}], SFAD[{{P/2 + q1, 0}, {0, 1}, 1}],
SFAD[{{k1 - P/2 + q2, 0}, {1, 1}, 1}],
SFAD[{{P/2 - q2, 0}, {1, 1}, 1}],
SFAD[{{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[
"fctopologyPFR8", {SFAD[{{q3, 0}, {1, 1}, 1}],
SFAD[{{q2, 0}, {0, 1}, 1}], SFAD[{{q1, 0}, {1, 1}, 1}],
SFAD[{{q1 + q2, 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[{{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}, {}]}}

Out[14]= {{GLI["fctopology", {1, 1, 1, 1, 1, 2, 1, 1, 1}] ->
FCGV["ApartHead"][-1,
GLI["fctopologyPFR5", {1, 1, 1, 1, 2, 1, 2, 1}]] +
FCGV["ApartHead"][2,
GLI["fctopologyPFR2", {1, 1, 1, 1, 2, 1, 2, 1}]]}, {FCTopology[
"fctopologyPFR5", {SFAD[{{q3, 0}, {1, 1}, 1}],
SFAD[{{q2, 0}, {0, 1}, 1}], SFAD[{{q1, 0}, {1, 1}, 1}],
SFAD[{{q1 + q2, 0}, {1, 1}, 1}], SFAD[{{P/2 + q1, 0}, {0, 1}, 1}],
SFAD[{{k1 - P/2 + q2, 0}, {1, 1}, 1}],
SFAD[{{P/2 - q2, 0}, {1, 1}, 1}],
SFAD[{{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[
"fctopologyPFR2", {SFAD[{{q3, 0}, {1, 1}, 1}],
SFAD[{{q1, 0}, {1, 1}, 1}], SFAD[{{q1 + q2, 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[{{P/2 - q2, 0}, {1, 1}, 1}],
SFAD[{{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}, {}]}}

Out[15]= False


"test2"
Quit[]
LaunchKernels[8];
$LoadAddOns={"FeynArts"};
<<FeynCalc`;
$FAVerbose=0;
$ParallelizeFeynCalc=True;
FCClearScalarProducts[];
x=1;
SPD[P]=4 x;
SPD[k1]=0;
SPD[P,k1]=2 x;
glis={GLI["fctopology",{1,1,1,1,1,2,1,1,1}]};
topos={FCTopology["fctopology",
{SFAD[{{q3,0},{1,1},1}],SFAD[{{q2,0},{0,1},1}],SFAD[{{q1,0},{1,1},1}],SFAD[{{q1+q2,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[{{P/2-q2,0},{1,1},1}],SFAD[{{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},{}]};
FCLoopCreatePartialFractioningRules[glis,topos,FCParallelize->True,"KeepApartHead"->True]//FCE
FCLoopCreatePartialFractioningRules[glis,topos,FCParallelize->False,"KeepApartHead"->True]//FCE
%%===%

During evaluation of In[1]:= FeynCalc 10.2.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!

During evaluation of In[1]:= FeynArts 3.12 (24 May 2024) patched for use with FeynCalc, for documentation see the manual or visit www.feynarts.de.

During evaluation of In[1]:= If you use FeynArts in your research, please cite

During evaluation of In[1]:= [Bullet] T. Hahn, Comput. Phys. Commun., 140, 418-431, 2001, arXiv:hep-ph/0012260

Out[13]= {{GLI["fctopology", {1, 1, 1, 1, 1, 2, 1, 1, 1}] ->
FCGV["ApartHead"][-1,
GLI["fctopologyPFR5", {1, 1, 1, 1, 2, 1, 2, 1}]] +
FCGV["ApartHead"][2,
GLI["fctopologyPFR2", {1, 1, 1, 1, 2, 1, 2, 1}]]}, {FCTopology[
"fctopologyPFR5", {SFAD[{{q3, 0}, {1, 1}, 1}],
SFAD[{{q2, 0}, {0, 1}, 1}], SFAD[{{q1, 0}, {1, 1}, 1}],
SFAD[{{q1 + q2, 0}, {1, 1}, 1}], SFAD[{{P/2 + q1, 0}, {0, 1}, 1}],
SFAD[{{k1 - P/2 + q2, 0}, {1, 1}, 1}],
SFAD[{{P/2 - q2, 0}, {1, 1}, 1}],
SFAD[{{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[
"fctopologyPFR2", {SFAD[{{q3, 0}, {1, 1}, 1}],
SFAD[{{q1, 0}, {1, 1}, 1}], SFAD[{{q1 + q2, 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[{{P/2 - q2, 0}, {1, 1}, 1}],
SFAD[{{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}, {}]}}

Out[14]= {{GLI["fctopology", {1, 1, 1, 1, 1, 2, 1, 1, 1}] ->
FCGV["ApartHead"][-1,
GLI["fctopologyPFR5", {1, 1, 1, 1, 2, 1, 2, 1}]] +
FCGV["ApartHead"][2,
GLI["fctopologyPFR2", {1, 1, 1, 1, 2, 1, 2, 1}]]}, {FCTopology[
"fctopologyPFR5", {SFAD[{{q3, 0}, {1, 1}, 1}],
SFAD[{{q2, 0}, {0, 1}, 1}], SFAD[{{q1, 0}, {1, 1}, 1}],
SFAD[{{q1 + q2, 0}, {1, 1}, 1}], SFAD[{{P/2 + q1, 0}, {0, 1}, 1}],
SFAD[{{k1 - P/2 + q2, 0}, {1, 1}, 1}],
SFAD[{{P/2 - q2, 0}, {1, 1}, 1}],
SFAD[{{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[
"fctopologyPFR2", {SFAD[{{q3, 0}, {1, 1}, 1}],
SFAD[{{q1, 0}, {1, 1}, 1}], SFAD[{{q1 + q2, 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[{{P/2 - q2, 0}, {1, 1}, 1}],
SFAD[{{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}, {}]}}

Out[15]= True


"test3"
Quit[]
LaunchKernels[8];
$LoadAddOns={"FeynArts"};
<<FeynCalc`;
$FAVerbose=0;
$ParallelizeFeynCalc=True;
FCClearScalarProducts[];
SPD[P]=4 x;
SPD[k1]=0;
SPD[P,k1]=2 x;
glis={GLI["fctopology",{1,1,1,1,1,2,1,1,1}]};
topos={FCTopology["fctopology",
{SFAD[{{q3,0},{1,1},1}],SFAD[{{q2,0},{0,1},1}],SFAD[{{q1,0},{1,1},1}],SFAD[{{q1+q2,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[{{P/2-q2,0},{1,1},1}],SFAD[{{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},{}]};
FCLoopCreatePartialFractioningRules[glis,topos,FCParallelize->True,"KeepApartHead"->True]//FCE
FCLoopCreatePartialFractioningRules[glis,topos,FCParallelize->False,"KeepApartHead"->True]//FCE
%%===%

During evaluation of In[1]:= FeynCalc 10.2.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!

During evaluation of In[1]:= FeynArts 3.12 (24 May 2024) patched for use with FeynCalc, for documentation see the manual or visit www.feynarts.de.

During evaluation of In[1]:= If you use FeynArts in your research, please cite

During evaluation of In[1]:= [Bullet] T. Hahn, Comput. Phys. Commun., 140, 418-431, 2001, arXiv:hep-ph/0012260

Out[12]= {{GLI["fctopology", {1, 1, 1, 1, 1, 2, 1, 1, 1}] ->
FCGV["ApartHead"][-(1/(-1 + x)),
GLI["fctopologyPFR2", {1, 1, 1, 1, 2, 1, 1, 1}]] +
FCGV["ApartHead"][1/(2 (-1 + x)),
GLI["fctopologyPFR5", {1, 1, 1, 1, 2, 1, 1, 1}]] +
FCGV["ApartHead"][1/(2 (-1 + x)),
GLI["fctopologyPFR8", {1, 1, 1, 1, 1, 2, 1, 1}]]}, {FCTopology[
"fctopologyPFR2", {SFAD[{{q3, 0}, {1, 1}, 1}],
SFAD[{{q1, 0}, {1, 1}, 1}], SFAD[{{q1 + q2, 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[{{P/2 - q2, 0}, {1, 1}, 1}],
SFAD[{{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[
"fctopologyPFR5", {SFAD[{{q3, 0}, {1, 1}, 1}],
SFAD[{{q2, 0}, {0, 1}, 1}], SFAD[{{q1, 0}, {1, 1}, 1}],
SFAD[{{q1 + q2, 0}, {1, 1}, 1}], SFAD[{{P/2 + q1, 0}, {0, 1}, 1}],
SFAD[{{k1 - P/2 + q2, 0}, {1, 1}, 1}],
SFAD[{{P/2 - q2, 0}, {1, 1}, 1}],
SFAD[{{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[
"fctopologyPFR8", {SFAD[{{q3, 0}, {1, 1}, 1}],
SFAD[{{q2, 0}, {0, 1}, 1}], SFAD[{{q1, 0}, {1, 1}, 1}],
SFAD[{{q1 + q2, 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[{{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}, {}]}}

Out[13]= {{GLI["fctopology", {1, 1, 1, 1, 1, 2, 1, 1, 1}] ->
FCGV["ApartHead"][-(1/(-1 + x)),
GLI["fctopologyPFR2", {1, 1, 1, 1, 2, 1, 1, 1}]] +
FCGV["ApartHead"][1/(2 (-1 + x)),
GLI["fctopologyPFR5", {1, 1, 1, 1, 2, 1, 1, 1}]] +
FCGV["ApartHead"][1/(2 (-1 + x)),
GLI["fctopologyPFR8", {1, 1, 1, 1, 1, 2, 1, 1}]]}, {FCTopology[
"fctopologyPFR2", {SFAD[{{q3, 0}, {1, 1}, 1}],
SFAD[{{q1, 0}, {1, 1}, 1}], SFAD[{{q1 + q2, 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[{{P/2 - q2, 0}, {1, 1}, 1}],
SFAD[{{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[
"fctopologyPFR5", {SFAD[{{q3, 0}, {1, 1}, 1}],
SFAD[{{q2, 0}, {0, 1}, 1}], SFAD[{{q1, 0}, {1, 1}, 1}],
SFAD[{{q1 + q2, 0}, {1, 1}, 1}], SFAD[{{P/2 + q1, 0}, {0, 1}, 1}],
SFAD[{{k1 - P/2 + q2, 0}, {1, 1}, 1}],
SFAD[{{P/2 - q2, 0}, {1, 1}, 1}],
SFAD[{{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[
"fctopologyPFR8", {SFAD[{{q3, 0}, {1, 1}, 1}],
SFAD[{{q2, 0}, {0, 1}, 1}], SFAD[{{q1, 0}, {1, 1}, 1}],
SFAD[{{q1 + q2, 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[{{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}, {}]}}

Out[14]= True

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions