123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452453454455456457458459460461462463464465466467468469470471472473474475476477478479480481482483484485486487488489490491492493494495496497498499500501502503504505506507508509510511512513514515516517518519520521522523524525526527528529530531532 |
- ///////////////////////////////////////////////////////////////
- // 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_MATH_DEBUG_ADAPTER_HPP
- #define BOOST_MATH_DEBUG_ADAPTER_HPP
- #include <boost/multiprecision/traits/extract_exponent_type.hpp>
- #include <boost/multiprecision/detail/integer_ops.hpp>
- namespace boost {
- namespace multiprecision {
- namespace backends {
- template <class Backend>
- 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<Backend, number_category<Backend>::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<std::ios_base::fmtflags>(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 <class T>
- debug_adaptor(const T& i, const typename std::enable_if<std::is_convertible<T, Backend>::value>::type* = 0)
- : m_value(i)
- {
- update_view();
- }
- template <class T>
- debug_adaptor(const T& i, const T& j)
- : m_value(i, j)
- {
- update_view();
- }
- template <class T>
- typename std::enable_if<boost::multiprecision::detail::is_arithmetic<T>::value || std::is_convertible<T, Backend>::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 <class T>
- int compare(const T& i) const
- {
- return m_value.compare(i);
- }
- Backend& value()
- {
- return m_value;
- }
- const Backend& value() const
- {
- return m_value;
- }
- template <class Archive>
- 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();
- }
- static unsigned default_precision() noexcept
- {
- return Backend::default_precision();
- }
- static void default_precision(unsigned v) noexcept
- {
- Backend::default_precision(v);
- }
- unsigned precision() const noexcept
- {
- return value().precision();
- }
- void precision(unsigned digits10) noexcept
- {
- value().precision(digits10);
- }
- };
- template <class Backend>
- inline Backend const& unwrap_debug_type(debug_adaptor<Backend> const& val)
- {
- return val.value();
- }
- template <class T>
- inline const T& unwrap_debug_type(const T& val)
- {
- return val;
- }
- #define NON_MEMBER_OP1(name, str) \
- template <class Backend> \
- inline void BOOST_JOIN(eval_, name)(debug_adaptor<Backend> & result) \
- { \
- using default_ops::BOOST_JOIN(eval_, name); \
- BOOST_JOIN(eval_, name) \
- (result.value()); \
- result.update_view(); \
- }
- #define NON_MEMBER_OP2(name, str) \
- template <class Backend, class T> \
- inline void BOOST_JOIN(eval_, name)(debug_adaptor<Backend> & 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 <class Backend> \
- inline void BOOST_JOIN(eval_, name)(debug_adaptor<Backend> & result, const debug_adaptor<Backend>& 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 <class Backend, class T, class U> \
- inline void BOOST_JOIN(eval_, name)(debug_adaptor<Backend> & 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 <class Backend, class T> \
- inline void BOOST_JOIN(eval_, name)(debug_adaptor<Backend> & result, const debug_adaptor<Backend>& 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 <class Backend, class T> \
- inline void BOOST_JOIN(eval_, name)(debug_adaptor<Backend> & result, const T& a, const debug_adaptor<Backend>& 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 <class Backend> \
- inline void BOOST_JOIN(eval_, name)(debug_adaptor<Backend> & result, const debug_adaptor<Backend>& a, const debug_adaptor<Backend>& 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 <class Backend, class T, class U, class V> \
- inline void BOOST_JOIN(eval_, name)(debug_adaptor<Backend> & 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 <class Backend, class T> \
- inline void BOOST_JOIN(eval_, name)(debug_adaptor<Backend> & result, const debug_adaptor<Backend>& a, const debug_adaptor<Backend>& 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 <class Backend, class T> \
- inline void BOOST_JOIN(eval_, name)(debug_adaptor<Backend> & result, const debug_adaptor<Backend>& a, const T& b, const debug_adaptor<Backend>& 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 <class Backend, class T> \
- inline void BOOST_JOIN(eval_, name)(debug_adaptor<Backend> & result, const T& a, const debug_adaptor<Backend>& b, const debug_adaptor<Backend>& 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 <class Backend> \
- inline void BOOST_JOIN(eval_, name)(debug_adaptor<Backend> & result, const debug_adaptor<Backend>& a, const debug_adaptor<Backend>& b, const debug_adaptor<Backend>& 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 <class Backend, class T, class U> \
- inline void BOOST_JOIN(eval_, name)(debug_adaptor<Backend> & result, const debug_adaptor<Backend>& 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 <class Backend, class R>
- inline void eval_convert_to(R* result, const debug_adaptor<Backend>& val)
- {
- using default_ops::eval_convert_to;
- eval_convert_to(result, val.value());
- }
- template <class Backend, class Exp>
- inline void eval_frexp(debug_adaptor<Backend>& result, const debug_adaptor<Backend>& arg, Exp* exp)
- {
- eval_frexp(result.value(), arg.value(), exp);
- result.update_view();
- }
- template <class Backend, class Exp>
- inline void eval_ldexp(debug_adaptor<Backend>& result, const debug_adaptor<Backend>& arg, Exp exp)
- {
- eval_ldexp(result.value(), arg.value(), exp);
- result.update_view();
- }
- template <class Backend, class Exp>
- inline void eval_scalbn(debug_adaptor<Backend>& result, const debug_adaptor<Backend>& arg, Exp exp)
- {
- using default_ops::eval_scalbn;
- eval_scalbn(result.value(), arg.value(), exp);
- result.update_view();
- }
- template <class Backend>
- inline typename Backend::exponent_type eval_ilogb(const debug_adaptor<Backend>& 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 <class Backend>
- inline int eval_fpclassify(const debug_adaptor<Backend>& 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 <class Backend>
- inline void eval_left_shift(debug_adaptor<Backend>& arg, std::size_t a)
- {
- using default_ops::eval_left_shift;
- eval_left_shift(arg.value(), a);
- arg.update_view();
- }
- template <class Backend>
- inline void eval_left_shift(debug_adaptor<Backend>& arg, const debug_adaptor<Backend>& a, std::size_t b)
- {
- using default_ops::eval_left_shift;
- eval_left_shift(arg.value(), a.value(), b);
- arg.update_view();
- }
- template <class Backend>
- inline void eval_right_shift(debug_adaptor<Backend>& arg, std::size_t a)
- {
- using default_ops::eval_right_shift;
- eval_right_shift(arg.value(), a);
- arg.update_view();
- }
- template <class Backend>
- inline void eval_right_shift(debug_adaptor<Backend>& arg, const debug_adaptor<Backend>& a, std::size_t b)
- {
- using default_ops::eval_right_shift;
- eval_right_shift(arg.value(), a.value(), b);
- arg.update_view();
- }
- template <class Backend, class T>
- inline unsigned eval_integer_modulus(const debug_adaptor<Backend>& arg, const T& a)
- {
- using default_ops::eval_integer_modulus;
- return eval_integer_modulus(arg.value(), a);
- }
- template <class Backend>
- inline unsigned eval_lsb(const debug_adaptor<Backend>& arg)
- {
- using default_ops::eval_lsb;
- return eval_lsb(arg.value());
- }
- template <class Backend>
- inline unsigned eval_msb(const debug_adaptor<Backend>& arg)
- {
- using default_ops::eval_msb;
- return eval_msb(arg.value());
- }
- template <class Backend>
- inline bool eval_bit_test(const debug_adaptor<Backend>& arg, unsigned a)
- {
- using default_ops::eval_bit_test;
- return eval_bit_test(arg.value(), a);
- }
- template <class Backend>
- inline void eval_bit_set(const debug_adaptor<Backend>& arg, unsigned a)
- {
- using default_ops::eval_bit_set;
- eval_bit_set(arg.value(), a);
- arg.update_view();
- }
- template <class Backend>
- inline void eval_bit_unset(const debug_adaptor<Backend>& arg, unsigned a)
- {
- using default_ops::eval_bit_unset;
- eval_bit_unset(arg.value(), a);
- arg.update_view();
- }
- template <class Backend>
- inline void eval_bit_flip(const debug_adaptor<Backend>& arg, unsigned 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_OP3(fmod, "fmod")
- NON_MEMBER_OP3(pow, "pow")
- NON_MEMBER_OP3(atan2, "atan2")
- template <class Backend>
- int eval_signbit(const debug_adaptor<Backend>& val)
- {
- return eval_signbit(val.value());
- }
- template <class Backend>
- std::size_t hash_value(const debug_adaptor<Backend>& val)
- {
- return hash_value(val.value());
- }
- } // namespace backends
- using backends::debug_adaptor;
- template <class Backend>
- struct number_category<backends::debug_adaptor<Backend> > : public number_category<Backend>
- {};
- }} // namespace boost::multiprecision
- namespace std {
- template <class Backend, boost::multiprecision::expression_template_option ExpressionTemplates>
- class numeric_limits<boost::multiprecision::number<boost::multiprecision::backends::debug_adaptor<Backend>, ExpressionTemplates> >
- : public std::numeric_limits<boost::multiprecision::number<Backend, ExpressionTemplates> >
- {
- using base_type = std::numeric_limits<boost::multiprecision::number<Backend, ExpressionTemplates> > ;
- using number_type = boost::multiprecision::number<boost::multiprecision::backends::debug_adaptor<Backend>, 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
- namespace boost {
- namespace math {
- namespace policies {
- template <class Backend, boost::multiprecision::expression_template_option ExpressionTemplates, class Policy>
- struct precision<boost::multiprecision::number<boost::multiprecision::debug_adaptor<Backend>, ExpressionTemplates>, Policy>
- : public precision<boost::multiprecision::number<Backend, ExpressionTemplates>, Policy>
- {};
- #undef NON_MEMBER_OP1
- #undef NON_MEMBER_OP2
- #undef NON_MEMBER_OP3
- #undef NON_MEMBER_OP4
- }
- }} // namespace boost::math::policies
- #endif
|