From d6542e62a3cc3880b08aa85258467abae7be7c6b Mon Sep 17 00:00:00 2001
From: Anthony Scemama <scemama@irsamc.ups-tlse.fr>
Date: Wed, 17 Nov 2021 09:02:26 +0100
Subject: [PATCH] Fix floating-point exception

---
 src/csf/sigma_vector.irp.f | 7 +++++--
 1 file changed, 5 insertions(+), 2 deletions(-)

diff --git a/src/csf/sigma_vector.irp.f b/src/csf/sigma_vector.irp.f
index 85ed5f84..77f6190b 100644
--- a/src/csf/sigma_vector.irp.f
+++ b/src/csf/sigma_vector.irp.f
@@ -1,9 +1,12 @@
  real*8 function logabsgamma(x)
  implicit none
  real*8, intent(in) :: x
- logabsgamma = log(abs(gamma(x)))
+ logabsgamma = 1.d32  ! Avoid floating point exception
+ if (x>0.d0) then
+   logabsgamma = log(abs(gamma(x)))
+ endif
  end function logabsgamma
-  
+
   BEGIN_PROVIDER [ integer, NSOMOMax]
  &BEGIN_PROVIDER [ integer, NCSFMax]
  &BEGIN_PROVIDER [ integer*8, NMO]