Commit eade9294 authored by ulrich_y's avatar ulrich_y
Browse files

Bug fixes (thanks to A. Gurgone)

parent 735e46cb
......@@ -45,7 +45,7 @@ def Psi(x):
3*trilog(1-x) - 2*Sot(1-x) + np.log(1-x)**3 - np.log(x)**3/6.
+ 3*np.log(x)**2*np.log(1-x)/2. - 3*np.log(x)*np.log(1-x)**2
- 3*dilog(1-x)*np.log(1-x) + 2*pymule.maths.zeta3
- 2*pymule.maths.zeta2 * np.log((1-x) / x)
- 3*pymule.maths.zeta2 * np.log((1-x) / x)
)
......
......@@ -401,9 +401,9 @@ def collintegrand(Ee, o):
xe = 2*Ee/Mmu
logx = np.log(1. - xe + (Mel/Mmu)**2)
if o == 0:
return pref3loop * (arbuzov.f3LL(xe) - 8 * logx**3)
return (2./Mmu) * pref3loop * (arbuzov.f3LL(xe) - 8 * logx**3)
else:
return pref3loop * (arbuzov.g3LL(xe) + 8 * logx**3)
return (2./Mmu) * pref3loop * (arbuzov.g3LL(xe) + 8 * logx**3)
f3LLsub = dointegral(lambda e: collintegrand(e, 0))
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment