From 445d48cfbf0305567b709af659e875ebf4f0cabe Mon Sep 17 00:00:00 2001 From: Andreas Suter Date: Tue, 5 Mar 2024 18:29:34 +0100 Subject: [PATCH] revert back to boost ellint_1, since on many compilers the c++17 std::comp_ellint_1 is not yet implemented. --- src/external/libBNMR/libLineProfile/LineProfile.cpp | 6 ++++-- 1 file changed, 4 insertions(+), 2 deletions(-) diff --git a/src/external/libBNMR/libLineProfile/LineProfile.cpp b/src/external/libBNMR/libLineProfile/LineProfile.cpp index eca2aa61..486198bb 100644 --- a/src/external/libBNMR/libLineProfile/LineProfile.cpp +++ b/src/external/libBNMR/libLineProfile/LineProfile.cpp @@ -23,6 +23,8 @@ #include "LineProfile.h" #include //for testing purposes +#include + //Implement helperfunctions Double_t GaussianShape(Double_t x, Double_t position, Double_t width) { @@ -80,11 +82,11 @@ Double_t IAsym(Double_t x, Double_t omega_center, Double_t omega_min,Double_t om Double_t IAsym_low(Double_t omega, Double_t omega_center, Double_t omega_min,Double_t omega_max){ Double_t m =(omega-omega_min)*(omega_max-omega_center)/((omega_max-omega)*(omega_center-omega_min)); - return 1./PI/pow(omega_max-omega,0.5)/pow(omega_center-omega_min,0.5)*std::comp_ellint_1( m ); + return 1./PI/pow(omega_max-omega,0.5)/pow(omega_center-omega_min,0.5)*boost::math::ellint_1( m ); } Double_t IAsym_high(Double_t omega, Double_t omega_center, Double_t omega_min,Double_t omega_max){ Double_t m =(omega_center-omega_min)*(omega_max-omega)/((omega_max-omega_center)*(omega-omega_min)); - return 1./PI/pow(omega-omega_min,0.5)/pow(omega_max-omega_center,0.5)*std::comp_ellint_1( m ); + return 1./PI/pow(omega-omega_min,0.5)/pow(omega_max-omega_center,0.5)*boost::math::ellint_1( m ); }