#include <algorithm>
#ifndef BOOST_MULTIPRECISION_MPFI_DEFAULT_PRECISION
-# define BOOST_MULTIPRECISION_MPFI_DEFAULT_PRECISION 20
+#define BOOST_MULTIPRECISION_MPFI_DEFAULT_PRECISION 20
#endif
-namespace boost{
-namespace multiprecision{
-namespace backends{
+namespace boost {
+namespace multiprecision {
+namespace backends {
template <unsigned digits10>
struct mpfi_float_backend;
} // namespace backends
template <unsigned digits10>
-struct number_category<backends::mpfi_float_backend<digits10> > : public mpl::int_<number_kind_floating_point>{};
+struct number_category<backends::mpfi_float_backend<digits10> > : public mpl::int_<number_kind_floating_point>
+{};
struct interval_error : public std::runtime_error
{
interval_error(const std::string& s) : std::runtime_error(s) {}
};
-namespace backends{
+namespace backends {
-namespace detail{
+namespace detail {
inline int mpfi_sgn(mpfi_srcptr p)
{
- if(mpfi_is_zero(p))
+ if (mpfi_is_zero(p))
return 0;
- if(mpfi_is_strictly_pos(p))
+ if (mpfi_is_strictly_pos(p))
return 1;
- if(mpfi_is_strictly_neg(p))
+ if (mpfi_is_strictly_neg(p))
return -1;
BOOST_THROW_EXCEPTION(interval_error("Sign of interval is ambiguous."));
}
struct mpfi_float_imp
{
#ifdef BOOST_HAS_LONG_LONG
- typedef mpl::list<long, boost::long_long_type> signed_types;
- typedef mpl::list<unsigned long, boost::ulong_long_type> unsigned_types;
+ typedef mpl::list<long, boost::long_long_type> signed_types;
+ typedef mpl::list<unsigned long, boost::ulong_long_type> unsigned_types;
#else
- typedef mpl::list<long> signed_types;
- typedef mpl::list<unsigned long> unsigned_types;
+ typedef mpl::list<long> signed_types;
+ typedef mpl::list<unsigned long> unsigned_types;
#endif
- typedef mpl::list<double, long double> float_types;
- typedef long exponent_type;
+ typedef mpl::list<double, long double> float_types;
+ typedef long exponent_type;
mpfi_float_imp()
{
mpfi_float_imp(const mpfi_float_imp& o)
{
- mpfi_init2(m_data, multiprecision::detail::digits10_2_2(digits10 ? digits10 : get_default_precision()));
- if(o.m_data[0].left._mpfr_d)
+ mpfi_init2(m_data, mpfi_get_prec(o.data()));
+ if (o.m_data[0].left._mpfr_d)
mpfi_set(m_data, o.m_data);
}
#ifndef BOOST_NO_CXX11_RVALUE_REFERENCES
mpfi_float_imp(mpfi_float_imp&& o) BOOST_NOEXCEPT
{
- m_data[0] = o.m_data[0];
+ m_data[0] = o.m_data[0];
o.m_data[0].left._mpfr_d = 0;
}
#endif
- mpfi_float_imp& operator = (const mpfi_float_imp& o)
+ mpfi_float_imp& operator=(const mpfi_float_imp& o)
{
- if(m_data[0].left._mpfr_d == 0)
- mpfi_init2(m_data, multiprecision::detail::digits10_2_2(digits10 ? digits10 : get_default_precision()));
- if(o.m_data[0].left._mpfr_d)
- mpfi_set(m_data, o.m_data);
+ if (m_data[0].left._mpfr_d == 0)
+ mpfi_init2(m_data, mpfi_get_prec(o.data()));
+ if (mpfi_get_prec(o.data()) != mpfi_get_prec(data()))
+ {
+ mpfi_float_imp t(mpfi_get_prec(o.data()));
+ t = o;
+ t.swap(*this);
+ }
+ else
+ {
+ if (o.m_data[0].left._mpfr_d)
+ mpfi_set(m_data, o.m_data);
+ }
return *this;
}
#ifndef BOOST_NO_CXX11_RVALUE_REFERENCES
- mpfi_float_imp& operator = (mpfi_float_imp&& o) BOOST_NOEXCEPT
+ mpfi_float_imp& operator=(mpfi_float_imp&& o) BOOST_NOEXCEPT
{
mpfi_swap(m_data, o.m_data);
return *this;
#endif
#ifdef BOOST_HAS_LONG_LONG
#ifdef _MPFR_H_HAVE_INTMAX_T
- mpfi_float_imp& operator = (boost::ulong_long_type i)
+ mpfi_float_imp& operator=(boost::ulong_long_type i)
{
- if(m_data[0].left._mpfr_d == 0)
+ if (m_data[0].left._mpfr_d == 0)
mpfi_init2(m_data, multiprecision::detail::digits10_2_2(digits10 ? digits10 : get_default_precision()));
mpfr_set_uj(left_data(), i, GMP_RNDD);
mpfr_set_uj(right_data(), i, GMP_RNDU);
return *this;
}
- mpfi_float_imp& operator = (boost::long_long_type i)
+ mpfi_float_imp& operator=(boost::long_long_type i)
{
- if(m_data[0].left._mpfr_d == 0)
+ if (m_data[0].left._mpfr_d == 0)
mpfi_init2(m_data, multiprecision::detail::digits10_2_2(digits10 ? digits10 : get_default_precision()));
mpfr_set_sj(left_data(), i, GMP_RNDD);
mpfr_set_sj(right_data(), i, GMP_RNDU);
return *this;
}
#else
- mpfi_float_imp& operator = (boost::ulong_long_type i)
+ mpfi_float_imp& operator=(boost::ulong_long_type i)
{
- if(m_data[0].left._mpfr_d == 0)
+ if (m_data[0].left._mpfr_d == 0)
mpfi_init2(m_data, multiprecision::detail::digits10_2_2(digits10 ? digits10 : get_default_precision()));
- boost::ulong_long_type mask = ((((1uLL << (std::numeric_limits<unsigned long>::digits - 1)) - 1) << 1) | 1u);
- unsigned shift = 0;
- mpfi_t t;
+ boost::ulong_long_type mask = ((((1uLL << (std::numeric_limits<unsigned long>::digits - 1)) - 1) << 1) | 1u);
+ unsigned shift = 0;
+ mpfi_t t;
mpfi_init2(t, (std::max)(static_cast<unsigned long>(std::numeric_limits<boost::ulong_long_type>::digits), static_cast<unsigned long>(multiprecision::detail::digits10_2_2(digits10))));
mpfi_set_ui(m_data, 0);
- while(i)
+ while (i)
{
mpfi_set_ui(t, static_cast<unsigned long>(i & mask));
- if(shift)
+ if (shift)
mpfi_mul_2exp(t, t, shift);
mpfi_add(m_data, m_data, t);
shift += std::numeric_limits<unsigned long>::digits;
mpfi_clear(t);
return *this;
}
- mpfi_float_imp& operator = (boost::long_long_type i)
+ mpfi_float_imp& operator=(boost::long_long_type i)
{
- if(m_data[0].left._mpfr_d == 0)
+ if (m_data[0].left._mpfr_d == 0)
mpfi_init2(m_data, multiprecision::detail::digits10_2_2(digits10 ? digits10 : get_default_precision()));
bool neg = i < 0;
- *this = boost::multiprecision::detail::unsigned_abs(i);
- if(neg)
+ *this = boost::multiprecision::detail::unsigned_abs(i);
+ if (neg)
mpfi_neg(m_data, m_data);
return *this;
}
#endif
#endif
- mpfi_float_imp& operator = (unsigned long i)
+ mpfi_float_imp& operator=(unsigned long i)
{
- if(m_data[0].left._mpfr_d == 0)
+ if (m_data[0].left._mpfr_d == 0)
mpfi_init2(m_data, multiprecision::detail::digits10_2_2(digits10 ? digits10 : get_default_precision()));
mpfi_set_ui(m_data, i);
return *this;
}
- mpfi_float_imp& operator = (long i)
+ mpfi_float_imp& operator=(long i)
{
- if(m_data[0].left._mpfr_d == 0)
+ if (m_data[0].left._mpfr_d == 0)
mpfi_init2(m_data, multiprecision::detail::digits10_2_2(digits10 ? digits10 : get_default_precision()));
mpfi_set_si(m_data, i);
return *this;
}
- mpfi_float_imp& operator = (double d)
+ mpfi_float_imp& operator=(double d)
{
- if(m_data[0].left._mpfr_d == 0)
+ if (m_data[0].left._mpfr_d == 0)
mpfi_init2(m_data, multiprecision::detail::digits10_2_2(digits10 ? digits10 : get_default_precision()));
mpfi_set_d(m_data, d);
return *this;
}
- mpfi_float_imp& operator = (long double a)
+ mpfi_float_imp& operator=(long double a)
{
- if(m_data[0].left._mpfr_d == 0)
+ if (m_data[0].left._mpfr_d == 0)
mpfi_init2(m_data, multiprecision::detail::digits10_2_2(digits10 ? digits10 : get_default_precision()));
mpfr_set_ld(left_data(), a, GMP_RNDD);
mpfr_set_ld(right_data(), a, GMP_RNDU);
return *this;
}
- mpfi_float_imp& operator = (const char* s)
+ mpfi_float_imp& operator=(const char* s)
{
using default_ops::eval_fpclassify;
-
- if(m_data[0].left._mpfr_d == 0)
+
+ if (m_data[0].left._mpfr_d == 0)
mpfi_init2(m_data, multiprecision::detail::digits10_2_2(digits10 ? digits10 : get_default_precision()));
- if(s && (*s == '{'))
+ if (s && (*s == '{'))
{
mpfr_float_backend<digits10> a, b;
- std::string part;
- const char* p = ++s;
- while(*p && (*p != ',') && (*p != '}'))
+ std::string part;
+ const char* p = ++s;
+ while (*p && (*p != ',') && (*p != '}'))
++p;
part.assign(s + 1, p);
a = part.c_str();
s = p;
- if(*p && (*p != '}'))
+ if (*p && (*p != '}'))
{
++p;
- while(*p && (*p != ',') && (*p != '}'))
+ while (*p && (*p != ',') && (*p != '}'))
++p;
part.assign(s + 1, p);
}
part.erase();
b = part.c_str();
- if(eval_fpclassify(a) == (int)FP_NAN)
+ if (eval_fpclassify(a) == (int)FP_NAN)
{
mpfi_set_fr(this->data(), a.data());
}
- else if(eval_fpclassify(b) == (int)FP_NAN)
+ else if (eval_fpclassify(b) == (int)FP_NAN)
{
mpfi_set_fr(this->data(), b.data());
}
else
{
- if(a.compare(b) > 0)
+ if (a.compare(b) > 0)
{
BOOST_THROW_EXCEPTION(std::runtime_error("Attempt to create interval with invalid range (start is greater than end)."));
}
mpfi_interv_fr(m_data, a.data(), b.data());
}
}
- else if(mpfi_set_str(m_data, s, 10) != 0)
+ else if (mpfi_set_str(m_data, s, 10) != 0)
{
BOOST_THROW_EXCEPTION(std::runtime_error(std::string("Unable to parse string \"") + s + std::string("\"as a valid floating point number.")));
}
{
mpfi_swap(m_data, o.m_data);
}
- std::string str(std::streamsize digits, std::ios_base::fmtflags f)const
+ std::string str(std::streamsize digits, std::ios_base::fmtflags f) const
{
BOOST_ASSERT(m_data[0].left._mpfr_d);
mpfi_get_left(a.data(), m_data);
mpfi_get_right(b.data(), m_data);
- if(a.compare(b) == 0)
+ if (a.compare(b) == 0)
return a.str(digits, f);
return "{" + a.str(digits, f) + "," + b.str(digits, f) + "}";
}
~mpfi_float_imp() BOOST_NOEXCEPT
{
- if(m_data[0].left._mpfr_d)
+ if (m_data[0].left._mpfr_d)
mpfi_clear(m_data);
}
void negate() BOOST_NOEXCEPT
BOOST_ASSERT(m_data[0].left._mpfr_d);
mpfi_neg(m_data, m_data);
}
- int compare(const mpfi_float_imp& o)const BOOST_NOEXCEPT
+ int compare(const mpfi_float_imp& o) const BOOST_NOEXCEPT
{
BOOST_ASSERT(m_data[0].left._mpfr_d && o.m_data[0].left._mpfr_d);
- if(mpfr_cmp(right_data(), o.left_data()) < 0)
+ if (mpfr_cmp(right_data(), o.left_data()) < 0)
return -1;
- if(mpfr_cmp(left_data(), o.right_data()) > 0)
+ if (mpfr_cmp(left_data(), o.right_data()) > 0)
return 1;
- if((mpfr_cmp(left_data(), o.left_data()) == 0) && (mpfr_cmp(right_data(), o.right_data()) == 0))
+ if ((mpfr_cmp(left_data(), o.left_data()) == 0) && (mpfr_cmp(right_data(), o.right_data()) == 0))
return 0;
BOOST_THROW_EXCEPTION(interval_error("Ambiguous comparison between two values."));
return 0;
}
template <class V>
- int compare(V v)const BOOST_NOEXCEPT
+ int compare(V v) const BOOST_NOEXCEPT
{
mpfi_float_imp d;
d = v;
BOOST_ASSERT(m_data[0].left._mpfr_d);
return m_data;
}
- const mpfi_t& data()const BOOST_NOEXCEPT
+ const mpfi_t& data() const BOOST_NOEXCEPT
{
BOOST_ASSERT(m_data[0].left._mpfr_d);
return m_data;
BOOST_ASSERT(m_data[0].left._mpfr_d);
return &(m_data[0].left);
}
- mpfr_srcptr left_data()const BOOST_NOEXCEPT
+ mpfr_srcptr left_data() const BOOST_NOEXCEPT
{
BOOST_ASSERT(m_data[0].left._mpfr_d);
return &(m_data[0].left);
BOOST_ASSERT(m_data[0].left._mpfr_d);
return &(m_data[0].right);
}
- mpfr_srcptr right_data()const BOOST_NOEXCEPT
+ mpfr_srcptr right_data() const BOOST_NOEXCEPT
{
BOOST_ASSERT(m_data[0].left._mpfr_d);
return &(m_data[0].right);
}
-protected:
- mpfi_t m_data;
+
+ protected:
+ mpfi_t m_data;
static unsigned& get_default_precision() BOOST_NOEXCEPT
{
static unsigned val = BOOST_MULTIPRECISION_MPFI_DEFAULT_PRECISION;
mpfi_float_backend() : detail::mpfi_float_imp<digits10>() {}
mpfi_float_backend(const mpfi_float_backend& o) : detail::mpfi_float_imp<digits10>(o) {}
#ifndef BOOST_NO_CXX11_RVALUE_REFERENCES
- mpfi_float_backend(mpfi_float_backend&& o) : detail::mpfi_float_imp<digits10>(static_cast<detail::mpfi_float_imp<digits10>&&>(o)) {}
+ mpfi_float_backend(mpfi_float_backend&& o) : detail::mpfi_float_imp<digits10>(static_cast<detail::mpfi_float_imp<digits10>&&>(o))
+ {}
#endif
template <unsigned D>
mpfi_float_backend(const mpfi_float_backend<D>& val, typename enable_if_c<D <= digits10>::type* = 0)
}
};
+template <unsigned D1, unsigned D2, mpfr_allocation_type AllocationType>
+void assign_components(mpfi_float_backend<D1>& result, const mpfr_float_backend<D2, AllocationType>& a, const mpfr_float_backend<D2, AllocationType>& b);
+
+template <unsigned Digits10, class V>
+typename enable_if_c<boost::is_constructible<number<mpfr_float_backend<Digits10, allocate_dynamic>, et_on>, V>::value || boost::is_convertible<V, const char*>::value>::type
+assign_components(mpfi_float_backend<Digits10>& result, const V& a, const V& b);
+
template <>
struct mpfi_float_backend<0> : public detail::mpfi_float_imp<0>
{
mpfi_float_backend() : detail::mpfi_float_imp<0>() {}
mpfi_float_backend(const mpfi_t val)
- : detail::mpfi_float_imp<0>(mpfi_get_prec(val))
+ : detail::mpfi_float_imp<0>(mpfi_get_prec(val))
{
mpfi_set(this->m_data, val);
}
mpfi_float_backend(const mpfi_float_backend& o) : detail::mpfi_float_imp<0>(o) {}
#ifndef BOOST_NO_CXX11_RVALUE_REFERENCES
- mpfi_float_backend(mpfi_float_backend&& o) BOOST_NOEXCEPT : detail::mpfi_float_imp<0>(static_cast<detail::mpfi_float_imp<0>&&>(o)) {}
+ mpfi_float_backend(mpfi_float_backend&& o) BOOST_NOEXCEPT : detail::mpfi_float_imp<0>(static_cast<detail::mpfi_float_imp<0>&&>(o))
+ {}
#endif
mpfi_float_backend(const mpfi_float_backend& o, unsigned digits10)
- : detail::mpfi_float_imp<0>(digits10)
+ : detail::mpfi_float_imp<0>(multiprecision::detail::digits10_2_2(digits10))
{
- *this = o;
+ mpfi_set(this->m_data, o.data());
}
+ template <class V>
+ mpfi_float_backend(const V& a, const V& b, unsigned digits10)
+ : detail::mpfi_float_imp<0>(multiprecision::detail::digits10_2_2(digits10))
+ {
+ assign_components(*this, a, b);
+ }
+
template <unsigned D>
mpfi_float_backend(const mpfi_float_backend<D>& val)
- : detail::mpfi_float_imp<0>(mpfi_get_prec(val.data()))
+ : detail::mpfi_float_imp<0>(mpfi_get_prec(val.data()))
{
mpfi_set(this->m_data, val.data());
}
#ifndef BOOST_NO_CXX11_RVALUE_REFERENCES
mpfi_float_backend& operator=(mpfi_float_backend&& o) BOOST_NOEXCEPT
{
- *static_cast<detail::mpfi_float_imp<0>*>(this) = static_cast<detail::mpfi_float_imp<0> &&>(o);
+ *static_cast<detail::mpfi_float_imp<0>*>(this) = static_cast<detail::mpfi_float_imp<0>&&>(o);
return *this;
}
#endif
{
get_default_precision() = v;
}
- unsigned precision()const BOOST_NOEXCEPT
+ unsigned precision() const BOOST_NOEXCEPT
{
return multiprecision::detail::digits2_2_10(mpfi_get_prec(this->m_data));
}
void precision(unsigned digits10) BOOST_NOEXCEPT
{
- mpfi_set_prec(this->m_data, multiprecision::detail::digits2_2_10((digits10)));
+ mpfi_float_backend t(*this, digits10);
+ this->swap(t);
}
};
template <unsigned D1, unsigned D2>
inline void eval_multiply(mpfi_float_backend<D1>& result, const mpfi_float_backend<D2>& o)
{
- if((void*)&result == (void*)&o)
+ if ((void*)&result == (void*)&o)
mpfi_sqr(result.data(), o.data());
else
mpfi_mul(result.data(), result.data(), o.data());
template <unsigned digits10>
inline void eval_add(mpfi_float_backend<digits10>& result, long i)
{
- if(i > 0)
+ if (i > 0)
mpfi_add_ui(result.data(), result.data(), i);
else
mpfi_sub_ui(result.data(), result.data(), boost::multiprecision::detail::unsigned_abs(i));
template <unsigned digits10>
inline void eval_subtract(mpfi_float_backend<digits10>& result, long i)
{
- if(i > 0)
+ if (i > 0)
mpfi_sub_ui(result.data(), result.data(), i);
else
mpfi_add_ui(result.data(), result.data(), boost::multiprecision::detail::unsigned_abs(i));
inline void eval_multiply(mpfi_float_backend<digits10>& result, long i)
{
mpfi_mul_ui(result.data(), result.data(), boost::multiprecision::detail::unsigned_abs(i));
- if(i < 0)
+ if (i < 0)
mpfi_neg(result.data(), result.data());
}
template <unsigned digits10>
inline void eval_divide(mpfi_float_backend<digits10>& result, long i)
{
mpfi_div_ui(result.data(), result.data(), boost::multiprecision::detail::unsigned_abs(i));
- if(i < 0)
+ if (i < 0)
mpfi_neg(result.data(), result.data());
}
//
template <unsigned D1, unsigned D2>
inline void eval_add(mpfi_float_backend<D1>& a, const mpfi_float_backend<D2>& x, long y)
{
- if(y < 0)
+ if (y < 0)
mpfi_sub_ui(a.data(), x.data(), boost::multiprecision::detail::unsigned_abs(y));
else
mpfi_add_ui(a.data(), x.data(), y);
template <unsigned D1, unsigned D2>
inline void eval_add(mpfi_float_backend<D1>& a, long x, const mpfi_float_backend<D2>& y)
{
- if(x < 0)
+ if (x < 0)
{
mpfi_ui_sub(a.data(), boost::multiprecision::detail::unsigned_abs(x), y.data());
mpfi_neg(a.data(), a.data());
template <unsigned D1, unsigned D2>
inline void eval_subtract(mpfi_float_backend<D1>& a, const mpfi_float_backend<D2>& x, long y)
{
- if(y < 0)
+ if (y < 0)
mpfi_add_ui(a.data(), x.data(), boost::multiprecision::detail::unsigned_abs(y));
else
mpfi_sub_ui(a.data(), x.data(), y);
template <unsigned D1, unsigned D2>
inline void eval_subtract(mpfi_float_backend<D1>& a, long x, const mpfi_float_backend<D2>& y)
{
- if(x < 0)
+ if (x < 0)
{
mpfi_add_ui(a.data(), y.data(), boost::multiprecision::detail::unsigned_abs(x));
mpfi_neg(a.data(), a.data());
template <unsigned D1, unsigned D2, unsigned D3>
inline void eval_multiply(mpfi_float_backend<D1>& a, const mpfi_float_backend<D2>& x, const mpfi_float_backend<D3>& y)
{
- if((void*)&x == (void*)&y)
+ if ((void*)&x == (void*)&y)
mpfi_sqr(a.data(), x.data());
else
mpfi_mul(a.data(), x.data(), y.data());
template <unsigned D1, unsigned D2>
inline void eval_multiply(mpfi_float_backend<D1>& a, const mpfi_float_backend<D2>& x, long y)
{
- if(y < 0)
+ if (y < 0)
{
mpfi_mul_ui(a.data(), x.data(), boost::multiprecision::detail::unsigned_abs(y));
a.negate();
template <unsigned D1, unsigned D2>
inline void eval_multiply(mpfi_float_backend<D1>& a, long x, const mpfi_float_backend<D2>& y)
{
- if(x < 0)
+ if (x < 0)
{
mpfi_mul_ui(a.data(), y.data(), boost::multiprecision::detail::unsigned_abs(x));
mpfi_neg(a.data(), a.data());
template <unsigned D1, unsigned D2>
inline void eval_divide(mpfi_float_backend<D1>& a, const mpfi_float_backend<D2>& x, long y)
{
- if(y < 0)
+ if (y < 0)
{
mpfi_div_ui(a.data(), x.data(), boost::multiprecision::detail::unsigned_abs(y));
a.negate();
template <unsigned D1, unsigned D2>
inline void eval_divide(mpfi_float_backend<D1>& a, long x, const mpfi_float_backend<D2>& y)
{
- if(x < 0)
+ if (x < 0)
{
mpfi_ui_div(a.data(), boost::multiprecision::detail::unsigned_abs(x), y.data());
mpfi_neg(a.data(), a.data());
inline void assign_components(mpfi_float_backend<D1>& result, const mpfr_float_backend<D2, AllocationType>& a, const mpfr_float_backend<D2, AllocationType>& b)
{
using default_ops::eval_fpclassify;
- if(eval_fpclassify(a) == (int)FP_NAN)
+ if (eval_fpclassify(a) == (int)FP_NAN)
{
mpfi_set_fr(result.data(), a.data());
}
- else if(eval_fpclassify(b) == (int)FP_NAN)
+ else if (eval_fpclassify(b) == (int)FP_NAN)
{
mpfi_set_fr(result.data(), b.data());
}
else
{
- if(a.compare(b) > 0)
+ if (a.compare(b) > 0)
{
BOOST_THROW_EXCEPTION(std::runtime_error("Attempt to create interval with invalid range (start is greater than end)."));
}
}
template <unsigned Digits10, class V>
-inline typename enable_if_c<is_convertible<V, number<mpfr_float_backend<Digits10, allocate_dynamic>, et_on> >::value >::type
- assign_components(mpfi_float_backend<Digits10>& result, const V& a, const V& b)
+inline typename enable_if_c<boost::is_constructible<number<mpfr_float_backend<Digits10, allocate_dynamic>, et_on>, V>::value || boost::is_convertible<V, const char*>::value>::type
+assign_components(mpfi_float_backend<Digits10>& result, const V& a, const V& b)
{
number<mpfr_float_backend<Digits10, allocate_dynamic>, et_on> x(a), y(b);
assign_components(result, x.backend(), y.backend());
mpfr_set(b.data(), val.right_data(), GMP_RNDN);
eval_ceil(a, a);
eval_ceil(b, b);
- if(a.compare(b) != 0)
+ if (a.compare(b) != 0)
{
BOOST_THROW_EXCEPTION(interval_error("Attempt to take the ceil of a value that straddles an integer boundary."));
}
mpfr_set(b.data(), val.right_data(), GMP_RNDN);
eval_floor(a, a);
eval_floor(b, b);
- if(a.compare(b) != 0)
+ if (a.compare(b) != 0)
{
BOOST_THROW_EXCEPTION(interval_error("Attempt to take the floor of a value that straddles an integer boundary."));
}
template <unsigned Digits10>
inline void eval_ldexp(mpfi_float_backend<Digits10>& result, const mpfi_float_backend<Digits10>& val, long e)
{
- if(e > 0)
+ if (e > 0)
mpfi_mul_2exp(result.data(), val.data(), e);
- else if(e < 0)
+ else if (e < 0)
mpfi_div_2exp(result.data(), val.data(), -e);
else
result = val;
typedef typename boost::multiprecision::detail::canonical<unsigned, mpfi_float_backend<Digits10> >::type ui_type;
using default_ops::eval_get_sign;
int s = eval_get_sign(b);
- if(s == 0)
+ if (s == 0)
{
- if(eval_get_sign(e) == 0)
+ if (eval_get_sign(e) == 0)
{
result = ui_type(1);
}
}
return;
}
- if(s < 0)
+ if (s < 0)
{
- if(eval_get_sign(e) < 0)
+ if (eval_get_sign(e) < 0)
{
mpfi_float_backend<Digits10> t1, t2;
t1 = e;
#endif
using default_ops::eval_convert_to;
eval_convert_to(&an, e);
- if(e.compare(an) == 0)
+ if (e.compare(an) == 0)
{
mpfi_float_backend<Digits10> pb(b);
pb.negate();
eval_pow(result, pb, e);
- if(an & 1u)
+ if (an & 1u)
result.negate();
return;
}
#ifndef BOOST_NO_EXCEPTIONS
}
- catch(const std::exception&)
+ catch (const std::exception&)
{
// conversion failed, just fall through, value is not an integer.
}
inline std::size_t hash_value(const mpfi_float_backend<Digits10>& val)
{
std::size_t result = 0;
- std::size_t len = val.left_data()[0]._mpfr_prec / mp_bits_per_limb;
- if(val.left_data()[0]._mpfr_prec % mp_bits_per_limb)
+ std::size_t len = val.left_data()[0]._mpfr_prec / mp_bits_per_limb;
+ if (val.left_data()[0]._mpfr_prec % mp_bits_per_limb)
++len;
- for(std::size_t i = 0; i < len; ++i)
+ for (std::size_t i = 0; i < len; ++i)
boost::hash_combine(result, val.left_data()[0]._mpfr_d[i]);
boost::hash_combine(result, val.left_data()[0]._mpfr_exp);
boost::hash_combine(result, val.left_data()[0]._mpfr_sign);
len = val.right_data()[0]._mpfr_prec / mp_bits_per_limb;
- if(val.right_data()[0]._mpfr_prec % mp_bits_per_limb)
+ if (val.right_data()[0]._mpfr_prec % mp_bits_per_limb)
++len;
- for(std::size_t i = 0; i < len; ++i)
+ for (std::size_t i = 0; i < len; ++i)
boost::hash_combine(result, val.right_data()[0]._mpfr_d[i]);
boost::hash_combine(result, val.right_data()[0]._mpfr_exp);
boost::hash_combine(result, val.right_data()[0]._mpfr_sign);
#ifdef BOOST_NO_SFINAE_EXPR
-namespace detail{
+namespace detail {
-template<unsigned D1, unsigned D2>
-struct is_explicitly_convertible<backends::mpfi_float_backend<D1>, backends::mpfi_float_backend<D2> > : public mpl::true_ {};
+template <unsigned D1, unsigned D2>
+struct is_explicitly_convertible<backends::mpfi_float_backend<D1>, backends::mpfi_float_backend<D2> > : public mpl::true_
+{};
-}
+} // namespace detail
#endif
-template<>
-struct number_category<detail::canonical<mpfi_t, backends::mpfi_float_backend<0> >::type> : public mpl::int_<number_kind_floating_point>{};
+namespace detail {
+template <>
+struct is_variable_precision<backends::mpfi_float_backend<0> > : public true_type
+{};
+} // namespace detail
+
+template <>
+struct number_category<detail::canonical<mpfi_t, backends::mpfi_float_backend<0> >::type> : public mpl::int_<number_kind_floating_point>
+{};
template <unsigned Digits10>
-struct is_interval_number<backends::mpfi_float_backend<Digits10> > : public mpl::true_ {};
+struct is_interval_number<backends::mpfi_float_backend<Digits10> > : public mpl::true_
+{};
using boost::multiprecision::backends::mpfi_float_backend;
-typedef number<mpfi_float_backend<50> > mpfi_float_50;
-typedef number<mpfi_float_backend<100> > mpfi_float_100;
-typedef number<mpfi_float_backend<500> > mpfi_float_500;
-typedef number<mpfi_float_backend<1000> > mpfi_float_1000;
-typedef number<mpfi_float_backend<0> > mpfi_float;
+typedef number<mpfi_float_backend<50> > mpfi_float_50;
+typedef number<mpfi_float_backend<100> > mpfi_float_100;
+typedef number<mpfi_float_backend<500> > mpfi_float_500;
+typedef number<mpfi_float_backend<1000> > mpfi_float_1000;
+typedef number<mpfi_float_backend<0> > mpfi_float;
//
// Special interval specific functions:
template <unsigned Digits10, expression_template_option ExpressionTemplates>
inline number<mpfr_float_backend<Digits10>, ExpressionTemplates> lower(const number<mpfi_float_backend<Digits10>, ExpressionTemplates>& val)
{
- number<mpfr_float_backend<Digits10> > result;
+ boost::multiprecision::detail::scoped_default_precision<number<mpfr_float_backend<Digits10>, ExpressionTemplates> > precision_guard(val);
+ number<mpfr_float_backend<Digits10> > result;
mpfr_set(result.backend().data(), val.backend().left_data(), GMP_RNDN);
return result;
}
template <unsigned Digits10, expression_template_option ExpressionTemplates>
inline number<mpfr_float_backend<Digits10>, ExpressionTemplates> upper(const number<mpfi_float_backend<Digits10>, ExpressionTemplates>& val)
{
- number<mpfr_float_backend<Digits10> > result;
+ boost::multiprecision::detail::scoped_default_precision<number<mpfr_float_backend<Digits10>, ExpressionTemplates> > precision_guard(val);
+ number<mpfr_float_backend<Digits10> > result;
mpfr_set(result.backend().data(), val.backend().right_data(), GMP_RNDN);
return result;
}
template <unsigned Digits10, expression_template_option ExpressionTemplates>
inline number<mpfr_float_backend<Digits10>, ExpressionTemplates> median(const number<mpfi_float_backend<Digits10>, ExpressionTemplates>& val)
{
- number<mpfr_float_backend<Digits10> > result;
+ boost::multiprecision::detail::scoped_default_precision<number<mpfr_float_backend<Digits10>, ExpressionTemplates> > precision_guard(val);
+ number<mpfr_float_backend<Digits10> > result;
mpfi_mid(result.backend().data(), val.backend().data());
return result;
}
template <unsigned Digits10, expression_template_option ExpressionTemplates>
inline number<mpfr_float_backend<Digits10>, ExpressionTemplates> width(const number<mpfi_float_backend<Digits10>, ExpressionTemplates>& val)
{
- number<mpfr_float_backend<Digits10> > result;
+ boost::multiprecision::detail::scoped_default_precision<number<mpfr_float_backend<Digits10>, ExpressionTemplates> > precision_guard(val);
+ number<mpfr_float_backend<Digits10> > result;
mpfi_diam_abs(result.backend().data(), val.backend().data());
return result;
}
template <unsigned Digits10, expression_template_option ExpressionTemplates>
-inline number<mpfi_float_backend<Digits10>, ExpressionTemplates> intersect(const number<mpfi_float_backend<Digits10>, ExpressionTemplates>& a, const number<mpfi_float_backend<Digits10>, ExpressionTemplates>& b)
+inline number<mpfi_float_backend<Digits10>, ExpressionTemplates> intersect(const number<mpfi_float_backend<Digits10>, ExpressionTemplates>& a, const number<mpfi_float_backend<Digits10>, ExpressionTemplates>& b)
{
- number<mpfi_float_backend<Digits10>, ExpressionTemplates> result;
+ boost::multiprecision::detail::scoped_default_precision<number<mpfi_float_backend<Digits10>, ExpressionTemplates> > precision_guard(a, b);
+ number<mpfi_float_backend<Digits10>, ExpressionTemplates> result;
mpfi_intersect(result.backend().data(), a.backend().data(), b.backend().data());
return result;
}
template <unsigned Digits10, expression_template_option ExpressionTemplates>
-inline number<mpfi_float_backend<Digits10>, ExpressionTemplates> hull(const number<mpfi_float_backend<Digits10>, ExpressionTemplates>& a, const number<mpfi_float_backend<Digits10>, ExpressionTemplates>& b)
+inline number<mpfi_float_backend<Digits10>, ExpressionTemplates> hull(const number<mpfi_float_backend<Digits10>, ExpressionTemplates>& a, const number<mpfi_float_backend<Digits10>, ExpressionTemplates>& b)
{
- number<mpfi_float_backend<Digits10>, ExpressionTemplates> result;
+ boost::multiprecision::detail::scoped_default_precision<number<mpfi_float_backend<Digits10>, ExpressionTemplates> > precision_guard(a, b);
+ number<mpfi_float_backend<Digits10>, ExpressionTemplates> result;
mpfi_union(result.backend().data(), a.backend().data(), b.backend().data());
return result;
}
template <unsigned Digits10, expression_template_option ExpressionTemplates>
-inline bool overlap(const number<mpfi_float_backend<Digits10>, ExpressionTemplates>& a, const number<mpfi_float_backend<Digits10>, ExpressionTemplates>& b)
+inline bool overlap(const number<mpfi_float_backend<Digits10>, ExpressionTemplates>& a, const number<mpfi_float_backend<Digits10>, ExpressionTemplates>& b)
{
- return (lower(a) <= lower(b) && lower(b) <= upper(a)) ||
- (lower(b) <= lower(a) && lower(a) <= upper(b));
+ return (lower(a) <= lower(b) && lower(b) <= upper(a)) ||
+ (lower(b) <= lower(a) && lower(a) <= upper(b));
}
template <unsigned Digits10, expression_template_option ExpressionTemplates1, expression_template_option ExpressionTemplates2>
-inline bool in(const number<mpfr_float_backend<Digits10>, ExpressionTemplates1>& a, const number<mpfi_float_backend<Digits10>, ExpressionTemplates2>& b)
+inline bool in(const number<mpfr_float_backend<Digits10>, ExpressionTemplates1>& a, const number<mpfi_float_backend<Digits10>, ExpressionTemplates2>& b)
{
- return mpfi_is_inside_fr(a.backend().data(), b.backend().data()) != 0;
+ return mpfi_is_inside_fr(a.backend().data(), b.backend().data()) != 0;
}
template <unsigned Digits10, expression_template_option ExpressionTemplates>
-inline bool zero_in(const number<mpfi_float_backend<Digits10>, ExpressionTemplates>& a)
+inline bool zero_in(const number<mpfi_float_backend<Digits10>, ExpressionTemplates>& a)
{
- return mpfi_has_zero(a.backend().data()) != 0;
+ return mpfi_has_zero(a.backend().data()) != 0;
}
template <unsigned Digits10, expression_template_option ExpressionTemplates>
-inline bool subset(const number<mpfi_float_backend<Digits10>, ExpressionTemplates>& a, const number<mpfi_float_backend<Digits10>, ExpressionTemplates>& b)
+inline bool subset(const number<mpfi_float_backend<Digits10>, ExpressionTemplates>& a, const number<mpfi_float_backend<Digits10>, ExpressionTemplates>& b)
{
- return mpfi_is_inside(a.backend().data(), b.backend().data()) != 0;
+ return mpfi_is_inside(a.backend().data(), b.backend().data()) != 0;
}
template <unsigned Digits10, expression_template_option ExpressionTemplates>
-inline bool proper_subset(const number<mpfi_float_backend<Digits10>, ExpressionTemplates>& a, const number<mpfi_float_backend<Digits10>, ExpressionTemplates>& b)
+inline bool proper_subset(const number<mpfi_float_backend<Digits10>, ExpressionTemplates>& a, const number<mpfi_float_backend<Digits10>, ExpressionTemplates>& b)
{
- return mpfi_is_strictly_inside(a.backend().data(), b.backend().data()) != 0;
+ return mpfi_is_strictly_inside(a.backend().data(), b.backend().data()) != 0;
}
template <unsigned Digits10, expression_template_option ExpressionTemplates>
-inline bool empty(const number<mpfi_float_backend<Digits10>, ExpressionTemplates>& a)
+inline bool empty(const number<mpfi_float_backend<Digits10>, ExpressionTemplates>& a)
{
- return mpfi_is_empty(a.backend().data()) != 0;
+ return mpfi_is_empty(a.backend().data()) != 0;
}
template <unsigned Digits10, expression_template_option ExpressionTemplates>
-inline bool singleton(const number<mpfi_float_backend<Digits10>, ExpressionTemplates>& a)
+inline bool singleton(const number<mpfi_float_backend<Digits10>, ExpressionTemplates>& a)
{
- return mpfr_cmp(a.backend().left_data(), a.backend().right_data()) == 0;
+ return mpfr_cmp(a.backend().left_data(), a.backend().right_data()) == 0;
}
template <unsigned Digits10, expression_template_option ExpressionTemplates>
template <unsigned Digits10, expression_template_option ExpressionTemplates>
inline boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> asinh BOOST_PREVENT_MACRO_SUBSTITUTION(const boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates>& arg)
{
+ boost::multiprecision::detail::scoped_default_precision<number<mpfi_float_backend<Digits10>, ExpressionTemplates> > precision_guard(arg);
+
boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> result;
mpfi_asinh(result.backend().data(), arg.backend().data());
- return BOOST_MP_MOVE(result);
+ return result;
}
template <unsigned Digits10, expression_template_option ExpressionTemplates>
inline boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> acosh BOOST_PREVENT_MACRO_SUBSTITUTION(const boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates>& arg)
{
+ boost::multiprecision::detail::scoped_default_precision<number<mpfi_float_backend<Digits10>, ExpressionTemplates> > precision_guard(arg);
+
boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> result;
mpfi_acosh(result.backend().data(), arg.backend().data());
- return BOOST_MP_MOVE(result);
+ return result;
}
template <unsigned Digits10, expression_template_option ExpressionTemplates>
inline boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> atanh BOOST_PREVENT_MACRO_SUBSTITUTION(const boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates>& arg)
{
+ boost::multiprecision::detail::scoped_default_precision<number<mpfi_float_backend<Digits10>, ExpressionTemplates> > precision_guard(arg);
+
boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> result;
mpfi_atanh(result.backend().data(), arg.backend().data());
- return BOOST_MP_MOVE(result);
+ return result;
}
template <unsigned Digits10, expression_template_option ExpressionTemplates>
inline boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> cbrt BOOST_PREVENT_MACRO_SUBSTITUTION(const boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates>& arg)
{
+ boost::multiprecision::detail::scoped_default_precision<number<mpfi_float_backend<Digits10>, ExpressionTemplates> > precision_guard(arg);
+
boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> result;
mpfi_cbrt(result.backend().data(), arg.backend().data());
- return BOOST_MP_MOVE(result);
+ return result;
}
template <unsigned Digits10, expression_template_option ExpressionTemplates>
inline boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> expm1 BOOST_PREVENT_MACRO_SUBSTITUTION(const boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates>& arg)
{
+ boost::multiprecision::detail::scoped_default_precision<number<mpfi_float_backend<Digits10>, ExpressionTemplates> > precision_guard(arg);
+
boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> result;
mpfi_expm1(result.backend().data(), arg.backend().data());
- return BOOST_MP_MOVE(result);
+ return result;
}
template <unsigned Digits10, expression_template_option ExpressionTemplates>
inline boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> log1p BOOST_PREVENT_MACRO_SUBSTITUTION(const boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates>& arg)
{
+ boost::multiprecision::detail::scoped_default_precision<number<mpfi_float_backend<Digits10>, ExpressionTemplates> > precision_guard(arg);
+
boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> result;
mpfi_log1p(result.backend().data(), arg.backend().data());
- return BOOST_MP_MOVE(result);
+ return result;
}
-
} // namespace multiprecision
-namespace math{
+namespace math {
-namespace tools{
+namespace tools {
template <>
inline int digits<boost::multiprecision::mpfi_float>()
#ifdef BOOST_MATH_NOEXCEPT
- BOOST_NOEXCEPT
+ BOOST_NOEXCEPT
#endif
{
- return multiprecision::detail::digits10_2_2(boost::multiprecision::mpfi_float::default_precision());
+ return multiprecision::detail::digits10_2_2(boost::multiprecision::mpfi_float::default_precision());
}
template <>
inline int digits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<0>, boost::multiprecision::et_off> >()
#ifdef BOOST_MATH_NOEXCEPT
- BOOST_NOEXCEPT
+ BOOST_NOEXCEPT
#endif
{
- return multiprecision::detail::digits10_2_2(boost::multiprecision::mpfi_float::default_precision());
+ return multiprecision::detail::digits10_2_2(boost::multiprecision::mpfi_float::default_precision());
}
template <>
}
// mpfi gets used with logged_adaptor fairly often, so specialize for that use case as well:
-typedef boost::multiprecision::number<boost::multiprecision::backends::logged_adaptor<boost::multiprecision::mpfi_float::backend_type>, boost::multiprecision::et_on> logged_type1;
+typedef boost::multiprecision::number<boost::multiprecision::backends::logged_adaptor<boost::multiprecision::mpfi_float::backend_type>, boost::multiprecision::et_on> logged_type1;
typedef boost::multiprecision::number<boost::multiprecision::backends::logged_adaptor<boost::multiprecision::mpfi_float::backend_type>, boost::multiprecision::et_off> logged_type2;
template <>
inline int digits<logged_type1>()
#ifdef BOOST_MATH_NOEXCEPT
-BOOST_NOEXCEPT
+ BOOST_NOEXCEPT
#endif
{
- return multiprecision::detail::digits10_2_2(logged_type1::default_precision());
+ return multiprecision::detail::digits10_2_2(logged_type1::default_precision());
}
template <>
-inline int digits<logged_type2 >()
+inline int digits<logged_type2>()
#ifdef BOOST_MATH_NOEXCEPT
-BOOST_NOEXCEPT
+ BOOST_NOEXCEPT
#endif
{
- return multiprecision::detail::digits10_2_2(logged_type1::default_precision());
+ return multiprecision::detail::digits10_2_2(logged_type1::default_precision());
}
template <>
template <>
inline logged_type2
-max_value<logged_type2 >()
+max_value<logged_type2>()
{
logged_type2 result(0.5);
mpfi_mul_2exp(result.backend().value().data(), result.backend().value().data(), mpfr_get_emax());
template <>
inline logged_type2
-min_value<logged_type2 >()
+min_value<logged_type2>()
{
logged_type2 result(0.5);
mpfi_div_2exp(result.backend().value().data(), result.backend().value().data(), -mpfr_get_emin());
}
} // namespace tools
-namespace constants{ namespace detail{
+namespace constants { namespace detail {
-template <class T> struct constant_pi;
-template <class T> struct constant_ln_two;
-template <class T> struct constant_euler;
-template <class T> struct constant_catalan;
+template <class T>
+struct constant_pi;
+template <class T>
+struct constant_ln_two;
+template <class T>
+struct constant_euler;
+template <class T>
+struct constant_catalan;
//
// Initializer: ensure all our constants are initialized prior to the first call of main:
boost::math::constants::euler<T>();
boost::math::constants::catalan<T>();
}
- void force_instantiate()const{}
+ void force_instantiate() const {}
};
static const init initializer;
- static void force_instantiate()
+ static void force_instantiate()
{
initializer.force_instantiate();
}
template <class T>
const typename mpfi_initializer<T>::init mpfi_initializer<T>::initializer;
-template<unsigned Digits10, boost::multiprecision::expression_template_option ExpressionTemplates>
+template <unsigned Digits10, boost::multiprecision::expression_template_option ExpressionTemplates>
struct constant_pi<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> >
{
typedef boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> result_type;
- template<int N>
+ template <int N>
static inline const result_type& get(const mpl::int_<N>&)
{
mpfi_initializer<result_type>::force_instantiate();
static result_type result;
- static bool has_init = false;
- if(!has_init)
+ static bool has_init = false;
+ if (!has_init)
{
has_init = true;
mpfi_const_pi(result.backend().data());
return result;
}
};
-template<unsigned Digits10, boost::multiprecision::expression_template_option ExpressionTemplates>
+template <unsigned Digits10, boost::multiprecision::expression_template_option ExpressionTemplates>
struct constant_ln_two<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> >
{
typedef boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> result_type;
- template<int N>
+ template <int N>
static inline const result_type& get(const mpl::int_<N>&)
{
mpfi_initializer<result_type>::force_instantiate();
static result_type result;
- static bool has_init = false;
- if(!has_init)
+ static bool has_init = false;
+ if (!has_init)
{
has_init = true;
mpfi_const_log2(result.backend().data());
return result;
}
};
-template<unsigned Digits10, boost::multiprecision::expression_template_option ExpressionTemplates>
+template <unsigned Digits10, boost::multiprecision::expression_template_option ExpressionTemplates>
struct constant_euler<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> >
{
typedef boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> result_type;
- template<int N>
+ template <int N>
static inline result_type const& get(const mpl::int_<N>&)
{
mpfi_initializer<result_type>::force_instantiate();
static result_type result;
- static bool has_init = false;
- if(!has_init)
+ static bool has_init = false;
+ if (!has_init)
{
has_init = true;
mpfi_const_euler(result.backend().data());
return result;
}
};
-template<unsigned Digits10, boost::multiprecision::expression_template_option ExpressionTemplates>
+template <unsigned Digits10, boost::multiprecision::expression_template_option ExpressionTemplates>
struct constant_catalan<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> >
{
typedef boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> result_type;
- template<int N>
+ template <int N>
static inline result_type const& get(const mpl::int_<N>&)
{
mpfi_initializer<result_type>::force_instantiate();
static result_type result;
- static bool has_init = false;
- if(!has_init)
+ static bool has_init = false;
+ if (!has_init)
{
has_init = true;
mpfi_const_catalan(result.backend().data());
}
};
-}} // namespaces
+}} // namespace constants::detail
-}} // namespaces
+} // namespace math
+} // namespace boost
-namespace std{
+namespace std {
//
// numeric_limits [partial] specializations for the types declared in this header:
//
-template<unsigned Digits10, boost::multiprecision::expression_template_option ExpressionTemplates>
+template <unsigned Digits10, boost::multiprecision::expression_template_option ExpressionTemplates>
class numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> >
{
typedef boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> number_type;
-public:
+
+ public:
BOOST_STATIC_CONSTEXPR bool is_specialized = true;
- static number_type (min)()
+ static number_type(min)()
{
initializer.do_nothing();
static std::pair<bool, number_type> value;
- if(!value.first)
+ if (!value.first)
{
- value.first = true;
+ value.first = true;
value.second = 0.5;
mpfi_div_2exp(value.second.backend().data(), value.second.backend().data(), -mpfr_get_emin());
}
return value.second;
}
- static number_type (max)()
+ static number_type(max)()
{
initializer.do_nothing();
static std::pair<bool, number_type> value;
- if(!value.first)
+ if (!value.first)
{
- value.first = true;
+ value.first = true;
value.second = 0.5;
mpfi_mul_2exp(value.second.backend().data(), value.second.backend().data(), mpfr_get_emax());
}
{
return -(max)();
}
- BOOST_STATIC_CONSTEXPR int digits = static_cast<int>((Digits10 * 1000L) / 301L + ((Digits10 * 1000L) % 301 ? 2 : 1));
+ BOOST_STATIC_CONSTEXPR int digits = static_cast<int>((Digits10 * 1000L) / 301L + ((Digits10 * 1000L) % 301 ? 2 : 1));
BOOST_STATIC_CONSTEXPR int digits10 = Digits10;
// Is this really correct???
- BOOST_STATIC_CONSTEXPR int max_digits10 = Digits10 + 3;
- BOOST_STATIC_CONSTEXPR bool is_signed = true;
- BOOST_STATIC_CONSTEXPR bool is_integer = false;
- BOOST_STATIC_CONSTEXPR bool is_exact = false;
- BOOST_STATIC_CONSTEXPR int radix = 2;
- static number_type epsilon()
+ BOOST_STATIC_CONSTEXPR int max_digits10 = Digits10 + 3;
+ BOOST_STATIC_CONSTEXPR bool is_signed = true;
+ BOOST_STATIC_CONSTEXPR bool is_integer = false;
+ BOOST_STATIC_CONSTEXPR bool is_exact = false;
+ BOOST_STATIC_CONSTEXPR int radix = 2;
+ static number_type epsilon()
{
initializer.do_nothing();
static std::pair<bool, number_type> value;
- if(!value.first)
+ if (!value.first)
{
- value.first = true;
+ value.first = true;
value.second = 1;
mpfi_div_2exp(value.second.backend().data(), value.second.backend().data(), std::numeric_limits<number_type>::digits - 1);
}
// returns epsilon/2
initializer.do_nothing();
static std::pair<bool, number_type> value;
- if(!value.first)
+ if (!value.first)
{
- value.first = true;
+ value.first = true;
value.second = 1;
mpfi_div_2exp(value.second.backend().data(), value.second.backend().data(), 1);
}
return value.second;
}
- BOOST_STATIC_CONSTEXPR long min_exponent = MPFR_EMIN_DEFAULT;
- BOOST_STATIC_CONSTEXPR long min_exponent10 = (MPFR_EMIN_DEFAULT / 1000) * 301L;
- BOOST_STATIC_CONSTEXPR long max_exponent = MPFR_EMAX_DEFAULT;
- BOOST_STATIC_CONSTEXPR long max_exponent10 = (MPFR_EMAX_DEFAULT / 1000) * 301L;
- BOOST_STATIC_CONSTEXPR bool has_infinity = true;
- BOOST_STATIC_CONSTEXPR bool has_quiet_NaN = true;
- BOOST_STATIC_CONSTEXPR bool has_signaling_NaN = false;
- BOOST_STATIC_CONSTEXPR float_denorm_style has_denorm = denorm_absent;
- BOOST_STATIC_CONSTEXPR bool has_denorm_loss = false;
- static number_type infinity()
+ BOOST_STATIC_CONSTEXPR long min_exponent = MPFR_EMIN_DEFAULT;
+ BOOST_STATIC_CONSTEXPR long min_exponent10 = (MPFR_EMIN_DEFAULT / 1000) * 301L;
+ BOOST_STATIC_CONSTEXPR long max_exponent = MPFR_EMAX_DEFAULT;
+ BOOST_STATIC_CONSTEXPR long max_exponent10 = (MPFR_EMAX_DEFAULT / 1000) * 301L;
+ BOOST_STATIC_CONSTEXPR bool has_infinity = true;
+ BOOST_STATIC_CONSTEXPR bool has_quiet_NaN = true;
+ BOOST_STATIC_CONSTEXPR bool has_signaling_NaN = false;
+ BOOST_STATIC_CONSTEXPR float_denorm_style has_denorm = denorm_absent;
+ BOOST_STATIC_CONSTEXPR bool has_denorm_loss = false;
+ static number_type infinity()
{
initializer.do_nothing();
static std::pair<bool, number_type> value;
- if(!value.first)
+ if (!value.first)
{
boost::multiprecision::mpfr_float_backend<Digits10> t;
mpfr_set_inf(t.data(), 1);
{
initializer.do_nothing();
static std::pair<bool, number_type> value;
- if(!value.first)
+ if (!value.first)
{
boost::multiprecision::mpfr_float_backend<Digits10> t;
mpfr_set_nan(t.data());
return number_type(0);
}
BOOST_STATIC_CONSTEXPR number_type denorm_min() { return number_type(0); }
- BOOST_STATIC_CONSTEXPR bool is_iec559 = false;
- BOOST_STATIC_CONSTEXPR bool is_bounded = true;
- BOOST_STATIC_CONSTEXPR bool is_modulo = false;
- BOOST_STATIC_CONSTEXPR bool traps = true;
- BOOST_STATIC_CONSTEXPR bool tinyness_before = false;
+ BOOST_STATIC_CONSTEXPR bool is_iec559 = false;
+ BOOST_STATIC_CONSTEXPR bool is_bounded = true;
+ BOOST_STATIC_CONSTEXPR bool is_modulo = false;
+ BOOST_STATIC_CONSTEXPR bool traps = true;
+ BOOST_STATIC_CONSTEXPR bool tinyness_before = false;
BOOST_STATIC_CONSTEXPR float_round_style round_style = round_to_nearest;
-private:
+ private:
struct data_initializer
{
data_initializer()
std::numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<digits10> > >::infinity();
std::numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<digits10> > >::quiet_NaN();
}
- void do_nothing()const{}
+ void do_nothing() const {}
};
static const data_initializer initializer;
};
-template<unsigned Digits10, boost::multiprecision::expression_template_option ExpressionTemplates>
+template <unsigned Digits10, boost::multiprecision::expression_template_option ExpressionTemplates>
const typename numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> >::data_initializer numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> >::initializer;
#ifndef BOOST_NO_INCLASS_MEMBER_INITIALIZATION
#endif
-
-template<boost::multiprecision::expression_template_option ExpressionTemplates>
+template <boost::multiprecision::expression_template_option ExpressionTemplates>
class numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<0>, ExpressionTemplates> >
{
typedef boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<0>, ExpressionTemplates> number_type;
-public:
+
+ public:
BOOST_STATIC_CONSTEXPR bool is_specialized = false;
- static number_type (min)() { return number_type(0); }
- static number_type (max)() { return number_type(0); }
- static number_type lowest() { return number_type(0); }
- BOOST_STATIC_CONSTEXPR int digits = 0;
- BOOST_STATIC_CONSTEXPR int digits10 = 0;
- BOOST_STATIC_CONSTEXPR int max_digits10 = 0;
- BOOST_STATIC_CONSTEXPR bool is_signed = false;
- BOOST_STATIC_CONSTEXPR bool is_integer = false;
- BOOST_STATIC_CONSTEXPR bool is_exact = false;
- BOOST_STATIC_CONSTEXPR int radix = 0;
- static number_type epsilon() { return number_type(0); }
- static number_type round_error() { return number_type(0); }
- BOOST_STATIC_CONSTEXPR int min_exponent = 0;
- BOOST_STATIC_CONSTEXPR int min_exponent10 = 0;
- BOOST_STATIC_CONSTEXPR int max_exponent = 0;
- BOOST_STATIC_CONSTEXPR int max_exponent10 = 0;
- BOOST_STATIC_CONSTEXPR bool has_infinity = false;
- BOOST_STATIC_CONSTEXPR bool has_quiet_NaN = false;
- BOOST_STATIC_CONSTEXPR bool has_signaling_NaN = false;
- BOOST_STATIC_CONSTEXPR float_denorm_style has_denorm = denorm_absent;
- BOOST_STATIC_CONSTEXPR bool has_denorm_loss = false;
- static number_type infinity() { return number_type(0); }
- static number_type quiet_NaN() { return number_type(0); }
- static number_type signaling_NaN() { return number_type(0); }
- static number_type denorm_min() { return number_type(0); }
- BOOST_STATIC_CONSTEXPR bool is_iec559 = false;
- BOOST_STATIC_CONSTEXPR bool is_bounded = false;
- BOOST_STATIC_CONSTEXPR bool is_modulo = false;
- BOOST_STATIC_CONSTEXPR bool traps = false;
- BOOST_STATIC_CONSTEXPR bool tinyness_before = false;
- BOOST_STATIC_CONSTEXPR float_round_style round_style = round_toward_zero;
+ static number_type(min)() { return number_type(0); }
+ static number_type(max)() { return number_type(0); }
+ static number_type lowest() { return number_type(0); }
+ BOOST_STATIC_CONSTEXPR int digits = 0;
+ BOOST_STATIC_CONSTEXPR int digits10 = 0;
+ BOOST_STATIC_CONSTEXPR int max_digits10 = 0;
+ BOOST_STATIC_CONSTEXPR bool is_signed = false;
+ BOOST_STATIC_CONSTEXPR bool is_integer = false;
+ BOOST_STATIC_CONSTEXPR bool is_exact = false;
+ BOOST_STATIC_CONSTEXPR int radix = 0;
+ static number_type epsilon() { return number_type(0); }
+ static number_type round_error() { return number_type(0); }
+ BOOST_STATIC_CONSTEXPR int min_exponent = 0;
+ BOOST_STATIC_CONSTEXPR int min_exponent10 = 0;
+ BOOST_STATIC_CONSTEXPR int max_exponent = 0;
+ BOOST_STATIC_CONSTEXPR int max_exponent10 = 0;
+ BOOST_STATIC_CONSTEXPR bool has_infinity = false;
+ BOOST_STATIC_CONSTEXPR bool has_quiet_NaN = false;
+ BOOST_STATIC_CONSTEXPR bool has_signaling_NaN = false;
+ BOOST_STATIC_CONSTEXPR float_denorm_style has_denorm = denorm_absent;
+ BOOST_STATIC_CONSTEXPR bool has_denorm_loss = false;
+ static number_type infinity() { return number_type(0); }
+ static number_type quiet_NaN() { return number_type(0); }
+ static number_type signaling_NaN() { return number_type(0); }
+ static number_type denorm_min() { return number_type(0); }
+ BOOST_STATIC_CONSTEXPR bool is_iec559 = false;
+ BOOST_STATIC_CONSTEXPR bool is_bounded = false;
+ BOOST_STATIC_CONSTEXPR bool is_modulo = false;
+ BOOST_STATIC_CONSTEXPR bool traps = false;
+ BOOST_STATIC_CONSTEXPR bool tinyness_before = false;
+ BOOST_STATIC_CONSTEXPR float_round_style round_style = round_toward_zero;
};
#ifndef BOOST_NO_INCLASS_MEMBER_INITIALIZATION