/////////////////////////////////////////////////////////////// // Copyright 2012 John Maddock. Distributed under the Boost // Software License, Version 1.0. (See accompanying file // LICENSE_1_0.txt or copy at https://www.boost.org/LICENSE_1_0.txt #ifndef BOOST_MP_DEBUG_ADAPTOR_HPP #define BOOST_MP_DEBUG_ADAPTOR_HPP #include #include #include namespace boost { namespace multiprecision { namespace backends { template struct debug_adaptor { using signed_types = typename Backend::signed_types ; using unsigned_types = typename Backend::unsigned_types; using float_types = typename Backend::float_types ; using exponent_type = typename extract_exponent_type::value>::type; private: std::string debug_value; Backend m_value; public: void update_view() { #ifndef BOOST_NO_EXCEPTIONS try { #endif debug_value = m_value.str(0, static_cast(0)); #ifndef BOOST_NO_EXCEPTIONS } catch (const std::exception& e) { debug_value = "String conversion failed with message: \""; debug_value += e.what(); debug_value += "\""; } #endif } debug_adaptor() { update_view(); } debug_adaptor(const debug_adaptor& o) : debug_value(o.debug_value), m_value(o.m_value) { } debug_adaptor& operator=(const debug_adaptor& o) { debug_value = o.debug_value; m_value = o.m_value; return *this; } template debug_adaptor(const T& i, const typename std::enable_if::value>::type* = nullptr) : m_value(i) { update_view(); } template debug_adaptor(const debug_adaptor& i, const typename std::enable_if::value>::type* = nullptr) : m_value(i.value()) { update_view(); } template debug_adaptor(const T& i, const U& j, typename std::enable_if::value>::type* = nullptr) : m_value(i, j) { update_view(); } template debug_adaptor(const B2& i, unsigned digits10, typename std::enable_if::value && std::is_constructible::value>::type* = nullptr) : m_value(i, digits10) { update_view(); } template typename std::enable_if::value || std::is_assignable::value, debug_adaptor&>::type operator=(const T& i) { m_value = i; update_view(); return *this; } debug_adaptor& operator=(const char* s) { m_value = s; update_view(); return *this; } void swap(debug_adaptor& o) { std::swap(m_value, o.value()); std::swap(debug_value, o.debug_value); } std::string str(std::streamsize digits, std::ios_base::fmtflags f) const { return m_value.str(digits, f); } void negate() { m_value.negate(); update_view(); } int compare(const debug_adaptor& o) const { return m_value.compare(o.value()); } template int compare(const T& i) const { return m_value.compare(i); } Backend& value() { return m_value; } const Backend& value() const { return m_value; } #ifndef BOOST_MP_STANDALONE template void serialize(Archive& ar, const unsigned int /*version*/) { ar & boost::make_nvp("value", m_value); using tag = typename Archive::is_loading; if (tag::value) update_view(); } #endif static unsigned default_precision() noexcept { return Backend::default_precision(); } static void default_precision(unsigned v) noexcept { Backend::default_precision(v); } static unsigned thread_default_precision() noexcept { return Backend::thread_default_precision(); } static void thread_default_precision(unsigned v) noexcept { Backend::thread_default_precision(v); } unsigned precision() const noexcept { return value().precision(); } void precision(unsigned digits10) noexcept { value().precision(digits10); } // // Variable precision options: // static constexpr variable_precision_options default_variable_precision_options()noexcept { return Backend::default_variable_precision_options(); } static constexpr variable_precision_options thread_default_variable_precision_options()noexcept { return Backend::thread_default_variable_precision_options(); } static BOOST_MP_CXX14_CONSTEXPR void default_variable_precision_options(variable_precision_options opts) { Backend::default_variable_precision_options(opts); } static BOOST_MP_CXX14_CONSTEXPR void thread_default_variable_precision_options(variable_precision_options opts) { Backend::thread_default_variable_precision_options(opts); } }; template inline Backend const& unwrap_debug_type(debug_adaptor const& val) { return val.value(); } template inline const T& unwrap_debug_type(const T& val) { return val; } template inline BOOST_MP_CXX14_CONSTEXPR void assign_components(debug_adaptor& result, const V& v1, const U& v2) { using default_ops::assign_components; assign_components(result.value(), unwrap_debug_type(v1), unwrap_debug_type(v2)); result.update_view(); } #define NON_MEMBER_OP1(name, str) \ template \ inline void BOOST_JOIN(eval_, name)(debug_adaptor & result) \ { \ using default_ops::BOOST_JOIN(eval_, name); \ BOOST_JOIN(eval_, name) \ (result.value()); \ result.update_view(); \ } #define NON_MEMBER_OP2(name, str) \ template \ inline void BOOST_JOIN(eval_, name)(debug_adaptor & result, const T& a) \ { \ using default_ops::BOOST_JOIN(eval_, name); \ BOOST_JOIN(eval_, name) \ (result.value(), unwrap_debug_type(a)); \ result.update_view(); \ } \ template \ inline void BOOST_JOIN(eval_, name)(debug_adaptor & result, const debug_adaptor& a) \ { \ using default_ops::BOOST_JOIN(eval_, name); \ BOOST_JOIN(eval_, name) \ (result.value(), unwrap_debug_type(a)); \ result.update_view(); \ } #define NON_MEMBER_OP3(name, str) \ template \ inline void BOOST_JOIN(eval_, name)(debug_adaptor & result, const T& a, const U& b) \ { \ using default_ops::BOOST_JOIN(eval_, name); \ BOOST_JOIN(eval_, name) \ (result.value(), unwrap_debug_type(a), unwrap_debug_type(b)); \ result.update_view(); \ } \ template \ inline void BOOST_JOIN(eval_, name)(debug_adaptor & result, const debug_adaptor& a, const T& b) \ { \ using default_ops::BOOST_JOIN(eval_, name); \ BOOST_JOIN(eval_, name) \ (result.value(), unwrap_debug_type(a), unwrap_debug_type(b)); \ result.update_view(); \ } \ template \ inline void BOOST_JOIN(eval_, name)(debug_adaptor & result, const T& a, const debug_adaptor& b) \ { \ using default_ops::BOOST_JOIN(eval_, name); \ BOOST_JOIN(eval_, name) \ (result.value(), unwrap_debug_type(a), unwrap_debug_type(b)); \ result.update_view(); \ } \ template \ inline void BOOST_JOIN(eval_, name)(debug_adaptor & result, const debug_adaptor& a, const debug_adaptor& b) \ { \ using default_ops::BOOST_JOIN(eval_, name); \ BOOST_JOIN(eval_, name) \ (result.value(), unwrap_debug_type(a), unwrap_debug_type(b)); \ result.update_view(); \ } #define NON_MEMBER_OP4(name, str) \ template \ inline void BOOST_JOIN(eval_, name)(debug_adaptor & result, const T& a, const U& b, const V& c) \ { \ using default_ops::BOOST_JOIN(eval_, name); \ BOOST_JOIN(eval_, name) \ (result.value(), unwrap_debug_type(a), unwrap_debug_type(b), unwrap_debug_type(c)); \ result.update_view(); \ } \ template \ inline void BOOST_JOIN(eval_, name)(debug_adaptor & result, const debug_adaptor& a, const debug_adaptor& b, const T& c) \ { \ using default_ops::BOOST_JOIN(eval_, name); \ BOOST_JOIN(eval_, name) \ (result.value(), unwrap_debug_type(a), unwrap_debug_type(b), unwrap_debug_type(c)); \ result.update_view(); \ } \ template \ inline void BOOST_JOIN(eval_, name)(debug_adaptor & result, const debug_adaptor& a, const T& b, const debug_adaptor& c) \ { \ using default_ops::BOOST_JOIN(eval_, name); \ BOOST_JOIN(eval_, name) \ (result.value(), unwrap_debug_type(a), unwrap_debug_type(b), unwrap_debug_type(c)); \ result.update_view(); \ } \ template \ inline void BOOST_JOIN(eval_, name)(debug_adaptor & result, const T& a, const debug_adaptor& b, const debug_adaptor& c) \ { \ using default_ops::BOOST_JOIN(eval_, name); \ BOOST_JOIN(eval_, name) \ (result.value(), unwrap_debug_type(a), unwrap_debug_type(b), unwrap_debug_type(c)); \ result.update_view(); \ } \ template \ inline void BOOST_JOIN(eval_, name)(debug_adaptor & result, const debug_adaptor& a, const debug_adaptor& b, const debug_adaptor& c) \ { \ using default_ops::BOOST_JOIN(eval_, name); \ BOOST_JOIN(eval_, name) \ (result.value(), unwrap_debug_type(a), unwrap_debug_type(b), unwrap_debug_type(c)); \ result.update_view(); \ } \ template \ inline void BOOST_JOIN(eval_, name)(debug_adaptor & result, const debug_adaptor& a, const T& b, const U& c) \ { \ using default_ops::BOOST_JOIN(eval_, name); \ BOOST_JOIN(eval_, name) \ (result.value(), unwrap_debug_type(a), unwrap_debug_type(b), unwrap_debug_type(c)); \ result.update_view(); \ } NON_MEMBER_OP2(add, "+=") NON_MEMBER_OP2(subtract, "-=") NON_MEMBER_OP2(multiply, "*=") NON_MEMBER_OP2(divide, "/=") template inline void eval_convert_to(R* result, const debug_adaptor& val) { using default_ops::eval_convert_to; eval_convert_to(result, val.value()); } template inline void eval_convert_to(debug_adaptor* result, const debug_adaptor& val) { using default_ops::eval_convert_to; eval_convert_to(&result->value(), val.value()); } template inline void eval_convert_to(debug_adaptor* result, const Backend& val) { using default_ops::eval_convert_to; eval_convert_to(&result->value(), val); } template inline void eval_convert_to(std::complex* result, const debug_adaptor& val) { using default_ops::eval_convert_to; eval_convert_to(result, val.value()); } template inline void eval_convert_to(std::complex* result, const debug_adaptor& val) { using default_ops::eval_convert_to; eval_convert_to(result, val.value()); } template inline void eval_convert_to(std::complex* result, const debug_adaptor& val) { using default_ops::eval_convert_to; eval_convert_to(result, val.value()); } template inline void eval_frexp(debug_adaptor& result, const debug_adaptor& arg, Exp* exp) { eval_frexp(result.value(), arg.value(), exp); result.update_view(); } template inline void eval_ldexp(debug_adaptor& result, const debug_adaptor& arg, Exp exp) { eval_ldexp(result.value(), arg.value(), exp); result.update_view(); } template inline void eval_scalbn(debug_adaptor& result, const debug_adaptor& arg, Exp exp) { using default_ops::eval_scalbn; eval_scalbn(result.value(), arg.value(), exp); result.update_view(); } template inline typename Backend::exponent_type eval_ilogb(const debug_adaptor& arg) { using default_ops::eval_ilogb; return eval_ilogb(arg.value()); } NON_MEMBER_OP2(floor, "floor") NON_MEMBER_OP2(ceil, "ceil") NON_MEMBER_OP2(sqrt, "sqrt") NON_MEMBER_OP2(logb, "logb") template inline int eval_fpclassify(const debug_adaptor& arg) { using default_ops::eval_fpclassify; return eval_fpclassify(arg.value()); } /********************************************************************* * * Optional arithmetic operations come next: * *********************************************************************/ NON_MEMBER_OP3(add, "+") NON_MEMBER_OP3(subtract, "-") NON_MEMBER_OP3(multiply, "*") NON_MEMBER_OP3(divide, "/") NON_MEMBER_OP3(multiply_add, "fused-multiply-add") NON_MEMBER_OP3(multiply_subtract, "fused-multiply-subtract") NON_MEMBER_OP4(multiply_add, "fused-multiply-add") NON_MEMBER_OP4(multiply_subtract, "fused-multiply-subtract") NON_MEMBER_OP1(increment, "increment") NON_MEMBER_OP1(decrement, "decrement") /********************************************************************* * * Optional integer operations come next: * *********************************************************************/ NON_MEMBER_OP2(modulus, "%=") NON_MEMBER_OP3(modulus, "%") NON_MEMBER_OP2(bitwise_or, "|=") NON_MEMBER_OP3(bitwise_or, "|") NON_MEMBER_OP2(bitwise_and, "&=") NON_MEMBER_OP3(bitwise_and, "&") NON_MEMBER_OP2(bitwise_xor, "^=") NON_MEMBER_OP3(bitwise_xor, "^") NON_MEMBER_OP4(qr, "quotient-and-remainder") NON_MEMBER_OP2(complement, "~") template inline void eval_left_shift(debug_adaptor& arg, std::size_t a) { using default_ops::eval_left_shift; eval_left_shift(arg.value(), a); arg.update_view(); } template inline void eval_left_shift(debug_adaptor& arg, const debug_adaptor& a, std::size_t b) { using default_ops::eval_left_shift; eval_left_shift(arg.value(), a.value(), b); arg.update_view(); } template inline void eval_right_shift(debug_adaptor& arg, std::size_t a) { using default_ops::eval_right_shift; eval_right_shift(arg.value(), a); arg.update_view(); } template inline void eval_right_shift(debug_adaptor& arg, const debug_adaptor& a, std::size_t b) { using default_ops::eval_right_shift; eval_right_shift(arg.value(), a.value(), b); arg.update_view(); } template inline T eval_integer_modulus(const debug_adaptor& arg, const T& a) { using default_ops::eval_integer_modulus; return eval_integer_modulus(arg.value(), a); } template inline std::size_t eval_lsb(const debug_adaptor& arg) { using default_ops::eval_lsb; return eval_lsb(arg.value()); } template inline std::size_t eval_msb(const debug_adaptor& arg) { using default_ops::eval_msb; return eval_msb(arg.value()); } template inline bool eval_bit_test(const debug_adaptor& arg, std::size_t a) { using default_ops::eval_bit_test; return eval_bit_test(arg.value(), a); } template inline void eval_bit_set(const debug_adaptor& arg, std::size_t a) { using default_ops::eval_bit_set; eval_bit_set(arg.value(), a); arg.update_view(); } template inline void eval_bit_unset(const debug_adaptor& arg, std::size_t a) { using default_ops::eval_bit_unset; eval_bit_unset(arg.value(), a); arg.update_view(); } template inline void eval_bit_flip(const debug_adaptor& arg, std::size_t a) { using default_ops::eval_bit_flip; eval_bit_flip(arg.value(), a); arg.update_view(); } NON_MEMBER_OP3(gcd, "gcd") NON_MEMBER_OP3(lcm, "lcm") NON_MEMBER_OP4(powm, "powm") /********************************************************************* * * abs/fabs: * *********************************************************************/ NON_MEMBER_OP2(abs, "abs") NON_MEMBER_OP2(fabs, "fabs") /********************************************************************* * * Floating point functions: * *********************************************************************/ NON_MEMBER_OP2(trunc, "trunc") NON_MEMBER_OP2(round, "round") NON_MEMBER_OP2(exp, "exp") NON_MEMBER_OP2(log, "log") NON_MEMBER_OP2(log10, "log10") NON_MEMBER_OP2(sin, "sin") NON_MEMBER_OP2(cos, "cos") NON_MEMBER_OP2(tan, "tan") NON_MEMBER_OP2(asin, "asin") NON_MEMBER_OP2(acos, "acos") NON_MEMBER_OP2(atan, "atan") NON_MEMBER_OP2(sinh, "sinh") NON_MEMBER_OP2(cosh, "cosh") NON_MEMBER_OP2(tanh, "tanh") NON_MEMBER_OP2(asinh, "asinh") NON_MEMBER_OP2(acosh, "acosh") NON_MEMBER_OP2(atanh, "atanh") NON_MEMBER_OP3(fmod, "fmod") NON_MEMBER_OP3(pow, "pow") NON_MEMBER_OP3(atan2, "atan2") NON_MEMBER_OP2(conj, "conj") template int eval_signbit(const debug_adaptor& val) { using default_ops::eval_signbit; return eval_signbit(val.value()); } template std::size_t hash_value(const debug_adaptor& val) { return hash_value(val.value()); } template inline typename std::enable_if::value == number_kind_rational, typename number, ExpressionTemplates>::value_type>::type numerator(const number, ExpressionTemplates>& arg) { number t(arg.backend().value()); return numerator(t).backend(); } template inline typename std::enable_if::value == number_kind_rational, typename number, ExpressionTemplates>::value_type>::type denominator(const number, ExpressionTemplates>& arg) { number t(arg.backend().value()); return denominator(t).backend(); } template inline BOOST_MP_CXX14_CONSTEXPR void eval_real(To& to, const debug_adaptor& from) { using default_ops::eval_real; eval_real(to, from.value()); } template inline BOOST_MP_CXX14_CONSTEXPR void eval_real(debug_adaptor& to, const debug_adaptor& from) { using default_ops::eval_real; eval_real(to.value(), from.value()); to.update_view(); } template inline BOOST_MP_CXX14_CONSTEXPR void eval_real(debug_adaptor& to, const From& from) { using default_ops::eval_real; eval_real(to.value(), from); to.update_view(); } template inline BOOST_MP_CXX14_CONSTEXPR void eval_imag(To& to, const debug_adaptor& from) { using default_ops::eval_imag; eval_imag(to, from.value()); } template inline BOOST_MP_CXX14_CONSTEXPR void eval_imag(debug_adaptor& to, const debug_adaptor& from) { using default_ops::eval_imag; eval_imag(to.value(), from.value()); to.update_view(); } template inline BOOST_MP_CXX14_CONSTEXPR void eval_imag(debug_adaptor& to, const From& from) { using default_ops::eval_imag; eval_imag(to.value(), from); to.update_view(); } template inline BOOST_MP_CXX14_CONSTEXPR void eval_set_real(To& to, const debug_adaptor& from) { using default_ops::eval_set_real; eval_set_real(to, from.value()); } template inline BOOST_MP_CXX14_CONSTEXPR void eval_set_real(debug_adaptor& to, const debug_adaptor& from) { using default_ops::eval_set_real; eval_set_real(to.value(), from.value()); to.update_view(); } template inline BOOST_MP_CXX14_CONSTEXPR void eval_set_real(debug_adaptor& to, const From& from) { using default_ops::eval_set_real; eval_set_real(to.value(), from); to.update_view(); } template inline BOOST_MP_CXX14_CONSTEXPR void eval_set_imag(To& to, const debug_adaptor& from) { using default_ops::eval_set_imag; eval_set_imag(to, from.value()); } template inline BOOST_MP_CXX14_CONSTEXPR void eval_set_imag(debug_adaptor& to, const debug_adaptor& from) { using default_ops::eval_set_imag; eval_set_imag(to.value(), from.value()); to.update_view(); } template inline BOOST_MP_CXX14_CONSTEXPR void eval_set_imag(debug_adaptor& to, const From& from) { using default_ops::eval_set_imag; eval_set_imag(to.value(), from); to.update_view(); } } // namespace backends namespace detail { template struct is_variable_precision > : public is_variable_precision {}; #ifdef BOOST_HAS_INT128 template struct is_convertible_arithmetic > : public is_convertible_arithmetic {}; template struct is_convertible_arithmetic > : public is_convertible_arithmetic {}; #endif #ifdef BOOST_HAS_FLOAT128 template struct is_convertible_arithmetic > : public is_convertible_arithmetic {}; #endif } // namespace detail template struct number_category > : public number_category {}; template struct component_type, ExpressionTemplates>> { // // We'll make the component_type also a debug_adaptor: // using base_component_type = typename component_type>::type; using base_component_backend = typename base_component_type::backend_type; using type = number, ExpressionTemplates>; }; template struct is_interval_number > : public is_interval_number {}; }} // namespace boost::multiprecision namespace std { template class numeric_limits, ExpressionTemplates> > : public std::numeric_limits > { using base_type = std::numeric_limits > ; using number_type = boost::multiprecision::number, ExpressionTemplates>; public: static number_type(min)() noexcept { return (base_type::min)(); } static number_type(max)() noexcept { return (base_type::max)(); } static number_type lowest() noexcept { return -(max)(); } static number_type epsilon() noexcept { return base_type::epsilon(); } static number_type round_error() noexcept { return epsilon() / 2; } static number_type infinity() noexcept { return base_type::infinity(); } static number_type quiet_NaN() noexcept { return base_type::quiet_NaN(); } static number_type signaling_NaN() noexcept { return base_type::signaling_NaN(); } static number_type denorm_min() noexcept { return base_type::denorm_min(); } }; } // namespace std #ifdef BOOST_MP_MATH_AVAILABLE namespace boost { namespace math { namespace policies { template struct precision, ExpressionTemplates>, Policy> : public precision, Policy> {}; #undef NON_MEMBER_OP1 #undef NON_MEMBER_OP2 #undef NON_MEMBER_OP3 #undef NON_MEMBER_OP4 } }} // namespace boost::math::policies #else #undef NON_MEMBER_OP1 #undef NON_MEMBER_OP2 #undef NON_MEMBER_OP3 #undef NON_MEMBER_OP4 #endif // BOOST_MP_MATH_AVAILABLE #endif