LCPQ/quantum_package
10
0
Fork 0
mirror of https://github.com/LCPQ/quantum_package synced 2024-12-23 04:43:50 +01:00
Code Issues Releases Wiki Activity

solve overflow error in int.f90 for the pseudo

Browse Source
This commit is contained in:
Thomas Applencourt 2015-05-05 11:56:37 +02:00
parent 8dc7ea8481
commit eeb60de46d
1 changed files with 12 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

19
src/Pseudo_integrals/int.f90
View File

@ -1879,13 +1879,18 @@ end
RETURN
END
double precision function coef_nk(n,k)
implicit none
integer n,k
double precision gam,dblefact,fact
gam=dblefact(2*(n+k)+1)
coef_nk=1.d0/(2.d0**k*fact(k)*gam)
end
double precision function coef_nk(n,k)
implicit none
integer n,k
double precision gam,dblefact,fact
if(k.GE.80) then
coef_nk = 0.d0
else
gam=dblefact(2*(n+k)+1)
coef_nk=1.d0/(2.d0**k*fact(k)*gam)
endif
end
!! Calculation of
!!