diff --git a/src/external/libBNMR/libLineProfile/LineProfile.cpp b/src/external/libBNMR/libLineProfile/LineProfile.cpp
index eca2aa616..486198bb4 100644
--- a/src/external/libBNMR/libLineProfile/LineProfile.cpp
+++ b/src/external/libBNMR/libLineProfile/LineProfile.cpp
@@ -23,6 +23,8 @@
 #include "LineProfile.h"
 #include <iostream> //for testing purposes
 
+#include <boost/math/special_functions/ellint_1.hpp>
+
 //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 );
     }