Name: Vladyslav Shtabovenko
Date: 06/10/17-09:33:09 PM Z
Hi,
sorry for the late reply. What about the following:
exp2 = Collect2[FCI[exp], Momentum, IsolateNames -> KK]
// Simplify
res = FVSeries[exp2, {q, 0, 4}]; // AbsoluteTiming
where exp is your original expression.
This requires 3.7 seconds on my machine with FeynCalc 9.2
Doing the expansion directly via
res = FVSeries[exp, {q, 0, 4}]; // AbsoluteTiming
takes 19.4 seconds, i.e. almost 5 times longer.
If 3.7 seconds are still too slow for you, I guess you would need
to go with FORM.
Cheers,
Vladyslav
Am 07.06.2017 um 13:04 schrieb Philipp:
> Hello,
>
> I have actually done a similar thing just that I have computed the
FourDivergences one after another to not derive already existing
expressions again. But that is probably as fast as it gets right?
Because I have expressions much longer than the said one like
>
*> -(CA*(2*FVD[p1, Lor1] - FVD[q,
Lor1])*(2*FVD[p1, Lor2] - FVD[q, Lor2])*SD[a,
b]*SMP[“g_s”]^6*(FAD[{p1, I*mE}, p2, {p1 - q,
I*mE}, {p1 - q, I*mE}]*((-(D*(psi[4] +
2*psi[13] - 2*psi[15])) - psi[26] -
2*psi[34] + CA^2*(omega[26] + D*(omega[4] -
2*(psi[4] + psi[13] - psi[15])) -
2*(psi[26] + psi[34] - psi[35])) -
2*CA*CF*(D*(psi[4] + 2*psi[13] -
2*psi[15]) + psi[26] + 2*psi[34] -
2*psi[35]) + 2*psi[35])*SPD[p1, p1] +
2*(D*(psi[4] + 2*psi[13] - 2*psi[15]) +
psi[26] + 2*psi[34] - CA^2*(omega[26] +
D*(omega[4] - 2*(psi[4] + psi[13] -
psi[15])) - 2*(psi[26] + psi[34] -
psi[35])) + 2*CA*CF*(D*(psi[4] + 2*psi[13]
- 2*psi[15]) + psi[26] + 2*psi[34] -
2*psi[35]) - 2*psi[35])*SPD[p1, q] +
CA^2*D*omega[5]*SPD[p2, p2] +
CA^2*omega[29]*SPD[p2, p2] -
D*psi[5]*SPD[p2, p2] -
2*CA^2*D*psi[5]*SPD[p2, p2] -
2*CA*CF*D*psi[5]*SPD[p2, p2] -
2*D*psi[16]*SPD[p2, p2] -
2*CA^2*D*psi[16]*SPD[p2, p2] -
4*CA*CF*D*psi[16]*SPD[p2, p2] +
2*D*psi[18]*SPD[p2, p2] +
2*CA^2*D*psi[18]*SPD[p2, p2] +
4*CA*CF*D*psi[18*
> ]*SPD[p2, p2] - psi[29]*SPD[p2, p2] -
2*CA^2*psi[29]*SPD[p2, p2] -
2*CA*CF*psi[29]*SPD[p2, p2] -
2*psi[42]*SPD[p2, p2] -
2*CA^2*psi[42]*SPD[p2, p2] -
4*CA*CF*psi[42]*SPD[p2, p2] +
2*psi[43]*SPD[p2, p2] +
2*CA^2*psi[43]*SPD[p2, p2] +
4*CA*CF*psi[43]*SPD[p2, p2] +
2*psi[44]*SPD[p2, p2] +
2*CA^2*psi[44]*SPD[p2, p2] +
4*CA*CF*psi[44]*SPD[p2, p2] +
CA^2*D*omega[4]*SPD[q, q] +
CA^2*omega[26]*SPD[q, q] -
D*psi[4]*SPD[q, q] -
2*CA^2*D*psi[4]*SPD[q, q] -
2*CA*CF*D*psi[4]*SPD[q, q] -
2*D*psi[13]*SPD[q, q] -
2*CA^2*D*psi[13]*SPD[q, q] -
4*CA*CF*D*psi[13]*SPD[q, q] +
2*D*psi[15]*SPD[q, q] +
2*CA^2*D*psi[15]*SPD[q, q] +
4*CA*CF*D*psi[15]*SPD[q, q] -
psi[26]*SPD[q, q] -
2*CA^2*psi[26]*SPD[q, q] -
2*CA*CF*psi[26]*SPD[q, q] -
2*psi[34]*SPD[q, q] -
2*CA^2*psi[34]*SPD[q, q] -
4*CA*CF*psi[34]*SPD[q, q] +
2*psi[35]*SPD[q, q] +
2*CA^2*psi[35]*SPD[q, q] +
4*CA*CF*psi[35]*SPD[q, q]) + (-1 +
GaugeXi)*FAD[{p1, I*mE}, p2, p2, {p1 - q, I*mE}, {p1 - q,
I*mE}]*((-psi[26] - 2*psi[34] +
CA^2*(omega[26] -
> 2*(psi[26] + psi[34] - psi[35])) -
2*CA*CF*(psi[26] + 2*psi[34] - 2*psi[35]) +
2*psi[35])*SPD[p1, p2]^2 + (-psi[4] -
2*psi[13] + CA^2*(omega[4] - 2*(psi[4] +
psi[13] - psi[15])) - 2*CA*CF*(psi[4] +
2*psi[13] - 2*psi[15]) +
2*psi[15])*SPD[p1, p1]*SPD[p2, p2] -
2*CA^2*omega[4]*SPD[p1, q]*SPD[p2, p2] +
2*psi[4]*SPD[p1, q]*SPD[p2, p2] +
4*CA^2*psi[4]*SPD[p1, q]*SPD[p2, p2] +
4*CA*CF*psi[4]*SPD[p1, q]*SPD[p2, p2] +
4*psi[13]*SPD[p1, q]*SPD[p2, p2] +
4*CA^2*psi[13]*SPD[p1, q]*SPD[p2, p2] +
8*CA*CF*psi[13]*SPD[p1, q]*SPD[p2, p2] -
4*psi[15]*SPD[p1, q]*SPD[p2, p2] -
4*CA^2*psi[15]*SPD[p1, q]*SPD[p2, p2] -
8*CA*CF*psi[15]*SPD[p1, q]*SPD[p2, p2] +
CA^2*omega[5]*SPD[p2, p2]^2 +
CA^2*omega[29]*SPD[p2, p2]^2 -
psi[5]*SPD[p2, p2]^2 -
2*CA^2*psi[5]*SPD[p2, p2]^2 -
2*CA*CF*psi[5]*SPD[p2, p2]^2 -
2*psi[16]*SPD[p2, p2]^2 -
2*CA^2*psi[16]*SPD[p2, p2]^2 -
4*CA*CF*psi[16]*SPD[p2, p2]^2 +
2*psi[18]*SPD[p2, p2]^2 +
2*CA^2*psi[18]*SPD[p2, p2]^2 +
4*CA*CF*psi[18]*SPD[p2, p2]^2 -
psi[29]*SPD[p2, p2]^
*> 2 - 2*CA^2*psi[29]*SPD[p2, p2]^2 -
2*CA*CF*psi[29]*SPD[p2, p2]^2 -
2*psi[42]*SPD[p2, p2]^2 -
2*CA^2*psi[42]*SPD[p2, p2]^2 -
4*CA*CF*psi[42]*SPD[p2, p2]^2 +
2*psi[43]*SPD[p2, p2]^2 +
2*CA^2*psi[43]*SPD[p2, p2]^2 +
4*CA*CF*psi[43]*SPD[p2, p2]^2 +
2*psi[44]*SPD[p2, p2]^2 +
2*CA^2*psi[44]*SPD[p2, p2]^2 +
4*CA*CF*psi[44]*SPD[p2, p2]^2 + 2*(psi[26]
- 2*psi[34] - CA^2*(omega[26] - 2*(psi[26] +
psi[34] - psi[35])) + 2*CA*CF*(psi[26] +
2*psi[34] - 2*psi[35]) -
2*psi[35])*SPD[p1, p2]*SPD[p2, q] +
CA^2*omega[26]*SPD[p2, q]^2 -
psi[26]*SPD[p2, q]^2 -
2*CA^2*psi[26]*SPD[p2, q]^2 -
2*CA*CF*psi[26]*SPD[p2, q]^2 -
2*psi[34]*SPD[p2, q]^2 -
2*CA^2*psi[34]*SPD[p2, q]^2 -
4*CA*CF*psi[34]*SPD[p2, q]^2 +
2*psi[35]*SPD[p2, q]^2 +
2*CA^2*psi[35]*SPD[p2, q]^2 +
4*CA*CF*psi[35]*SPD[p2, q]^2 +
CA^2*omega[4]*SPD[p2, p2]*SPD[q, q] -
psi[4]*SPD[p2, p2]*SPD[q, q] -
2*CA^2*psi[4]*SPD[p2, p2]*SPD[q, q] -
2*CA*CF*psi[4]*SPD[p2, p2]*SPD[q, q] -
2*psi[13]*SPD[p2, p2]*SPD[q, q] -
2*CA^2*psi[13]*SPD[p2, p2]*SPD[*
> q, q] - 4*CA*CF*psi[13]*SPD[p2,
p2]*SPD[q, q] + 2*psi[15]*SPD[p2,
p2]*SPD[q, q] + 2*CA^2*psi[15]*SPD[p2,
p2]*SPD[q, q] + 4*CA*CF*psi[15]*SPD[p2,
p2]*SPD[q, q])))/2
>
> or even longer. Those will take up significant time but I haven’t
yet figured out how to speed that up or simplify it somehow.
>
> Cheers,
> Philipp
>