File o2scl-exp-overflow.patch of Package o2scl

Overview Repositories Revisions Requests Users Attributes Meta

File o2scl-exp-overflow.patch of Package o2scl (Revision e5f7dc2d2893e7fae4a6b11598f47caf)

Currently displaying revision e5f7dc2d2893e7fae4a6b11598f47caf , Show latest

From 78c3ce2c44030059e208c62c62217e594eca02ff Mon Sep 17 00:00:00 2001
From: "Andrew W. Steiner" <awsteiner@mykolab.com>
Date: Fri, 13 Nov 2020 22:23:34 -0800
Subject: [PATCH] Fix for overflows from GSL exp.c on opensuse.

---
 src/base/prev_commit.h |  4 ++--
 src/other/polylog.h    | 43 ++++++++++++++++++++++++++++++++++++------
 2 files changed, 39 insertions(+), 8 deletions(-)

diff --git a/src/other/polylog.h b/src/other/polylog.h
index 55c8cd6d..097f2e85 100644
--- a/src/other/polylog.h
+++ b/src/other/polylog.h
@@ -43,6 +43,7 @@
 #include <o2scl/lib_settings.h>
 #include <o2scl/inte_adapt_cern.h>
 #include <o2scl/inte_double_exp_boost.h>
+#include <o2scl/exception.h>
 
 #ifndef DOXYGEN_NO_O2NS
 namespace o2scl {
@@ -229,31 +230,61 @@ namespace o2scl {
     /** \brief Fermi-Dirac integral of order \f$ 1/2 \f$
      */
     double calc_1o2(double y) {
-      return gsl_sf_fermi_dirac_half(y)*sqrt(o2scl_const::pi)/2.0;
+      double ret;
+      try {
+	ret=gsl_sf_fermi_dirac_half(y)*sqrt(o2scl_const::pi)/2.0;
+      } catch (const o2scl::exc_overflow_error &e) {
+	return std::numeric_limits<double>::infinity();
+      }
+      return ret;
     }
     
     /** \brief Fermi-Dirac integral of order \f$ -1/2 \f$
      */
     double calc_m1o2(double y) {
-      return gsl_sf_fermi_dirac_mhalf(y)*sqrt(o2scl_const::pi);
+      double ret;
+      try {
+	ret=gsl_sf_fermi_dirac_mhalf(y)*sqrt(o2scl_const::pi);
+      } catch (const o2scl::exc_overflow_error &e) {
+	return std::numeric_limits<double>::infinity();
+      }
+      return ret;
     }
     
     /** \brief Fermi-Dirac integral of order \f$ 3/2 \f$
      */
     double calc_3o2(double y) {
-      return gsl_sf_fermi_dirac_3half(y)*sqrt(o2scl_const::pi)*0.75;
+      double ret;
+      try {
+	ret=gsl_sf_fermi_dirac_3half(y)*sqrt(o2scl_const::pi)*0.75;
+      } catch (const o2scl::exc_overflow_error &e) {
+	return std::numeric_limits<double>::infinity();
+      }
+      return ret;
     }
-
+    
     /** \brief Fermi-Dirac integral of order \f$ 2 \f$
      */
     double calc_2(double y) {
-      return gsl_sf_fermi_dirac_int(2,y)*2.0;
+      double ret;
+      try {
+	ret=gsl_sf_fermi_dirac_int(2,y)*2.0;
+      } catch (const o2scl::exc_overflow_error &e) {
+	return std::numeric_limits<double>::infinity();
+      }
+      return ret;
     }
     
     /** \brief Fermi-Dirac integral of order \f$ 3 \f$
      */
     double calc_3(double y) {
-      return gsl_sf_fermi_dirac_int(3,y)*6.0;
+      double ret;
+      try {
+	ret=gsl_sf_fermi_dirac_int(3,y)*6.0;
+      } catch (const o2scl::exc_overflow_error &e) {
+	return std::numeric_limits<double>::infinity();
+      }
+      return ret;
     }
     
   };

Places

File o2scl-exp-overflow.patch of Package o2scl (Revision e5f7dc2d2893e7fae4a6b11598f47caf)

Places