From 57527c94a27491107e49d9623c0124c516976d34 Mon Sep 17 00:00:00 2001
From: ydamour <yann.damour@hotmail.fr>
Date: Thu, 5 May 2022 17:26:41 +0200
Subject: [PATCH] update binom_func to avoid .999...

---
 src/utils/util.irp.f | 4 ++++
 1 file changed, 4 insertions(+)

diff --git a/src/utils/util.irp.f b/src/utils/util.irp.f
index ef846bdb..127d5433 100644
--- a/src/utils/util.irp.f
+++ b/src/utils/util.irp.f
@@ -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