9
1
mirror of https://github.com/QuantumPackage/qp2.git synced 2024-11-07 22:13:38 +01:00

update binom_func to avoid .999...

This commit is contained in:
Yann Damour 2022-05-05 17:26:41 +02:00
parent c7f120f35e
commit 57527c94a2

View File

@ -37,6 +37,10 @@ double precision function binom_func(i,j)
else
binom_func = dexp( logfact(i)-logfact(j)-logfact(i-j) )
endif
! To avoid .999999 numbers
binom_func = binom_func + 0.5d0
end