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;