Use correct cut-off for polylog

src/globals.f90

@@ -16,7 +16,7 @@ MODULE globals
 
   ! 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 ::  Lidelta = 0._prec       ! like MPLdelta but for polylogs
+  real(kind=prec), protected ::  Lidelta = zero          ! like MPLdelta but for polylogs
   real(kind=prec), protected :: HoelderCircle = 1.1_prec ! when to apply Hoelder convolution?
   integer, parameter :: PolyLogCacheSize(2) = (/ 5, 100 /)
         ! = (/ mmax, n /). At most n polylogs with weight mmax will be cached

src/maths_functions.f90

@@ -183,7 +183,7 @@ CONTAINS
     ! All remaining terms
     m = n + 1
     del = 1._prec
-    do while (abs(del) > zero)
+    do while (abs(del) > LiDelta)
       zetamn = (-1)**(m-n) * bernoullinumber(1+m-n) / (1+m-n)
       fac = fac * m * (m-1)
       del = zetamn / fac * logz**m