Commit c74ee500 authored by Luca Naterop's avatar Luca Naterop

get simply a log if two args but not convergent

parent 5fd2633a
...@@ -104,9 +104,17 @@ CONTAINS ...@@ -104,9 +104,17 @@ CONTAINS
call print_G(z_flat,y) call print_G(z_flat,y)
! only two arguments? -> we get a logarithm
! need make convergent? ! need make convergent?
if(.not. is_convergent(z_flat,y)) then if(.not. is_convergent(z_flat,y)) then
print*, 'need to make convergent' print*, 'need to make convergent'
if(size(z_flat) == 1) then
print*, 'we get simply a log'
res = log(y-z_flat(1)) - log(-z_flat(1))
return
end if
res = reduce_to_convergent(z_flat, y) res = reduce_to_convergent(z_flat, y)
return return
end if end if
......
...@@ -18,7 +18,7 @@ PROGRAM TEST ...@@ -18,7 +18,7 @@ PROGRAM TEST
! call do_GPL_tests() ! call do_GPL_tests()
! call do_shuffle_tests() ! put this somewhere else ! call do_shuffle_tests() ! put this somewhere else
res = G_flat(cmplx((/1.0,10.0/)),cmplx(2.0)) res = G_flat(cmplx((/1.0,3.0/)),cmplx(2.0))
print*, res print*, res
! if(tests_successful) then ! if(tests_successful) then
......
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