diff --git a/include/boost/math/special_functions/detail/bessel_ik.hpp b/include/boost/math/special_functions/detail/bessel_ik.hpp
index 4db8f14a7f..1c9be54a36 100644
--- a/include/boost/math/special_functions/detail/bessel_ik.hpp
+++ b/include/boost/math/special_functions/detail/bessel_ik.hpp
@@ -327,7 +327,7 @@ int bessel_ik(T v, T x, T* result_I, T* result_K, int kind, const Policy& pol)
     BOOST_MATH_INSTRUMENT_VARIABLE(n);
     BOOST_MATH_INSTRUMENT_VARIABLE(u);
 
-    BOOST_ASSERT(x > 0); // Error handling for x <= 0 handled in cyl_bessel_i and cyl_bessel_k
+    BOOST_MATH_ASSERT(x > 0); // Error handling for x <= 0 handled in cyl_bessel_i and cyl_bessel_k
 
     // x is positive until reflection
     W = 1 / x;                                 // Wronskian
diff --git a/include/boost/math/special_functions/detail/bessel_j0.hpp b/include/boost/math/special_functions/detail/bessel_j0.hpp
index 2a40407170..9a0b26fe6b 100644
--- a/include/boost/math/special_functions/detail/bessel_j0.hpp
+++ b/include/boost/math/special_functions/detail/bessel_j0.hpp
@@ -131,7 +131,7 @@ T bessel_j0(T x)
     using namespace boost::math::tools;
     using namespace boost::math::constants;
 
-    BOOST_ASSERT(x >= 0); // reflection handled elsewhere.
+    BOOST_MATH_ASSERT(x >= 0); // reflection handled elsewhere.
 
     if (x == 0)
     {
diff --git a/include/boost/math/special_functions/detail/bessel_j1.hpp b/include/boost/math/special_functions/detail/bessel_j1.hpp
index 696dbc83d4..6d354dcce7 100644
--- a/include/boost/math/special_functions/detail/bessel_j1.hpp
+++ b/include/boost/math/special_functions/detail/bessel_j1.hpp
@@ -196,7 +196,7 @@ T bessel_j1(T x)
         value = factor * (rc * (sx - cx) + y * rs * (sx + cx));
     }
 
-    BOOST_ASSERT(x >= 0);  // Negative values handled by the caller.
+    BOOST_MATH_ASSERT(x >= 0);  // Negative values handled by the caller.
 
     return value;
 }