Function neg instead of imone

src/gpl_module.f90

@@ -103,7 +103,7 @@ CONTAINS
       !res = pending_integral(p,2,[sub_ieps(g(1))]) - pending_integral(p,2,[cmplx(0.0)]) &
       !  + G_flat(p(2:size(p)), p(1)) * log(-sub_ieps(g(1)))
       res = pending_integral(p,2,[g(1)]) - pending_integral(p,2,[izero]) &
-        + G_flat(p(2:size(p)), p(1)) * log(imone*g(1))
+        + G_flat(p(2:size(p)), p(1)) * log(neg(g(1)))
       return
     end if
   

src/ieps.f90

@@ -48,6 +48,12 @@ MODULE ieps
   end interface real
 CONTAINS
 
+  FUNCTION NEG(n1)
+  implicit none
+  type(inum), intent(in) :: n1
+  type(inum) :: neg
+  neg = inum(-n1%c,-n1%i0)
+  END FUNCTION
 
   FUNCTION MULTINUMSS(n1, n2)
   implicit none

src/maths_functions.f90

@@ -312,9 +312,8 @@ CONTAINS
       res = -(1. - 2.**(1-m))*zeta(m)
       return
     else if (abs(x) .gt. 1) then
-      print*,imone*x
       res = (-1)**(m-1)*polylog(m,ione/x) &
-          - cmplx(0,2*pi)**m * bernoulli_polynomial(m, 0.5-cmplx(0.,1.)*log(imone*x)/2/pi) / factorial(m)
+          - cmplx(0,2*pi)**m * bernoulli_polynomial(m, 0.5-cmplx(0.,1.)*log(neg(x))/2/pi) / factorial(m)
       return
     endif