Commit df3ee22c authored by ulrich_y's avatar ulrich_y

Used new polylog function

parent e769b3af
......@@ -142,7 +142,6 @@ CONTAINS
subroutine do_poly_tests()
implicit none
complex(kind=prec) :: ref, z
tol = 1e-9 ! TODO way too loose
z = 0.9999_prec * exp((0.,1._prec))
ref = ( 0.324133535392567026495406418480254855764_prec,1.013852051950427381643871737213288702950_prec)
......
......@@ -271,6 +271,8 @@ CONTAINS
if(abs(aimag(x)) < zero ) then
res = Li2(real(x,kind=prec))
else if ( (0.5_prec .lt. abs(x)) .and. (abs(x) .lt. 2._prec) ) then
res = logz_polylog(2,x)
else
res = naive_polylog(2,x)
endif
......@@ -446,6 +448,8 @@ CONTAINS
complex(kind=prec) :: x
if(abs(aimag(x)) < zero ) then
res = Li3(real(x,kind=prec))
else if ( (0.5_prec .lt. abs(x)) .and. (abs(x) .lt. 2._prec) ) then
res = logz_polylog(3,x)
else
res = naive_polylog(3,x)
endif
......@@ -555,6 +559,8 @@ CONTAINS
res = dilog(x)
else if(m == 3) then
res = trilog(x)
else if ( (0.5_prec .lt. abs(x)) .and. (abs(x) .lt. 2._prec) ) then
res = logz_polylog(m,x)
else
res = naive_polylog(m,x)
end if
......
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