Commit e85716b5 authored by ulrich_y's avatar ulrich_y

Fixed some TODOs

parent f652eea6
......@@ -20,7 +20,6 @@ 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
......@@ -127,9 +126,8 @@ CONTAINS
#endif
!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)))
!TODO
res = pending_integral(p,2,[g(1)]) - pending_integral(p,2,[izero]) &
+ G_flat(p(2:size(p)), p(1)) * (log(g(1)%c) + p(1)%i0 * pi * (0.,1.))
res = pending_integral(p,2,[inum( g(1)%c,-g(1)%i0 ) ]) - pending_integral(p,2,[izero]) &
+ G_flat(p(2:size(p)), p(1)) * (log(g(1)%c) + p(1)%i0 * pi * i_)
return
end if
......@@ -399,9 +397,8 @@ CONTAINS
! requires Hoelder convolution?
if( any(1.0 <= abs(z_flat%c/y%c) .and. abs(z_flat%c/y%c) <= HoelderCircle) ) then
! Here we just *assume* that y is positive and doesn't mess up the
! Here we *assume* that y is positive and doesn't mess up the
! ieps, which is what GiNaC does (l. 1013)
! TODO
do j=1,size(z_flat)
znorm(j) = inum(z_flat(j)%c/y%c, z_flat(j)%i0)
enddo
......@@ -479,7 +476,7 @@ CONTAINS
res = G_flat(z_flat,y)
return
end if
!TODO is that okay?
! This is okay because in the MPLs the ieps doesn't matter anymore
x(1) = y%c/z(1)%c
do i = 2,k
x(i) = z(i-1)%c/z(i)%c
......
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