scaddins/source/analysis/bessel.cxx |    7 +++----
 1 file changed, 3 insertions(+), 4 deletions(-)

New commits:
commit 2acf3fae5967fb2002621ffa2280ed51311e27e9
Author:     Jesus Solis <sol...@alum.mit.edu>
AuthorDate: Tue Aug 19 19:35:24 2025 -0500
Commit:     David Gilbert <freedesk...@treblig.org>
CommitDate: Tue Aug 26 01:14:23 2025 +0200

    tdf#145759 replace magic constant with symbolic constant
    
    Change-Id: Ic2148533b1e1fdea66c8b6f575474dcad4069419
    Reviewed-on: https://gerrit.libreoffice.org/c/core/+/189940
    Tested-by: Jenkins
    Reviewed-by: David Gilbert <freedesk...@treblig.org>

diff --git a/scaddins/source/analysis/bessel.cxx 
b/scaddins/source/analysis/bessel.cxx
index 44b79e798f85..5120a3747096 100644
--- a/scaddins/source/analysis/bessel.cxx
+++ b/scaddins/source/analysis/bessel.cxx
@@ -19,6 +19,7 @@
 
 #include "bessel.hxx"
 #include <cmath>
+#include <numbers>
 #include <rtl/math.hxx>
 
 #include <com/sun/star/lang/IllegalArgumentException.hpp>
@@ -327,8 +328,7 @@ static double Bessely0( double fX )
         return sqrt(1/M_PI/fX)
                 *(std::sin(fX)-std::cos(fX));
     const double epsilon = 1.0e-15;
-    const double EulerGamma = 0.57721566490153286060;
-    double alpha = log(fX/2.0)+EulerGamma;
+    double alpha = log(fX/2.0)+std::numbers::egamma;
     double u = alpha;
 
     double k = 1.0;
@@ -381,12 +381,11 @@ static double Bessely1( double fX )
         return - sqrt(1/M_PI/fX)
                 *(std::sin(fX)+std::cos(fX));
     const double epsilon = 1.0e-15;
-    const double EulerGamma = 0.57721566490153286060;
     double alpha = 1.0/fX;
     double f_bar = -1.0;
     double u = alpha;
     double k = 1.0;
-    alpha = 1.0 - EulerGamma - log(fX/2.0);
+    alpha = 1.0 - std::numbers::egamma - log(fX/2.0);
     double g_bar_delta_u = -alpha;
     double g_bar = -2.0 / fX;
     double delta_u = g_bar_delta_u / g_bar;

Reply via email to