From 9452013782815966c6a6f3c0504e878ea1ea37eb Mon Sep 17 00:00:00 2001 From: Thomas Applencourt Date: Tue, 5 May 2015 20:58:41 +0200 Subject: [PATCH] More clarity --- src/Pseudo_integrals/int.f90 | 12 +++--------- 1 file changed, 3 insertions(+), 9 deletions(-) diff --git a/src/Pseudo_integrals/int.f90 b/src/Pseudo_integrals/int.f90 index 6f0d4988..bade7eb2 100644 --- a/src/Pseudo_integrals/int.f90 +++ b/src/Pseudo_integrals/int.f90 @@ -1828,12 +1828,6 @@ double precision function coef_nk(n,k) double precision gam,dble_fact,fact double precision, save :: store_result(0:2,0:10) - integer, save :: ifirst - - if (ifirst == 0) then - ifirst =1 - - endif store_result(0, 0) = 1.00000000000000d0 store_result(0, 1) = 0.166666666666667d0 @@ -1878,7 +1872,7 @@ double precision function coef_nk(n,k) endif endif - if ((n+k).GE.30) then + if ((n+k).GE.60) then coef_nk = 0.d0 return endif @@ -1945,10 +1939,10 @@ double precision function int_prod_bessel(l,gam,n,m,a,b,arg) sum=0.d0 temp_array_a(q) = coef_nk(n,q)*(a**(n+2*q)) - temp_array_b(m+2*q) = coef_nk(m,q)*b**(m+2*q) + temp_array_b(q) = coef_nk(m,q)*b**(m+2*q) do k=0,q - sum=sum+temp_array_a(k)*temp_array_b(m-2*k+2*q) + sum=sum+temp_array_a(k)*temp_array_b(q-2*k) enddo int=int+sum*crochet(2*q+n+m+l,gam)