Correct compilation in bessel boost version

Inria-Empenn · Jul 13, 2018 · 2482cff · 2482cff
1 parent 0a7e99c
commit 2482cff
Showing 1 changed file with 4 additions and 1 deletion.
diff --git a/Anima/math-tools/special_functions/animaBesselFunctions.cxx b/Anima/math-tools/special_functions/animaBesselFunctions.cxx
@@ -46,7 +46,10 @@ double scaled_bessel_i(unsigned int N, double x)
 #else
 
     boost::multiprecision::cpp_dec_float_100 y = x;
-    return std::exp(-y) * boost::math::cyl_bessel_i<boost::multiprecision::cpp_dec_float_100>(N, y);
+    boost::multiprecision::cpp_dec_float_100 expX = std::exp(-x);
+
+    boost::multiprecision::cpp_dec_float_100 resVal = expX * boost::math::cyl_bessel_i<boost::multiprecision::cpp_dec_float_100>(N, y);
+    return resVal.convert_to<double>();
 
 #endif
 }