Commit 9c315951 authored by ulrich_y's avatar ulrich_y
Browse files

Minor bug-fixes

parent e360aa6a
...@@ -262,7 +262,7 @@ CONTAINS ...@@ -262,7 +262,7 @@ CONTAINS
- G_flat([a(i+1)],sr) * G_flat([a(1:i-1),a(i+1:size(a))],y2) - G_flat([a(i+1)],sr) * G_flat([a(1:i-1),a(i+1:size(a))],y2)
END FUNCTION make_convergent END FUNCTION make_convergent
FUNCTION improve_convergence(z) result(res) RECURSIVE FUNCTION improve_convergence(z) result(res)
! improves the convergence by applying the Hoelder convolution to G(z1,...zk,1) ! improves the convergence by applying the Hoelder convolution to G(z1,...zk,1)
complex(kind=prec) :: z(:),oneminusz(size(z)), res complex(kind=prec) :: z(:),oneminusz(size(z)), res
complex(kind=prec), parameter :: p = 2.0 complex(kind=prec), parameter :: p = 2.0
...@@ -427,8 +427,9 @@ CONTAINS ...@@ -427,8 +427,9 @@ CONTAINS
return return
end if end if
do i = 1,k x(1) = y/z(1)
x(i) = merge(y/z(1), z(i-1)/z(i),i == 1) do i = 2,k
x(i) = z(i-1)/z(i)
end do end do
! print*, 'computed using Li with ' ! print*, 'computed using Li with '
! print*, 'm = ', m ! print*, 'm = ', m
......
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