Commit 71d3319a authored by Luca Naterop's avatar Luca Naterop

minor

parent ae79dbd1
......@@ -160,7 +160,7 @@ CONTAINS
if(verb >= 30) print*, 'case higher depth, sr at beginning'
if(verb >= 50) then
print*, 'map to (using 68)'
print*, 'map to (using 21)'
print*, 'PI with p=',real(p),'i=',0,'g =',tocmplx(zeroes(0) )
call print_G([izero,a],y2)
print*, 'PI with p=',real([p,y2]),'i=',0,'g =',tocmplx(zeroes(0) )
......@@ -193,11 +193,11 @@ CONTAINS
if(verb >= 30) print*, 's_r in the middle under PI'
#endif
res = +pending_integral(p,1,zeroes(0)) * G_flat([a(1:i-1),izero,a(i:size(a))],y2) & !64
- pending_integral([p,a(i-1)],i-1,[a(1:i-2),a(i:size(a)),y2]) & ! 67a
+ pending_integral([p,a(i-1)],1,zeroes(0)) * G_flat(a,y2) & ! 67c
+ pending_integral([p,a(i)], i, [a(1:i-1), a(i+1:size(a)),y2]) & !67b
- pending_integral([p,a(i)],1,zeroes(0)) * G_flat(a,y2) !67a
res = +pending_integral(p,1,zeroes(0)) * G_flat([a(1:i-1),izero,a(i:size(a))],y2) &
- pending_integral([p,a(i-1)],i-1,[a(1:i-2),a(i:size(a)),y2]) &
+ pending_integral([p,a(i-1)],1,zeroes(0)) * G_flat(a,y2) &
+ pending_integral([p,a(i)], i, [a(1:i-1), a(i+1:size(a)),y2]) &
- pending_integral([p,a(i)],1,zeroes(0)) * G_flat(a,y2)
END FUNCTION pending_integral
RECURSIVE FUNCTION remove_sr_from_last_place_in_G(a,y2,m,sr) result(res)
......
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