diff --git a/include/boost/decimal/detail/concepts.hpp b/include/boost/decimal/detail/concepts.hpp
new file mode 100644
index 000000000..1c9cab5c2
--- /dev/null
+++ b/include/boost/decimal/detail/concepts.hpp
@@ -0,0 +1,380 @@
+// Copyright 2022 - 2023 Matt Borland
+// Distributed under the Boost Software License, Version 1.0.
+// https://www.boost.org/LICENSE_1_0.txt
+
+#ifndef BOOST_DECIMAL_DETAIL_CONCEPTS
+#define BOOST_DECIMAL_DETAIL_CONCEPTS
+
+#if (__cplusplus >= 202002L || _MSVC_LANG >= 202002L) && !defined(BOOST_MATH_DISABLE_CONCEPTS)
+
+#if __has_include(<concepts>)
+
+#include <utility>
+#include <algorithm>
+#include <concepts>
+#include <functional>
+#include <type_traits>
+#include <limits>
+#include <iterator>
+#include <complex>
+#include <boost/decimal/detail/type_traits.hpp>
+
+namespace boost::decimal::concepts {
+
+namespace detail {
+
+template <typename X, typename Y, typename Op>
+struct op_valid_impl
+{
+    template <typename U, typename L, typename R>
+    static constexpr auto test(int) -> decltype(std::declval<U>()(std::declval<L>(), std::declval<R>()),
+                                                void(), std::true_type());
+
+    template <typename U, typename L, typename R>
+    static constexpr auto test(...) -> std::false_type;
+
+    using type = decltype(test<Op, X, Y>(0));
+};
+
+template <typename X, typename Y, typename Op>
+using op_valid_t = typename op_valid_impl<X, Y, Op>::type;
+
+template <typename X, typename Y, typename Op>
+inline constexpr bool op_valid_v = op_valid_t<X, Y, Op>::value;
+
+// Detector for class member functions
+struct nonesuch
+{
+    nonesuch(const nonesuch&) = delete;
+    ~nonesuch() = delete;
+    void operator=(const nonesuch&) = delete;
+};
+
+template <typename Default, typename AlwaysVoid, template<typename...> typename Op, typename... Args>
+struct detector
+{
+    using value_t = std::false_type;
+    using type = Default;
+};
+
+template <typename Default, template<typename...> typename Op, typename... Args>
+struct detector<Default, std::void_t<Op<Args...>>, Op, Args...>
+{
+    using value_t = std::true_type;
+    using type = Op<Args...>;
+};
+
+template <template<typename...> typename Op, typename... Args>
+using is_detected = typename detector<nonesuch, void, Op, Args...>::value_t;
+
+template <template<typename...> typename Op, typename... Args>
+using detected_t = typename detector<nonesuch, void, Op, Args...>::type;
+
+#define BOOST_DECIMAL_HAS_MEMBER_FUNCTION(member)                                      \
+template <typename T>                                                               \
+using has_##member##_t = decltype(std::declval<T&>().member());                     \
+                                                                                    \
+template <typename T>                                                               \
+inline constexpr bool has_##member##_v = is_detected<has_##member##_t, T>::value;
+
+BOOST_DECIMAL_HAS_MEMBER_FUNCTION(begin)
+BOOST_DECIMAL_HAS_MEMBER_FUNCTION(end)
+BOOST_DECIMAL_HAS_MEMBER_FUNCTION(real)
+BOOST_DECIMAL_HAS_MEMBER_FUNCTION(imag)
+
+} // Namespace detail
+
+template <typename T>
+concept integral = std::is_integral_v<T>;
+
+template <typename T>
+concept signed_integral = integral<T> && std::is_signed_v<T>;
+
+template <typename T>
+concept unsigned_integral = integral<T> && std::is_unsigned_v<T>;
+
+template <typename T>
+concept real = std::is_floating_point_v<T>;
+
+template <typename T>
+concept complex = std::is_same_v<T, std::complex<float>>
+                  || std::is_same_v<T, std::complex<double>>
+                  #ifndef BOOST_MATH_NO_LONG_DOUBLE_MATH_FUNCTIONS
+                  || std::is_same_v<T, std::complex<long double>>
+                  #endif
+                  ;
+
+template <typename T>
+concept real_or_complex = real<T> || complex<T>;
+
+template <typename T>
+concept arithmetic = std::is_arithmetic_v<T>;
+
+template <typename T>
+concept numerical = arithmetic<T> || real_or_complex<T>;
+
+template <typename T>
+concept signed_arithmetic = arithmetic<T> && std::is_signed_v<T>;
+
+template <typename T>
+concept unsigned_arithmetic = arithmetic<T> && std::is_unsigned_v<T>;
+
+template <typename T>
+concept arbitrary_unsigned_arithmetic_type = unsigned_arithmetic<T> ||
+                                             (detail::op_valid_v<T, T, std::equal_to<>> &&
+                                              detail::op_valid_v<T, T, std::not_equal_to<>> &&
+                                              detail::op_valid_v<T, T, std::greater<>> &&
+                                              detail::op_valid_v<T, T, std::less<>> &&
+                                              detail::op_valid_v<T, T, std::greater_equal<>> &&
+                                              detail::op_valid_v<T, T, std::less_equal<>> &&
+                                              detail::op_valid_v<T, T, std::plus<>> &&
+                                              detail::op_valid_v<T, T, std::minus<>> &&
+                                              detail::op_valid_v<T, T, std::multiplies<>> &&
+                                              detail::op_valid_v<T, T, std::divides<>>);
+
+template <typename T>
+concept arbitrary_signed_arithmetic_type = signed_arithmetic<T> ||
+                                           (arbitrary_unsigned_arithmetic_type<T> &&
+                                            (detail::op_valid_v<T, T, std::negate<>> ||
+                                             std::numeric_limits<T>::is_signed));
+
+template <typename T>
+concept arbitrary_arithmetic_type = arbitrary_unsigned_arithmetic_type<T> ||
+                                    arbitrary_signed_arithmetic_type<T>;
+
+template <typename T>
+concept arbitrary_unsigned_integer_type = arbitrary_unsigned_arithmetic_type<T> &&
+                                           std::numeric_limits<T>::is_integer;
+
+template <typename T>
+concept arbitrary_signed_integer_type = arbitrary_signed_arithmetic_type<T> &&
+                                         std::numeric_limits<T>::is_integer;
+
+template <typename T>
+concept arbitrary_integer_type = arbitrary_unsigned_integer_type<T> ||
+                                  arbitrary_signed_integer_type<T>;
+
+template <typename T>
+concept arbitrary_real_type = arbitrary_arithmetic_type<T> &&
+                               !std::numeric_limits<T>::is_integer;
+
+template <typename T>
+concept arbitrary_complex_type = complex<T> ||
+                                 (detail::has_real_v<T> &&
+                                  detail::has_imag_v<T>);
+
+template <typename T>
+concept arbitrary_real_or_complex_type = arbitrary_real_type<T> ||
+                                         arbitrary_complex_type<T>;
+
+template <typename T>
+concept arbitrary_numerical_type = arbitrary_real_or_complex_type<T> ||
+                                   arbitrary_arithmetic_type<T>;
+
+template <typename Derived, typename Base>
+concept derived_from = std::is_base_of_v<Base, Derived> &&
+                       std::is_convertible_v<const volatile Derived*, const volatile Base*>;
+
+template <typename T>
+concept forward_iterator = derived_from<typename std::iterator_traits<T>::iterator_category, std::forward_iterator_tag>;
+
+template <typename T>
+concept bidirectional_iterator = derived_from<typename std::iterator_traits<T>::iterator_category, std::bidirectional_iterator_tag>;
+
+template <typename T>
+concept random_access_iterator = derived_from<typename std::iterator_traits<T>::iterator_category, std::random_access_iterator_tag>;
+
+template <typename I, typename T>
+concept output_iterator = derived_from<typename std::iterator_traits<I>::iterator_category, std::input_iterator_tag> &&
+                          derived_from<typename std::iterator_traits<T>::iterator_category, std::output_iterator_tag>;
+
+template <typename T>
+concept is_container = detail::has_begin_v<T> &&
+                       detail::has_end_v<T>;
+
+template <typename T>
+concept random_access_container = is_container<T> &&
+                                  boost::decimal::concepts::random_access_iterator<typename T::iterator>;
+
+template <typename T>
+concept decimal_floating_point_type = detail::is_decimal_floating_point_v<T>;
+
+} // boost::decimal::concepts
+
+#define BOOST_DECIMAL_INTEGRAL boost::decimal::concepts::integral
+#define BOOST_DECIMAL_SIGNED_INTEGRAL boost::decimal::concepts::signed_integral
+#define BOOST_DECIMAL_UNSIGNED_INTEGRAL boost::decimal::concepts::unsigned_integral
+#define BOOST_DECIMAL_REAL boost::decimal::concepts::real
+#define BOOST_DECIMAL_COMPLEX boost::decimal::concepts::complex
+#define BOOST_DECIMAL_REAL_OR_COMPLEX boost::decimal::concepts::real_or_complex
+#define BOOST_DECIMAL_ARITHMETIC boost::decimal::concepts::arithmetic
+#define BOOST_DECIMAL_NUMERICAL boost::decimal::concepts::numerical
+#define BOOST_DECIMAL_SIGNED_ARITHMETIC boost::decimal::concepts::signed_arithmetic
+#define BOOST_DECIMAL_UNSIGNED_ARITHMETIC boost::decimal::concepts::unsigned_arithmetic
+#define BOOST_DECIMAL_ARBITRARY_UNSIGNED_ARITHMETIC boost::decimal::concepts::arbitrary_unsigned_arithmetic_type
+#define BOOST_DECIMAL_ARBITRARY_SIGNED_ARITHMETIC boost::decimal::concepts::arbitrary_signed_arithmetic_type
+#define BOOST_DECIMAL_ARBITRARY_ARITHMETIC boost::decimal::concepts::arbitrary_arithmetic_type
+#define BOOST_DECIMAL_ARBITRARY_UNSIGNED_INTEGER boost::decimal::concepts::arbitrary_unsigned_integer_type
+#define BOOST_DECIMAL_ARBITRARY_SIGNED_INTEGER boost::decimal::concepts::arbitrary_signed_integer_type
+#define BOOST_DECIMAL_ARBITRARY_INTEGER boost::decimal::concepts::arbitrary_integer_type
+#define BOOST_DECIMAL_ARBITRARY_REAL boost::decimal::concepts::arbitrary_real_type
+#define BOOST_DECIMAL_ARBITRARY_COMPLEX boost::decimal::concepts::arbitrary_complex_type
+#define BOOST_DECIMAL_ARBITRARY_REAL_OR_COMPLEX boost::decimal::concepts::arbitrary_real_or_complex_type
+#define BOOST_DECIMAL_ARBITRARY_NUMERICAL boost::decimal::concepts::arbitrary_numerical_type
+#define BOOST_DECIMAL_DECIMAL_FLOATING_TYPE boost::decimal::concepts::decimal_floating_point_type
+
+#define BOOST_DECIMAL_CONTAINER boost::decimal::concepts::is_container
+#define BOOST_DECIMAL_RANDOM_ACCESS_CONTAINER boost::decimal::concepts::random_access_container
+
+#define BOOST_DECIMAL_FORWARD_ITER boost::decimal::concepts::forward_iterator
+#define BOOST_DECIMAL_BIDIRECTIONAL_ITER boost::decimal::concepts::bidirectional_iterator
+#define BOOST_DECIMAL_RANDOM_ACCESS_ITER boost::decimal::concepts::random_access_iterator
+#define BOOST_DECIMAL_OUTPUT_ITER(I, T) boost::decimal::concepts::output_iterator<I, T>
+#define BOOST_DECIMAL_REQUIRES_ITER(X) requires X
+
+#define BOOST_DECIMAL_REQUIRES(X, T) requires X<T>
+
+#ifdef BOOST_DECIMAL_EXEC_COMPATIBLE
+#include <execution>
+
+namespace boost::decimal::concepts {
+
+template <typename T>
+concept execution_policy = std::is_execution_policy_v<std::remove_cvref_t<T>>;
+
+} // Namespace boost::decimal::concepts
+
+#define BOOST_DECIMAL_EXECUTION_POLICY boost::decimal::concepts::execution_policy
+
+#endif // Has <execution>
+
+#endif // Has <concepts>
+#endif // C++20
+
+// If concepts are unavailable replace them with typename for compatibility
+
+#ifndef BOOST_DECIMAL_INTEGRAL
+#  define BOOST_DECIMAL_INTEGRAL typename
+#endif
+
+#ifndef BOOST_DECIMAL_SIGNED_INTEGRAL
+#  define BOOST_DECIMAL_SIGNED_INTEGRAL typename
+#endif
+
+#ifndef BOOST_DECIMAL_UNSIGNED_INTEGRAL
+#  define BOOST_DECIMAL_UNSIGNED_INTEGRAL typename
+#endif
+
+#ifndef BOOST_DECIMAL_REAL
+#  define BOOST_DECIMAL_REAL typename
+#endif
+
+#ifndef BOOST_DECIMAL_COMPLEX
+#  define BOOST_DECIMAL_COMPLEX typename
+#endif
+
+#ifndef BOOST_DECIMAL_REAL_OR_COMPLEX
+#  define BOOST_DECIMAL_REAL_OR_COMPLEX typename
+#endif
+
+#ifndef BOOST_DECIMAL_ARITHMETIC
+#  define BOOST_DECIMAL_ARITHMETIC typename
+#endif
+
+#ifndef BOOST_DECIMAL_NUMERICAL
+#  define BOOST_DECIMAL_NUMERICAL typename
+#endif
+
+#ifndef BOOST_DECIMAL_SIGNED_ARITHMETIC
+#  define BOOST_DECIMAL_SIGNED_ARITHMETIC typename
+#endif
+
+#ifndef BOOST_DECIMAL_UNSIGNED_ARITHMETIC
+#  define BOOST_DECIMAL_UNSIGNED_ARITHMETIC typename
+#endif
+
+#ifndef BOOST_DECIMAL_ARBITRARY_UNSIGNED_ARITHMETIC
+#  define BOOST_DECIMAL_ARBITRARY_UNSIGNED_ARITHMETIC typename
+#endif
+
+#ifndef BOOST_DECIMAL_ARBITRARY_SIGNED_ARITHMETIC
+#  define BOOST_DECIMAL_ARBITRARY_SIGNED_ARITHMETIC typename
+#endif
+
+#ifndef BOOST_DECIMAL_ARBITRARY_ARITHMETIC
+#  define BOOST_DECIMAL_ARBITRARY_ARITHMETIC typename
+#endif
+
+#ifndef BOOST_DECIMAL_ARBITRARY_UNSIGNED_INTEGER
+#  define BOOST_DECIMAL_ARBITRARY_UNSIGNED_INTEGER typename
+#endif
+
+#ifndef BOOST_DECIMAL_ARBITRARY_SIGNED_INTEGER
+#  define BOOST_DECIMAL_ARBITRARY_SIGNED_INTEGER typename
+#endif
+
+#ifndef BOOST_DECIMAL_ARBITRARY_INTEGER
+#  define BOOST_DECIMAL_ARBITRARY_INTEGER typename
+#endif
+
+#ifndef BOOST_DECIMAL_ARBITRARY_REAL
+#  define BOOST_DECIMAL_ARBITRARY_REAL typename
+#endif
+
+#ifndef BOOST_DECIMAL_ARBITRARY_COMPLEX
+#  define BOOST_DECIMAL_ARBITRARY_COMPLEX typename
+#endif
+
+#ifndef BOOST_DECIMAL_ARBITRARY_REAL_OR_COMPLEX
+#  define BOOST_DECIMAL_ARBITRARY_REAL_OR_COMPLEX typename
+#endif
+
+#ifndef BOOST_DECIMAL_ARBITRARY_NUMERICAL
+#  define BOOST_DECIMAL_ARBITRARY_NUMERICAL typename
+#endif
+
+#ifndef BOOST_DECIMAL_POLICY
+#  define BOOST_DECIMAL_POLICY typename
+#endif
+
+#ifndef BOOST_DECIMAL_FORWARD_ITER
+#  define BOOST_DECIMAL_FORWARD_ITER typename
+#endif
+
+#ifndef BOOST_DECIMAL_BIDIRECTIONAL_ITER
+#  define BOOST_DECIMAL_BIDIRECTIONAL_ITER typename
+#endif
+
+#ifndef BOOST_DECIMAL_RANDOM_ACCESS_ITER
+#  define BOOST_DECIMAL_RANDOM_ACCESS_ITER typename
+#endif
+
+#ifndef BOOST_DECIMAL_DECIMAL_FLOATING_TYPE
+#  define BOOST_DECIMAL_DECIMAL_FLOATING_TYPE typename
+#endif
+
+#ifndef BOOST_DECIMAL_OUTPUT_ITER
+#  define BOOST_DECIMAL_OUTPUT_ITER(I, T)
+#endif
+
+#ifndef BOOST_DECIMAL_REQUIRES_ITER
+#  define BOOST_DECIMAL_REQUIRES_ITER(X)
+#endif
+
+#ifndef BOOST_DECIMAL_CONTAINER
+#  define BOOST_DECIMAL_CONTAINER typename
+#endif
+
+#ifndef BOOST_DECIMAL_RANDOM_ACCESS_CONTAINER
+#  define BOOST_DECIMAL_RANDOM_ACCESS_CONTAINER typename
+#endif
+
+#ifndef BOOST_DECIMAL_EXECUTION_POLICY
+#  define BOOST_DECIMAL_EXECUTION_POLICY typename
+#endif
+
+#ifndef BOOST_DECIMAL_REQUIRES
+#  define BOOST_DECIMAL_REQUIRES(X, T)
+#endif
+
+#endif //BOOST_DECIMAL_DETAIL_CONCEPTS