Name: Vladyslav Shtabovenko (email_not_shown)
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
>



This archive was generated by hypermail 2b29 : 11/23/17-07:20:01 PM Z CET