Commit cd3ea4cf authored by ulrich_y's avatar ulrich_y

Reactivated LiDelta

parent 39767146
...@@ -16,7 +16,7 @@ MODULE globals ...@@ -16,7 +16,7 @@ MODULE globals
! The following parameters control the accuracy of the evaluation ! The following parameters control the accuracy of the evaluation
real(kind=prec), protected :: MPLdelta = zero ! if the MPL sum changes less then del it is truncated. real(kind=prec), protected :: MPLdelta = zero ! if the MPL sum changes less then del it is truncated.
real(kind=prec), protected :: Lidelta = zero ! like MPLdelta but for polylogs real(kind=prec), protected :: Lidelta = 0._prec ! like MPLdelta but for polylogs
real(kind=prec), protected :: HoelderCircle = 1.1_prec ! when to apply Hoelder convolution? real(kind=prec), protected :: HoelderCircle = 1.1_prec ! when to apply Hoelder convolution?
integer, parameter :: PolyLogCacheSize(2) = (/ 5, 100 /) integer, parameter :: PolyLogCacheSize(2) = (/ 5, 100 /)
! = (/ mmax, n /). At most n polylogs with weight mmax will be cached ! = (/ mmax, n /). At most n polylogs with weight mmax will be cached
......
...@@ -41,6 +41,7 @@ CONTAINS ...@@ -41,6 +41,7 @@ CONTAINS
if(i**m.lt.0) return ! roll over if(i**m.lt.0) return ! roll over
if(abs(x**i).lt.1.e-250_prec) return if(abs(x**i).lt.1.e-250_prec) return
del = x**i/i**m del = x**i/i**m
if (abs(del) < LiDelta) return
res = res+del res = res+del
i = i+1 i = i+1
end do end do
......
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