Commit a360ebc8 authored by ulrich_y's avatar ulrich_y

Bug fix re 9174f364 (GPL_zero_zi had the wrong branch-cut for neg. arguments)

parent 0672aaaa
......@@ -18,11 +18,12 @@ CONTAINS
integer :: l
type(inum) :: y
complex(kind=prec) :: GPL_zero_zi
!TODO
if (abs(aimag(y)).lt.zero) then
if (real(y).gt.0) then
GPL_zero_zi = 1.0_prec/factorial(l) * log(real(y)) ** l
else
GPL_zero_zi = 1.0_prec/factorial(l) * (log(-real(y))-cmplx(0,y%i0*pi)) ** l
GPL_zero_zi = 1.0_prec/factorial(l) * (log(-real(y))+cmplx(0,y%i0*pi)) ** l
endif
else
GPL_zero_zi = 1.0_prec/factorial(l) * log(y%c) ** l
......
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