Name: Philipp
Date: 06/07/17-01:04:11 PM Z
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