Name: Vladyslav Shtabovenko Date: 06/06/17-06:30:44 PM Z
Hi,
could you provide a 2-loop example where it takes too much time?
For my own purposes I used to employ FVSeries
exp = -(CA^3*
FAD[{p1, I*mE}, {p2, I*mE}, {p1 - q, I*mE}, {p1 - q,
I*mE}]*(FVD[p1, Lor1] + FVD[p1 - q,
Lor1])*(FVD[p1, Lor2] +
FVD[p1 - q, Lor2])*(2*omega[21] -
5*psi[19])*SD[a, b]*
SMP[“g_s”]^6*(SPD[p2, p2] + SPD[p1 - q, p1 -
q]))
ClearAll[FVSeries];
Options[FVSeries] = {List -> False, Dimension -> D};
FVSeries[expr_, {q_, q0_, n_}, OptionsPattern[]]
:=
Block[{fvs, ex, res, dim}, dim = OptionValue[Dimension];
ex = FCI[expr];
res = Join[{EpsEvaluate[(ex /. q -> q0)]},
Table[fvs =
Table[FCI[
Pair[Momentum[q, dim],
LorentzIndex[Unique[], dim]]], {j, 1,
i}]; ((1/i!) FourDivergence[ex, Sequence @@ fvs,
Contract -> False]) // ReplaceAll[#, q -> q0]
& //
Contract[# Times @@ fvs] &, {i, 1, n}]];
If[! OptionValue[List], res = Total@res];
res]
The following requires around 12 second on my machine with FeynCalc 9.3
exp = FVSeries[exp, {q, 0, 6}]; // AbsoluteTiming
I guess one could also make it a bit faster, but I’m not sure how long
time does it require in your 2-loop calculation.
Cheers,
Vladyslav
Am 06.06.2017 um 16:30 schrieb Philipp:
> Hi,
>
> I wanted to expand the following expression in terms of the external
momentum q or (q/mE)
>
*> -(CA^3*FAD[{p1, I*mE}, {p2, I*mE}, {p1 - q, I*mE}, {p1 - q,
I*mE}]*(FVD[p1, Lor1] + FVD[p1 - q,
Lor1])*(FVD[p1, Lor2] + FVD[p1 - q,
Lor2])*(2*omega[21] - 5*psi[19])*SD[a,
b]*SMP[“g_s”]^6*(SPD[p2, p2] + SPD[p1 - q, p1