1 ///////////////////////////////////////////////////////////////////////////////
2 // Copyright 2011 John Maddock. Distributed under the Boost
3 // Software License, Version 1.0. (See accompanying file
4 // LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)
6 #ifndef BOOST_MATH_BN_MPFI_HPP
7 #define BOOST_MATH_BN_MPFI_HPP
9 #include <boost/multiprecision/number.hpp>
10 #include <boost/math/special_functions/fpclassify.hpp>
11 #include <boost/cstdint.hpp>
12 #include <boost/multiprecision/detail/big_lanczos.hpp>
13 #include <boost/multiprecision/detail/digits.hpp>
14 #include <boost/multiprecision/mpfr.hpp>
15 #include <boost/multiprecision/logged_adaptor.hpp>
16 #include <boost/math/constants/constants.hpp>
17 #include <boost/functional/hash_fwd.hpp>
22 #ifndef BOOST_MULTIPRECISION_MPFI_DEFAULT_PRECISION
23 #define BOOST_MULTIPRECISION_MPFI_DEFAULT_PRECISION 20
27 namespace multiprecision {
30 template <unsigned digits10>
31 struct mpfi_float_backend;
33 } // namespace backends
35 template <unsigned digits10>
36 struct number_category<backends::mpfi_float_backend<digits10> > : public mpl::int_<number_kind_floating_point>
39 struct interval_error : public std::runtime_error
41 interval_error(const std::string& s) : std::runtime_error(s) {}
48 inline int mpfi_sgn(mpfi_srcptr p)
52 if (mpfi_is_strictly_pos(p))
54 if (mpfi_is_strictly_neg(p))
56 BOOST_THROW_EXCEPTION(interval_error("Sign of interval is ambiguous."));
59 template <unsigned digits10>
60 struct mpfi_float_imp;
62 template <unsigned digits10>
65 #ifdef BOOST_HAS_LONG_LONG
66 typedef mpl::list<long, boost::long_long_type> signed_types;
67 typedef mpl::list<unsigned long, boost::ulong_long_type> unsigned_types;
69 typedef mpl::list<long> signed_types;
70 typedef mpl::list<unsigned long> unsigned_types;
72 typedef mpl::list<double, long double> float_types;
73 typedef long exponent_type;
77 mpfi_init2(m_data, multiprecision::detail::digits10_2_2(digits10 ? digits10 : get_default_precision()));
78 mpfi_set_ui(m_data, 0u);
80 mpfi_float_imp(unsigned prec)
82 mpfi_init2(m_data, prec);
83 mpfi_set_ui(m_data, 0u);
86 mpfi_float_imp(const mpfi_float_imp& o)
88 mpfi_init2(m_data, mpfi_get_prec(o.data()));
89 if (o.m_data[0].left._mpfr_d)
90 mpfi_set(m_data, o.m_data);
92 #ifndef BOOST_NO_CXX11_RVALUE_REFERENCES
93 mpfi_float_imp(mpfi_float_imp&& o) BOOST_NOEXCEPT
95 m_data[0] = o.m_data[0];
96 o.m_data[0].left._mpfr_d = 0;
99 mpfi_float_imp& operator=(const mpfi_float_imp& o)
101 if (m_data[0].left._mpfr_d == 0)
102 mpfi_init2(m_data, mpfi_get_prec(o.data()));
103 if (mpfi_get_prec(o.data()) != mpfi_get_prec(data()))
105 mpfi_float_imp t(mpfi_get_prec(o.data()));
111 if (o.m_data[0].left._mpfr_d)
112 mpfi_set(m_data, o.m_data);
116 #ifndef BOOST_NO_CXX11_RVALUE_REFERENCES
117 mpfi_float_imp& operator=(mpfi_float_imp&& o) BOOST_NOEXCEPT
119 mpfi_swap(m_data, o.m_data);
123 #ifdef BOOST_HAS_LONG_LONG
124 #ifdef _MPFR_H_HAVE_INTMAX_T
125 mpfi_float_imp& operator=(boost::ulong_long_type i)
127 if (m_data[0].left._mpfr_d == 0)
128 mpfi_init2(m_data, multiprecision::detail::digits10_2_2(digits10 ? digits10 : get_default_precision()));
129 mpfr_set_uj(left_data(), i, GMP_RNDD);
130 mpfr_set_uj(right_data(), i, GMP_RNDU);
133 mpfi_float_imp& operator=(boost::long_long_type i)
135 if (m_data[0].left._mpfr_d == 0)
136 mpfi_init2(m_data, multiprecision::detail::digits10_2_2(digits10 ? digits10 : get_default_precision()));
137 mpfr_set_sj(left_data(), i, GMP_RNDD);
138 mpfr_set_sj(right_data(), i, GMP_RNDU);
142 mpfi_float_imp& operator=(boost::ulong_long_type i)
144 if (m_data[0].left._mpfr_d == 0)
145 mpfi_init2(m_data, multiprecision::detail::digits10_2_2(digits10 ? digits10 : get_default_precision()));
146 boost::ulong_long_type mask = ((((1uLL << (std::numeric_limits<unsigned long>::digits - 1)) - 1) << 1) | 1u);
149 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))));
150 mpfi_set_ui(m_data, 0);
153 mpfi_set_ui(t, static_cast<unsigned long>(i & mask));
155 mpfi_mul_2exp(t, t, shift);
156 mpfi_add(m_data, m_data, t);
157 shift += std::numeric_limits<unsigned long>::digits;
158 i >>= std::numeric_limits<unsigned long>::digits;
163 mpfi_float_imp& operator=(boost::long_long_type i)
165 if (m_data[0].left._mpfr_d == 0)
166 mpfi_init2(m_data, multiprecision::detail::digits10_2_2(digits10 ? digits10 : get_default_precision()));
168 *this = boost::multiprecision::detail::unsigned_abs(i);
170 mpfi_neg(m_data, m_data);
175 mpfi_float_imp& operator=(unsigned long i)
177 if (m_data[0].left._mpfr_d == 0)
178 mpfi_init2(m_data, multiprecision::detail::digits10_2_2(digits10 ? digits10 : get_default_precision()));
179 mpfi_set_ui(m_data, i);
182 mpfi_float_imp& operator=(long i)
184 if (m_data[0].left._mpfr_d == 0)
185 mpfi_init2(m_data, multiprecision::detail::digits10_2_2(digits10 ? digits10 : get_default_precision()));
186 mpfi_set_si(m_data, i);
189 mpfi_float_imp& operator=(double d)
191 if (m_data[0].left._mpfr_d == 0)
192 mpfi_init2(m_data, multiprecision::detail::digits10_2_2(digits10 ? digits10 : get_default_precision()));
193 mpfi_set_d(m_data, d);
196 mpfi_float_imp& operator=(long double a)
198 if (m_data[0].left._mpfr_d == 0)
199 mpfi_init2(m_data, multiprecision::detail::digits10_2_2(digits10 ? digits10 : get_default_precision()));
200 mpfr_set_ld(left_data(), a, GMP_RNDD);
201 mpfr_set_ld(right_data(), a, GMP_RNDU);
204 mpfi_float_imp& operator=(const char* s)
206 using default_ops::eval_fpclassify;
208 if (m_data[0].left._mpfr_d == 0)
209 mpfi_init2(m_data, multiprecision::detail::digits10_2_2(digits10 ? digits10 : get_default_precision()));
211 if (s && (*s == '{'))
213 mpfr_float_backend<digits10> a, b;
216 while (*p && (*p != ',') && (*p != '}'))
218 part.assign(s + 1, p);
221 if (*p && (*p != '}'))
224 while (*p && (*p != ',') && (*p != '}'))
226 part.assign(s + 1, p);
232 if (eval_fpclassify(a) == (int)FP_NAN)
234 mpfi_set_fr(this->data(), a.data());
236 else if (eval_fpclassify(b) == (int)FP_NAN)
238 mpfi_set_fr(this->data(), b.data());
242 if (a.compare(b) > 0)
244 BOOST_THROW_EXCEPTION(std::runtime_error("Attempt to create interval with invalid range (start is greater than end)."));
246 mpfi_interv_fr(m_data, a.data(), b.data());
249 else if (mpfi_set_str(m_data, s, 10) != 0)
251 BOOST_THROW_EXCEPTION(std::runtime_error(std::string("Unable to parse string \"") + s + std::string("\"as a valid floating point number.")));
255 void swap(mpfi_float_imp& o) BOOST_NOEXCEPT
257 mpfi_swap(m_data, o.m_data);
259 std::string str(std::streamsize digits, std::ios_base::fmtflags f) const
261 BOOST_ASSERT(m_data[0].left._mpfr_d);
263 mpfr_float_backend<digits10> a, b;
265 mpfi_get_left(a.data(), m_data);
266 mpfi_get_right(b.data(), m_data);
268 if (a.compare(b) == 0)
269 return a.str(digits, f);
271 return "{" + a.str(digits, f) + "," + b.str(digits, f) + "}";
273 ~mpfi_float_imp() BOOST_NOEXCEPT
275 if (m_data[0].left._mpfr_d)
278 void negate() BOOST_NOEXCEPT
280 BOOST_ASSERT(m_data[0].left._mpfr_d);
281 mpfi_neg(m_data, m_data);
283 int compare(const mpfi_float_imp& o) const BOOST_NOEXCEPT
285 BOOST_ASSERT(m_data[0].left._mpfr_d && o.m_data[0].left._mpfr_d);
286 if (mpfr_cmp(right_data(), o.left_data()) < 0)
288 if (mpfr_cmp(left_data(), o.right_data()) > 0)
290 if ((mpfr_cmp(left_data(), o.left_data()) == 0) && (mpfr_cmp(right_data(), o.right_data()) == 0))
292 BOOST_THROW_EXCEPTION(interval_error("Ambiguous comparison between two values."));
296 int compare(V v) const BOOST_NOEXCEPT
302 mpfi_t& data() BOOST_NOEXCEPT
304 BOOST_ASSERT(m_data[0].left._mpfr_d);
307 const mpfi_t& data() const BOOST_NOEXCEPT
309 BOOST_ASSERT(m_data[0].left._mpfr_d);
312 mpfr_ptr left_data() BOOST_NOEXCEPT
314 BOOST_ASSERT(m_data[0].left._mpfr_d);
315 return &(m_data[0].left);
317 mpfr_srcptr left_data() const BOOST_NOEXCEPT
319 BOOST_ASSERT(m_data[0].left._mpfr_d);
320 return &(m_data[0].left);
322 mpfr_ptr right_data() BOOST_NOEXCEPT
324 BOOST_ASSERT(m_data[0].left._mpfr_d);
325 return &(m_data[0].right);
327 mpfr_srcptr right_data() const BOOST_NOEXCEPT
329 BOOST_ASSERT(m_data[0].left._mpfr_d);
330 return &(m_data[0].right);
335 static unsigned& get_default_precision() BOOST_NOEXCEPT
337 static unsigned val = BOOST_MULTIPRECISION_MPFI_DEFAULT_PRECISION;
342 } // namespace detail
344 template <unsigned digits10>
345 struct mpfi_float_backend : public detail::mpfi_float_imp<digits10>
347 mpfi_float_backend() : detail::mpfi_float_imp<digits10>() {}
348 mpfi_float_backend(const mpfi_float_backend& o) : detail::mpfi_float_imp<digits10>(o) {}
349 #ifndef BOOST_NO_CXX11_RVALUE_REFERENCES
350 mpfi_float_backend(mpfi_float_backend&& o) : detail::mpfi_float_imp<digits10>(static_cast<detail::mpfi_float_imp<digits10>&&>(o))
353 template <unsigned D>
354 mpfi_float_backend(const mpfi_float_backend<D>& val, typename enable_if_c<D <= digits10>::type* = 0)
355 : detail::mpfi_float_imp<digits10>()
357 mpfi_set(this->m_data, val.data());
359 template <unsigned D>
360 explicit mpfi_float_backend(const mpfi_float_backend<D>& val, typename disable_if_c<D <= digits10>::type* = 0)
361 : detail::mpfi_float_imp<digits10>()
363 mpfi_set(this->m_data, val.data());
365 mpfi_float_backend(const mpfi_t val)
366 : detail::mpfi_float_imp<digits10>()
368 mpfi_set(this->m_data, val);
370 mpfi_float_backend& operator=(const mpfi_float_backend& o)
372 *static_cast<detail::mpfi_float_imp<digits10>*>(this) = static_cast<detail::mpfi_float_imp<digits10> const&>(o);
375 template <unsigned D>
376 mpfi_float_backend(const mpfr_float_backend<D>& val, typename enable_if_c<D <= digits10>::type* = 0)
377 : detail::mpfi_float_imp<digits10>()
379 mpfi_set_fr(this->m_data, val.data());
381 template <unsigned D>
382 mpfi_float_backend& operator=(const mpfr_float_backend<D>& val)
384 mpfi_set_fr(this->m_data, val.data());
387 template <unsigned D>
388 explicit mpfi_float_backend(const mpfr_float_backend<D>& val, typename disable_if_c<D <= digits10>::type* = 0)
389 : detail::mpfi_float_imp<digits10>()
391 mpfi_set_fr(this->m_data, val.data());
393 #ifndef BOOST_NO_CXX11_RVALUE_REFERENCES
394 mpfi_float_backend& operator=(mpfi_float_backend&& o) BOOST_NOEXCEPT
396 *static_cast<detail::mpfi_float_imp<digits10>*>(this) = static_cast<detail::mpfi_float_imp<digits10>&&>(o);
401 mpfi_float_backend& operator=(const V& v)
403 *static_cast<detail::mpfi_float_imp<digits10>*>(this) = v;
406 mpfi_float_backend& operator=(const mpfi_t val)
408 mpfi_set(this->m_data, val);
411 // We don't change our precision here, this is a fixed precision type:
412 template <unsigned D>
413 mpfi_float_backend& operator=(const mpfi_float_backend<D>& val)
415 mpfi_set(this->m_data, val.data());
420 template <unsigned D1, unsigned D2, mpfr_allocation_type AllocationType>
421 void assign_components(mpfi_float_backend<D1>& result, const mpfr_float_backend<D2, AllocationType>& a, const mpfr_float_backend<D2, AllocationType>& b);
423 template <unsigned Digits10, class V>
424 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
425 assign_components(mpfi_float_backend<Digits10>& result, const V& a, const V& b);
428 struct mpfi_float_backend<0> : public detail::mpfi_float_imp<0>
430 mpfi_float_backend() : detail::mpfi_float_imp<0>() {}
431 mpfi_float_backend(const mpfi_t val)
432 : detail::mpfi_float_imp<0>(mpfi_get_prec(val))
434 mpfi_set(this->m_data, val);
436 mpfi_float_backend(const mpfi_float_backend& o) : detail::mpfi_float_imp<0>(o) {}
437 #ifndef BOOST_NO_CXX11_RVALUE_REFERENCES
438 mpfi_float_backend(mpfi_float_backend&& o) BOOST_NOEXCEPT : detail::mpfi_float_imp<0>(static_cast<detail::mpfi_float_imp<0>&&>(o))
441 mpfi_float_backend(const mpfi_float_backend& o, unsigned digits10)
442 : detail::mpfi_float_imp<0>(multiprecision::detail::digits10_2_2(digits10))
444 mpfi_set(this->m_data, o.data());
447 mpfi_float_backend(const V& a, const V& b, unsigned digits10)
448 : detail::mpfi_float_imp<0>(multiprecision::detail::digits10_2_2(digits10))
450 assign_components(*this, a, b);
453 template <unsigned D>
454 mpfi_float_backend(const mpfi_float_backend<D>& val)
455 : detail::mpfi_float_imp<0>(mpfi_get_prec(val.data()))
457 mpfi_set(this->m_data, val.data());
459 mpfi_float_backend& operator=(const mpfi_float_backend& o)
461 mpfi_set_prec(this->m_data, mpfi_get_prec(o.data()));
462 mpfi_set(this->m_data, o.data());
465 #ifndef BOOST_NO_CXX11_RVALUE_REFERENCES
466 mpfi_float_backend& operator=(mpfi_float_backend&& o) BOOST_NOEXCEPT
468 *static_cast<detail::mpfi_float_imp<0>*>(this) = static_cast<detail::mpfi_float_imp<0>&&>(o);
473 mpfi_float_backend& operator=(const V& v)
475 *static_cast<detail::mpfi_float_imp<0>*>(this) = v;
478 mpfi_float_backend& operator=(const mpfi_t val)
480 mpfi_set_prec(this->m_data, mpfi_get_prec(val));
481 mpfi_set(this->m_data, val);
484 template <unsigned D>
485 mpfi_float_backend& operator=(const mpfi_float_backend<D>& val)
487 mpfi_set_prec(this->m_data, mpfi_get_prec(val.data()));
488 mpfi_set(this->m_data, val.data());
491 static unsigned default_precision() BOOST_NOEXCEPT
493 return get_default_precision();
495 static void default_precision(unsigned v) BOOST_NOEXCEPT
497 get_default_precision() = v;
499 unsigned precision() const BOOST_NOEXCEPT
501 return multiprecision::detail::digits2_2_10(mpfi_get_prec(this->m_data));
503 void precision(unsigned digits10) BOOST_NOEXCEPT
505 mpfi_float_backend t(*this, digits10);
510 template <unsigned digits10, class T>
511 inline typename enable_if<is_arithmetic<T>, bool>::type eval_eq(const mpfi_float_backend<digits10>& a, const T& b) BOOST_NOEXCEPT
513 return a.compare(b) == 0;
515 template <unsigned digits10, class T>
516 inline typename enable_if<is_arithmetic<T>, bool>::type eval_lt(const mpfi_float_backend<digits10>& a, const T& b) BOOST_NOEXCEPT
518 return a.compare(b) < 0;
520 template <unsigned digits10, class T>
521 inline typename enable_if<is_arithmetic<T>, bool>::type eval_gt(const mpfi_float_backend<digits10>& a, const T& b) BOOST_NOEXCEPT
523 return a.compare(b) > 0;
526 template <unsigned D1, unsigned D2>
527 inline void eval_add(mpfi_float_backend<D1>& result, const mpfi_float_backend<D2>& o)
529 mpfi_add(result.data(), result.data(), o.data());
531 template <unsigned D1, unsigned D2>
532 inline void eval_subtract(mpfi_float_backend<D1>& result, const mpfi_float_backend<D2>& o)
534 mpfi_sub(result.data(), result.data(), o.data());
536 template <unsigned D1, unsigned D2>
537 inline void eval_multiply(mpfi_float_backend<D1>& result, const mpfi_float_backend<D2>& o)
539 if ((void*)&result == (void*)&o)
540 mpfi_sqr(result.data(), o.data());
542 mpfi_mul(result.data(), result.data(), o.data());
544 template <unsigned D1, unsigned D2>
545 inline void eval_divide(mpfi_float_backend<D1>& result, const mpfi_float_backend<D2>& o)
547 mpfi_div(result.data(), result.data(), o.data());
549 template <unsigned digits10>
550 inline void eval_add(mpfi_float_backend<digits10>& result, unsigned long i)
552 mpfi_add_ui(result.data(), result.data(), i);
554 template <unsigned digits10>
555 inline void eval_subtract(mpfi_float_backend<digits10>& result, unsigned long i)
557 mpfi_sub_ui(result.data(), result.data(), i);
559 template <unsigned digits10>
560 inline void eval_multiply(mpfi_float_backend<digits10>& result, unsigned long i)
562 mpfi_mul_ui(result.data(), result.data(), i);
564 template <unsigned digits10>
565 inline void eval_divide(mpfi_float_backend<digits10>& result, unsigned long i)
567 mpfi_div_ui(result.data(), result.data(), i);
569 template <unsigned digits10>
570 inline void eval_add(mpfi_float_backend<digits10>& result, long i)
573 mpfi_add_ui(result.data(), result.data(), i);
575 mpfi_sub_ui(result.data(), result.data(), boost::multiprecision::detail::unsigned_abs(i));
577 template <unsigned digits10>
578 inline void eval_subtract(mpfi_float_backend<digits10>& result, long i)
581 mpfi_sub_ui(result.data(), result.data(), i);
583 mpfi_add_ui(result.data(), result.data(), boost::multiprecision::detail::unsigned_abs(i));
585 template <unsigned digits10>
586 inline void eval_multiply(mpfi_float_backend<digits10>& result, long i)
588 mpfi_mul_ui(result.data(), result.data(), boost::multiprecision::detail::unsigned_abs(i));
590 mpfi_neg(result.data(), result.data());
592 template <unsigned digits10>
593 inline void eval_divide(mpfi_float_backend<digits10>& result, long i)
595 mpfi_div_ui(result.data(), result.data(), boost::multiprecision::detail::unsigned_abs(i));
597 mpfi_neg(result.data(), result.data());
600 // Specialised 3 arg versions of the basic operators:
602 template <unsigned D1, unsigned D2, unsigned D3>
603 inline void eval_add(mpfi_float_backend<D1>& a, const mpfi_float_backend<D2>& x, const mpfi_float_backend<D3>& y)
605 mpfi_add(a.data(), x.data(), y.data());
607 template <unsigned D1, unsigned D2>
608 inline void eval_add(mpfi_float_backend<D1>& a, const mpfi_float_backend<D2>& x, unsigned long y)
610 mpfi_add_ui(a.data(), x.data(), y);
612 template <unsigned D1, unsigned D2>
613 inline void eval_add(mpfi_float_backend<D1>& a, const mpfi_float_backend<D2>& x, long y)
616 mpfi_sub_ui(a.data(), x.data(), boost::multiprecision::detail::unsigned_abs(y));
618 mpfi_add_ui(a.data(), x.data(), y);
620 template <unsigned D1, unsigned D2>
621 inline void eval_add(mpfi_float_backend<D1>& a, unsigned long x, const mpfi_float_backend<D2>& y)
623 mpfi_add_ui(a.data(), y.data(), x);
625 template <unsigned D1, unsigned D2>
626 inline void eval_add(mpfi_float_backend<D1>& a, long x, const mpfi_float_backend<D2>& y)
630 mpfi_ui_sub(a.data(), boost::multiprecision::detail::unsigned_abs(x), y.data());
631 mpfi_neg(a.data(), a.data());
634 mpfi_add_ui(a.data(), y.data(), x);
636 template <unsigned D1, unsigned D2, unsigned D3>
637 inline void eval_subtract(mpfi_float_backend<D1>& a, const mpfi_float_backend<D2>& x, const mpfi_float_backend<D3>& y)
639 mpfi_sub(a.data(), x.data(), y.data());
641 template <unsigned D1, unsigned D2>
642 inline void eval_subtract(mpfi_float_backend<D1>& a, const mpfi_float_backend<D2>& x, unsigned long y)
644 mpfi_sub_ui(a.data(), x.data(), y);
646 template <unsigned D1, unsigned D2>
647 inline void eval_subtract(mpfi_float_backend<D1>& a, const mpfi_float_backend<D2>& x, long y)
650 mpfi_add_ui(a.data(), x.data(), boost::multiprecision::detail::unsigned_abs(y));
652 mpfi_sub_ui(a.data(), x.data(), y);
654 template <unsigned D1, unsigned D2>
655 inline void eval_subtract(mpfi_float_backend<D1>& a, unsigned long x, const mpfi_float_backend<D2>& y)
657 mpfi_ui_sub(a.data(), x, y.data());
659 template <unsigned D1, unsigned D2>
660 inline void eval_subtract(mpfi_float_backend<D1>& a, long x, const mpfi_float_backend<D2>& y)
664 mpfi_add_ui(a.data(), y.data(), boost::multiprecision::detail::unsigned_abs(x));
665 mpfi_neg(a.data(), a.data());
668 mpfi_ui_sub(a.data(), x, y.data());
671 template <unsigned D1, unsigned D2, unsigned D3>
672 inline void eval_multiply(mpfi_float_backend<D1>& a, const mpfi_float_backend<D2>& x, const mpfi_float_backend<D3>& y)
674 if ((void*)&x == (void*)&y)
675 mpfi_sqr(a.data(), x.data());
677 mpfi_mul(a.data(), x.data(), y.data());
679 template <unsigned D1, unsigned D2>
680 inline void eval_multiply(mpfi_float_backend<D1>& a, const mpfi_float_backend<D2>& x, unsigned long y)
682 mpfi_mul_ui(a.data(), x.data(), y);
684 template <unsigned D1, unsigned D2>
685 inline void eval_multiply(mpfi_float_backend<D1>& a, const mpfi_float_backend<D2>& x, long y)
689 mpfi_mul_ui(a.data(), x.data(), boost::multiprecision::detail::unsigned_abs(y));
693 mpfi_mul_ui(a.data(), x.data(), y);
695 template <unsigned D1, unsigned D2>
696 inline void eval_multiply(mpfi_float_backend<D1>& a, unsigned long x, const mpfi_float_backend<D2>& y)
698 mpfi_mul_ui(a.data(), y.data(), x);
700 template <unsigned D1, unsigned D2>
701 inline void eval_multiply(mpfi_float_backend<D1>& a, long x, const mpfi_float_backend<D2>& y)
705 mpfi_mul_ui(a.data(), y.data(), boost::multiprecision::detail::unsigned_abs(x));
706 mpfi_neg(a.data(), a.data());
709 mpfi_mul_ui(a.data(), y.data(), x);
712 template <unsigned D1, unsigned D2, unsigned D3>
713 inline void eval_divide(mpfi_float_backend<D1>& a, const mpfi_float_backend<D2>& x, const mpfi_float_backend<D3>& y)
715 mpfi_div(a.data(), x.data(), y.data());
717 template <unsigned D1, unsigned D2>
718 inline void eval_divide(mpfi_float_backend<D1>& a, const mpfi_float_backend<D2>& x, unsigned long y)
720 mpfi_div_ui(a.data(), x.data(), y);
722 template <unsigned D1, unsigned D2>
723 inline void eval_divide(mpfi_float_backend<D1>& a, const mpfi_float_backend<D2>& x, long y)
727 mpfi_div_ui(a.data(), x.data(), boost::multiprecision::detail::unsigned_abs(y));
731 mpfi_div_ui(a.data(), x.data(), y);
733 template <unsigned D1, unsigned D2>
734 inline void eval_divide(mpfi_float_backend<D1>& a, unsigned long x, const mpfi_float_backend<D2>& y)
736 mpfi_ui_div(a.data(), x, y.data());
738 template <unsigned D1, unsigned D2>
739 inline void eval_divide(mpfi_float_backend<D1>& a, long x, const mpfi_float_backend<D2>& y)
743 mpfi_ui_div(a.data(), boost::multiprecision::detail::unsigned_abs(x), y.data());
744 mpfi_neg(a.data(), a.data());
747 mpfi_ui_div(a.data(), x, y.data());
750 template <unsigned digits10>
751 inline bool eval_is_zero(const mpfi_float_backend<digits10>& val) BOOST_NOEXCEPT
753 return 0 != mpfi_is_zero(val.data());
755 template <unsigned digits10>
756 inline int eval_get_sign(const mpfi_float_backend<digits10>& val)
758 return detail::mpfi_sgn(val.data());
761 template <unsigned digits10>
762 inline void eval_convert_to(unsigned long* result, const mpfi_float_backend<digits10>& val)
764 mpfr_float_backend<digits10> t;
765 mpfi_mid(t.data(), val.data());
766 eval_convert_to(result, t);
768 template <unsigned digits10>
769 inline void eval_convert_to(long* result, const mpfi_float_backend<digits10>& val)
771 mpfr_float_backend<digits10> t;
772 mpfi_mid(t.data(), val.data());
773 eval_convert_to(result, t);
775 #ifdef _MPFR_H_HAVE_INTMAX_T
776 template <unsigned digits10>
777 inline void eval_convert_to(boost::ulong_long_type* result, const mpfi_float_backend<digits10>& val)
779 mpfr_float_backend<digits10> t;
780 mpfi_mid(t.data(), val.data());
781 eval_convert_to(result, t);
783 template <unsigned digits10>
784 inline void eval_convert_to(boost::long_long_type* result, const mpfi_float_backend<digits10>& val)
786 mpfr_float_backend<digits10> t;
787 mpfi_mid(t.data(), val.data());
788 eval_convert_to(result, t);
791 template <unsigned digits10>
792 inline void eval_convert_to(double* result, const mpfi_float_backend<digits10>& val) BOOST_NOEXCEPT
794 *result = mpfi_get_d(val.data());
796 template <unsigned digits10>
797 inline void eval_convert_to(long double* result, const mpfi_float_backend<digits10>& val) BOOST_NOEXCEPT
799 mpfr_float_backend<digits10> t;
800 mpfi_mid(t.data(), val.data());
801 eval_convert_to(result, t);
804 template <unsigned D1, unsigned D2, mpfr_allocation_type AllocationType>
805 inline void assign_components(mpfi_float_backend<D1>& result, const mpfr_float_backend<D2, AllocationType>& a, const mpfr_float_backend<D2, AllocationType>& b)
807 using default_ops::eval_fpclassify;
808 if (eval_fpclassify(a) == (int)FP_NAN)
810 mpfi_set_fr(result.data(), a.data());
812 else if (eval_fpclassify(b) == (int)FP_NAN)
814 mpfi_set_fr(result.data(), b.data());
818 if (a.compare(b) > 0)
820 BOOST_THROW_EXCEPTION(std::runtime_error("Attempt to create interval with invalid range (start is greater than end)."));
822 mpfi_interv_fr(result.data(), a.data(), b.data());
826 template <unsigned Digits10, class V>
827 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
828 assign_components(mpfi_float_backend<Digits10>& result, const V& a, const V& b)
830 number<mpfr_float_backend<Digits10, allocate_dynamic>, et_on> x(a), y(b);
831 assign_components(result, x.backend(), y.backend());
835 // Native non-member operations:
837 template <unsigned Digits10>
838 inline void eval_sqrt(mpfi_float_backend<Digits10>& result, const mpfi_float_backend<Digits10>& val)
840 mpfi_sqrt(result.data(), val.data());
843 template <unsigned Digits10>
844 inline void eval_abs(mpfi_float_backend<Digits10>& result, const mpfi_float_backend<Digits10>& val)
846 mpfi_abs(result.data(), val.data());
849 template <unsigned Digits10>
850 inline void eval_fabs(mpfi_float_backend<Digits10>& result, const mpfi_float_backend<Digits10>& val)
852 mpfi_abs(result.data(), val.data());
854 template <unsigned Digits10>
855 inline void eval_ceil(mpfi_float_backend<Digits10>& result, const mpfi_float_backend<Digits10>& val)
857 mpfr_float_backend<Digits10> a, b;
858 mpfr_set(a.data(), val.left_data(), GMP_RNDN);
859 mpfr_set(b.data(), val.right_data(), GMP_RNDN);
862 if (a.compare(b) != 0)
864 BOOST_THROW_EXCEPTION(interval_error("Attempt to take the ceil of a value that straddles an integer boundary."));
866 mpfi_set_fr(result.data(), a.data());
868 template <unsigned Digits10>
869 inline void eval_floor(mpfi_float_backend<Digits10>& result, const mpfi_float_backend<Digits10>& val)
871 mpfr_float_backend<Digits10> a, b;
872 mpfr_set(a.data(), val.left_data(), GMP_RNDN);
873 mpfr_set(b.data(), val.right_data(), GMP_RNDN);
876 if (a.compare(b) != 0)
878 BOOST_THROW_EXCEPTION(interval_error("Attempt to take the floor of a value that straddles an integer boundary."));
880 mpfi_set_fr(result.data(), a.data());
882 template <unsigned Digits10>
883 inline void eval_ldexp(mpfi_float_backend<Digits10>& result, const mpfi_float_backend<Digits10>& val, long e)
886 mpfi_mul_2exp(result.data(), val.data(), e);
888 mpfi_div_2exp(result.data(), val.data(), -e);
892 template <unsigned Digits10>
893 inline void eval_frexp(mpfi_float_backend<Digits10>& result, const mpfi_float_backend<Digits10>& val, int* e)
895 mpfr_float_backend<Digits10> t, rt;
896 mpfi_mid(t.data(), val.data());
897 eval_frexp(rt, t, e);
898 eval_ldexp(result, val, -*e);
900 template <unsigned Digits10>
901 inline void eval_frexp(mpfi_float_backend<Digits10>& result, const mpfi_float_backend<Digits10>& val, long* e)
903 mpfr_float_backend<Digits10> t, rt;
904 mpfi_mid(t.data(), val.data());
905 eval_frexp(rt, t, e);
906 eval_ldexp(result, val, -*e);
909 template <unsigned Digits10>
910 inline int eval_fpclassify(const mpfi_float_backend<Digits10>& val) BOOST_NOEXCEPT
912 return mpfi_inf_p(val.data()) ? FP_INFINITE : mpfi_nan_p(val.data()) ? FP_NAN : mpfi_is_zero(val.data()) ? FP_ZERO : FP_NORMAL;
915 template <unsigned Digits10>
916 inline void eval_pow(mpfi_float_backend<Digits10>& result, const mpfi_float_backend<Digits10>& b, const mpfi_float_backend<Digits10>& e)
918 typedef typename boost::multiprecision::detail::canonical<unsigned, mpfi_float_backend<Digits10> >::type ui_type;
919 using default_ops::eval_get_sign;
920 int s = eval_get_sign(b);
923 if (eval_get_sign(e) == 0)
935 if (eval_get_sign(e) < 0)
937 mpfi_float_backend<Digits10> t1, t2;
942 eval_divide(result, t1, t2);
945 typename boost::multiprecision::detail::canonical<boost::uintmax_t, mpfi_float_backend<Digits10> >::type an;
946 #ifndef BOOST_NO_EXCEPTIONS
950 using default_ops::eval_convert_to;
951 eval_convert_to(&an, e);
952 if (e.compare(an) == 0)
954 mpfi_float_backend<Digits10> pb(b);
956 eval_pow(result, pb, e);
961 #ifndef BOOST_NO_EXCEPTIONS
963 catch (const std::exception&)
965 // conversion failed, just fall through, value is not an integer.
968 result = std::numeric_limits<number<mpfi_float_backend<Digits10>, et_on> >::quiet_NaN().backend();
971 mpfi_log(result.data(), b.data());
972 mpfi_mul(result.data(), result.data(), e.data());
973 mpfi_exp(result.data(), result.data());
976 template <unsigned Digits10>
977 inline void eval_exp(mpfi_float_backend<Digits10>& result, const mpfi_float_backend<Digits10>& arg)
979 mpfi_exp(result.data(), arg.data());
982 template <unsigned Digits10>
983 inline void eval_exp2(mpfi_float_backend<Digits10>& result, const mpfi_float_backend<Digits10>& arg)
985 mpfi_exp2(result.data(), arg.data());
988 template <unsigned Digits10>
989 inline void eval_log(mpfi_float_backend<Digits10>& result, const mpfi_float_backend<Digits10>& arg)
991 mpfi_log(result.data(), arg.data());
994 template <unsigned Digits10>
995 inline void eval_log10(mpfi_float_backend<Digits10>& result, const mpfi_float_backend<Digits10>& arg)
997 mpfi_log10(result.data(), arg.data());
1000 template <unsigned Digits10>
1001 inline void eval_sin(mpfi_float_backend<Digits10>& result, const mpfi_float_backend<Digits10>& arg)
1003 mpfi_sin(result.data(), arg.data());
1006 template <unsigned Digits10>
1007 inline void eval_cos(mpfi_float_backend<Digits10>& result, const mpfi_float_backend<Digits10>& arg)
1009 mpfi_cos(result.data(), arg.data());
1012 template <unsigned Digits10>
1013 inline void eval_tan(mpfi_float_backend<Digits10>& result, const mpfi_float_backend<Digits10>& arg)
1015 mpfi_tan(result.data(), arg.data());
1018 template <unsigned Digits10>
1019 inline void eval_asin(mpfi_float_backend<Digits10>& result, const mpfi_float_backend<Digits10>& arg)
1021 mpfi_asin(result.data(), arg.data());
1024 template <unsigned Digits10>
1025 inline void eval_acos(mpfi_float_backend<Digits10>& result, const mpfi_float_backend<Digits10>& arg)
1027 mpfi_acos(result.data(), arg.data());
1030 template <unsigned Digits10>
1031 inline void eval_atan(mpfi_float_backend<Digits10>& result, const mpfi_float_backend<Digits10>& arg)
1033 mpfi_atan(result.data(), arg.data());
1036 template <unsigned Digits10>
1037 inline void eval_atan2(mpfi_float_backend<Digits10>& result, const mpfi_float_backend<Digits10>& arg1, const mpfi_float_backend<Digits10>& arg2)
1039 mpfi_atan2(result.data(), arg1.data(), arg2.data());
1042 template <unsigned Digits10>
1043 inline void eval_sinh(mpfi_float_backend<Digits10>& result, const mpfi_float_backend<Digits10>& arg)
1045 mpfi_sinh(result.data(), arg.data());
1048 template <unsigned Digits10>
1049 inline void eval_cosh(mpfi_float_backend<Digits10>& result, const mpfi_float_backend<Digits10>& arg)
1051 mpfi_cosh(result.data(), arg.data());
1054 template <unsigned Digits10>
1055 inline void eval_tanh(mpfi_float_backend<Digits10>& result, const mpfi_float_backend<Digits10>& arg)
1057 mpfi_tanh(result.data(), arg.data());
1060 template <unsigned Digits10>
1061 inline void eval_log2(mpfi_float_backend<Digits10>& result, const mpfi_float_backend<Digits10>& arg)
1063 mpfi_log2(result.data(), arg.data());
1066 template <unsigned Digits10>
1067 inline std::size_t hash_value(const mpfi_float_backend<Digits10>& val)
1069 std::size_t result = 0;
1070 std::size_t len = val.left_data()[0]._mpfr_prec / mp_bits_per_limb;
1071 if (val.left_data()[0]._mpfr_prec % mp_bits_per_limb)
1073 for (std::size_t i = 0; i < len; ++i)
1074 boost::hash_combine(result, val.left_data()[0]._mpfr_d[i]);
1075 boost::hash_combine(result, val.left_data()[0]._mpfr_exp);
1076 boost::hash_combine(result, val.left_data()[0]._mpfr_sign);
1078 len = val.right_data()[0]._mpfr_prec / mp_bits_per_limb;
1079 if (val.right_data()[0]._mpfr_prec % mp_bits_per_limb)
1081 for (std::size_t i = 0; i < len; ++i)
1082 boost::hash_combine(result, val.right_data()[0]._mpfr_d[i]);
1083 boost::hash_combine(result, val.right_data()[0]._mpfr_exp);
1084 boost::hash_combine(result, val.right_data()[0]._mpfr_sign);
1088 template <class To, unsigned D>
1089 void generic_interconvert(To& to, const mpfi_float_backend<D>& from, const mpl::int_<number_kind_integer>& to_type, const mpl::int_<number_kind_floating_point>& from_type)
1091 using boost::multiprecision::detail::generic_interconvert;
1092 mpfr_float_backend<D> t;
1093 mpfi_mid(t.data(), from.data());
1094 generic_interconvert(to, t, to_type, from_type);
1097 template <class To, unsigned D>
1098 void generic_interconvert(To& to, const mpfi_float_backend<D>& from, const mpl::int_<number_kind_rational>& to_type, const mpl::int_<number_kind_floating_point>& from_type)
1100 using boost::multiprecision::detail::generic_interconvert;
1101 mpfr_float_backend<D> t;
1102 mpfi_mid(t.data(), from.data());
1103 generic_interconvert(to, t, to_type, from_type);
1106 template <class To, unsigned D>
1107 void generic_interconvert(To& to, const mpfi_float_backend<D>& from, const mpl::int_<number_kind_floating_point>& to_type, const mpl::int_<number_kind_floating_point>& from_type)
1109 using boost::multiprecision::detail::generic_interconvert;
1110 mpfr_float_backend<D> t;
1111 mpfi_mid(t.data(), from.data());
1112 generic_interconvert(to, t, to_type, from_type);
1115 } // namespace backends
1117 #ifdef BOOST_NO_SFINAE_EXPR
1121 template <unsigned D1, unsigned D2>
1122 struct is_explicitly_convertible<backends::mpfi_float_backend<D1>, backends::mpfi_float_backend<D2> > : public mpl::true_
1125 } // namespace detail
1130 struct is_variable_precision<backends::mpfi_float_backend<0> > : public true_type
1132 } // namespace detail
1135 struct number_category<detail::canonical<mpfi_t, backends::mpfi_float_backend<0> >::type> : public mpl::int_<number_kind_floating_point>
1137 template <unsigned Digits10>
1138 struct is_interval_number<backends::mpfi_float_backend<Digits10> > : public mpl::true_
1141 using boost::multiprecision::backends::mpfi_float_backend;
1143 typedef number<mpfi_float_backend<50> > mpfi_float_50;
1144 typedef number<mpfi_float_backend<100> > mpfi_float_100;
1145 typedef number<mpfi_float_backend<500> > mpfi_float_500;
1146 typedef number<mpfi_float_backend<1000> > mpfi_float_1000;
1147 typedef number<mpfi_float_backend<0> > mpfi_float;
1150 // Special interval specific functions:
1152 template <unsigned Digits10, expression_template_option ExpressionTemplates>
1153 inline number<mpfr_float_backend<Digits10>, ExpressionTemplates> lower(const number<mpfi_float_backend<Digits10>, ExpressionTemplates>& val)
1155 boost::multiprecision::detail::scoped_default_precision<number<mpfr_float_backend<Digits10>, ExpressionTemplates> > precision_guard(val);
1156 number<mpfr_float_backend<Digits10> > result;
1157 mpfr_set(result.backend().data(), val.backend().left_data(), GMP_RNDN);
1161 template <unsigned Digits10, expression_template_option ExpressionTemplates>
1162 inline number<mpfr_float_backend<Digits10>, ExpressionTemplates> upper(const number<mpfi_float_backend<Digits10>, ExpressionTemplates>& val)
1164 boost::multiprecision::detail::scoped_default_precision<number<mpfr_float_backend<Digits10>, ExpressionTemplates> > precision_guard(val);
1165 number<mpfr_float_backend<Digits10> > result;
1166 mpfr_set(result.backend().data(), val.backend().right_data(), GMP_RNDN);
1170 template <unsigned Digits10, expression_template_option ExpressionTemplates>
1171 inline number<mpfr_float_backend<Digits10>, ExpressionTemplates> median(const number<mpfi_float_backend<Digits10>, ExpressionTemplates>& val)
1173 boost::multiprecision::detail::scoped_default_precision<number<mpfr_float_backend<Digits10>, ExpressionTemplates> > precision_guard(val);
1174 number<mpfr_float_backend<Digits10> > result;
1175 mpfi_mid(result.backend().data(), val.backend().data());
1179 template <unsigned Digits10, expression_template_option ExpressionTemplates>
1180 inline number<mpfr_float_backend<Digits10>, ExpressionTemplates> width(const number<mpfi_float_backend<Digits10>, ExpressionTemplates>& val)
1182 boost::multiprecision::detail::scoped_default_precision<number<mpfr_float_backend<Digits10>, ExpressionTemplates> > precision_guard(val);
1183 number<mpfr_float_backend<Digits10> > result;
1184 mpfi_diam_abs(result.backend().data(), val.backend().data());
1188 template <unsigned Digits10, expression_template_option ExpressionTemplates>
1189 inline number<mpfi_float_backend<Digits10>, ExpressionTemplates> intersect(const number<mpfi_float_backend<Digits10>, ExpressionTemplates>& a, const number<mpfi_float_backend<Digits10>, ExpressionTemplates>& b)
1191 boost::multiprecision::detail::scoped_default_precision<number<mpfi_float_backend<Digits10>, ExpressionTemplates> > precision_guard(a, b);
1192 number<mpfi_float_backend<Digits10>, ExpressionTemplates> result;
1193 mpfi_intersect(result.backend().data(), a.backend().data(), b.backend().data());
1197 template <unsigned Digits10, expression_template_option ExpressionTemplates>
1198 inline number<mpfi_float_backend<Digits10>, ExpressionTemplates> hull(const number<mpfi_float_backend<Digits10>, ExpressionTemplates>& a, const number<mpfi_float_backend<Digits10>, ExpressionTemplates>& b)
1200 boost::multiprecision::detail::scoped_default_precision<number<mpfi_float_backend<Digits10>, ExpressionTemplates> > precision_guard(a, b);
1201 number<mpfi_float_backend<Digits10>, ExpressionTemplates> result;
1202 mpfi_union(result.backend().data(), a.backend().data(), b.backend().data());
1206 template <unsigned Digits10, expression_template_option ExpressionTemplates>
1207 inline bool overlap(const number<mpfi_float_backend<Digits10>, ExpressionTemplates>& a, const number<mpfi_float_backend<Digits10>, ExpressionTemplates>& b)
1209 return (lower(a) <= lower(b) && lower(b) <= upper(a)) ||
1210 (lower(b) <= lower(a) && lower(a) <= upper(b));
1213 template <unsigned Digits10, expression_template_option ExpressionTemplates1, expression_template_option ExpressionTemplates2>
1214 inline bool in(const number<mpfr_float_backend<Digits10>, ExpressionTemplates1>& a, const number<mpfi_float_backend<Digits10>, ExpressionTemplates2>& b)
1216 return mpfi_is_inside_fr(a.backend().data(), b.backend().data()) != 0;
1219 template <unsigned Digits10, expression_template_option ExpressionTemplates>
1220 inline bool zero_in(const number<mpfi_float_backend<Digits10>, ExpressionTemplates>& a)
1222 return mpfi_has_zero(a.backend().data()) != 0;
1225 template <unsigned Digits10, expression_template_option ExpressionTemplates>
1226 inline bool subset(const number<mpfi_float_backend<Digits10>, ExpressionTemplates>& a, const number<mpfi_float_backend<Digits10>, ExpressionTemplates>& b)
1228 return mpfi_is_inside(a.backend().data(), b.backend().data()) != 0;
1231 template <unsigned Digits10, expression_template_option ExpressionTemplates>
1232 inline bool proper_subset(const number<mpfi_float_backend<Digits10>, ExpressionTemplates>& a, const number<mpfi_float_backend<Digits10>, ExpressionTemplates>& b)
1234 return mpfi_is_strictly_inside(a.backend().data(), b.backend().data()) != 0;
1237 template <unsigned Digits10, expression_template_option ExpressionTemplates>
1238 inline bool empty(const number<mpfi_float_backend<Digits10>, ExpressionTemplates>& a)
1240 return mpfi_is_empty(a.backend().data()) != 0;
1243 template <unsigned Digits10, expression_template_option ExpressionTemplates>
1244 inline bool singleton(const number<mpfi_float_backend<Digits10>, ExpressionTemplates>& a)
1246 return mpfr_cmp(a.backend().left_data(), a.backend().right_data()) == 0;
1249 template <unsigned Digits10, expression_template_option ExpressionTemplates>
1250 struct component_type<number<mpfi_float_backend<Digits10>, ExpressionTemplates> >
1252 typedef number<mpfr_float_backend<Digits10>, ExpressionTemplates> type;
1256 // Overloaded special functions which call native mpfr routines:
1258 template <unsigned Digits10, expression_template_option ExpressionTemplates>
1259 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)
1261 boost::multiprecision::detail::scoped_default_precision<number<mpfi_float_backend<Digits10>, ExpressionTemplates> > precision_guard(arg);
1263 boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> result;
1264 mpfi_asinh(result.backend().data(), arg.backend().data());
1267 template <unsigned Digits10, expression_template_option ExpressionTemplates>
1268 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)
1270 boost::multiprecision::detail::scoped_default_precision<number<mpfi_float_backend<Digits10>, ExpressionTemplates> > precision_guard(arg);
1272 boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> result;
1273 mpfi_acosh(result.backend().data(), arg.backend().data());
1276 template <unsigned Digits10, expression_template_option ExpressionTemplates>
1277 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)
1279 boost::multiprecision::detail::scoped_default_precision<number<mpfi_float_backend<Digits10>, ExpressionTemplates> > precision_guard(arg);
1281 boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> result;
1282 mpfi_atanh(result.backend().data(), arg.backend().data());
1285 template <unsigned Digits10, expression_template_option ExpressionTemplates>
1286 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)
1288 boost::multiprecision::detail::scoped_default_precision<number<mpfi_float_backend<Digits10>, ExpressionTemplates> > precision_guard(arg);
1290 boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> result;
1291 mpfi_cbrt(result.backend().data(), arg.backend().data());
1294 template <unsigned Digits10, expression_template_option ExpressionTemplates>
1295 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)
1297 boost::multiprecision::detail::scoped_default_precision<number<mpfi_float_backend<Digits10>, ExpressionTemplates> > precision_guard(arg);
1299 boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> result;
1300 mpfi_expm1(result.backend().data(), arg.backend().data());
1303 template <unsigned Digits10, expression_template_option ExpressionTemplates>
1304 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)
1306 boost::multiprecision::detail::scoped_default_precision<number<mpfi_float_backend<Digits10>, ExpressionTemplates> > precision_guard(arg);
1308 boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> result;
1309 mpfi_log1p(result.backend().data(), arg.backend().data());
1313 } // namespace multiprecision
1320 inline int digits<boost::multiprecision::mpfi_float>()
1321 #ifdef BOOST_MATH_NOEXCEPT
1325 return multiprecision::detail::digits10_2_2(boost::multiprecision::mpfi_float::default_precision());
1328 inline int digits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<0>, boost::multiprecision::et_off> >()
1329 #ifdef BOOST_MATH_NOEXCEPT
1333 return multiprecision::detail::digits10_2_2(boost::multiprecision::mpfi_float::default_precision());
1337 inline boost::multiprecision::mpfi_float
1338 max_value<boost::multiprecision::mpfi_float>()
1340 boost::multiprecision::mpfi_float result(0.5);
1341 mpfi_mul_2exp(result.backend().data(), result.backend().data(), mpfr_get_emax());
1342 //BOOST_ASSERT(mpfi_number_p(result.backend().data()));
1347 inline boost::multiprecision::mpfi_float
1348 min_value<boost::multiprecision::mpfi_float>()
1350 boost::multiprecision::mpfi_float result(0.5);
1351 mpfi_div_2exp(result.backend().data(), result.backend().data(), -mpfr_get_emin());
1352 //BOOST_ASSERT(mpfi_number_p(result.backend().data()));
1357 inline boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<0>, boost::multiprecision::et_off>
1358 max_value<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<0>, boost::multiprecision::et_off> >()
1360 boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<0>, boost::multiprecision::et_off> result(0.5);
1361 mpfi_mul_2exp(result.backend().data(), result.backend().data(), mpfr_get_emax());
1362 //BOOST_ASSERT(mpfi_number_p(result.backend().data()));
1367 inline boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<0>, boost::multiprecision::et_off>
1368 min_value<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<0>, boost::multiprecision::et_off> >()
1370 boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<0>, boost::multiprecision::et_off> result(0.5);
1371 mpfi_div_2exp(result.backend().data(), result.backend().data(), -mpfr_get_emin());
1372 //BOOST_ASSERT(mpfi_number_p(result.backend().data()));
1376 // mpfi gets used with logged_adaptor fairly often, so specialize for that use case as well:
1377 typedef boost::multiprecision::number<boost::multiprecision::backends::logged_adaptor<boost::multiprecision::mpfi_float::backend_type>, boost::multiprecision::et_on> logged_type1;
1378 typedef boost::multiprecision::number<boost::multiprecision::backends::logged_adaptor<boost::multiprecision::mpfi_float::backend_type>, boost::multiprecision::et_off> logged_type2;
1381 inline int digits<logged_type1>()
1382 #ifdef BOOST_MATH_NOEXCEPT
1386 return multiprecision::detail::digits10_2_2(logged_type1::default_precision());
1389 inline int digits<logged_type2>()
1390 #ifdef BOOST_MATH_NOEXCEPT
1394 return multiprecision::detail::digits10_2_2(logged_type1::default_precision());
1399 max_value<logged_type1>()
1401 logged_type1 result(0.5);
1402 mpfi_mul_2exp(result.backend().value().data(), result.backend().value().data(), mpfr_get_emax());
1403 //BOOST_ASSERT(mpfi_number_p(result.backend().data()));
1409 min_value<logged_type1>()
1411 logged_type1 result(0.5);
1412 mpfi_div_2exp(result.backend().value().data(), result.backend().value().data(), -mpfr_get_emin());
1413 //BOOST_ASSERT(mpfi_number_p(result.backend().data()));
1419 max_value<logged_type2>()
1421 logged_type2 result(0.5);
1422 mpfi_mul_2exp(result.backend().value().data(), result.backend().value().data(), mpfr_get_emax());
1423 //BOOST_ASSERT(mpfi_number_p(result.backend().data()));
1429 min_value<logged_type2>()
1431 logged_type2 result(0.5);
1432 mpfi_div_2exp(result.backend().value().data(), result.backend().value().data(), -mpfr_get_emin());
1433 //BOOST_ASSERT(mpfi_number_p(result.backend().data()));
1436 } // namespace tools
1438 namespace constants { namespace detail {
1443 struct constant_ln_two;
1445 struct constant_euler;
1447 struct constant_catalan;
1450 // Initializer: ensure all our constants are initialized prior to the first call of main:
1453 struct mpfi_initializer
1459 boost::math::constants::pi<T>();
1460 boost::math::constants::ln_two<T>();
1461 boost::math::constants::euler<T>();
1462 boost::math::constants::catalan<T>();
1464 void force_instantiate() const {}
1466 static const init initializer;
1467 static void force_instantiate()
1469 initializer.force_instantiate();
1474 const typename mpfi_initializer<T>::init mpfi_initializer<T>::initializer;
1476 template <unsigned Digits10, boost::multiprecision::expression_template_option ExpressionTemplates>
1477 struct constant_pi<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> >
1479 typedef boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> result_type;
1481 static inline const result_type& get(const mpl::int_<N>&)
1483 mpfi_initializer<result_type>::force_instantiate();
1484 static result_type result;
1485 static bool has_init = false;
1489 mpfi_const_pi(result.backend().data());
1493 static inline result_type get(const mpl::int_<0>&)
1496 mpfi_const_pi(result.backend().data());
1500 template <unsigned Digits10, boost::multiprecision::expression_template_option ExpressionTemplates>
1501 struct constant_ln_two<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> >
1503 typedef boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> result_type;
1505 static inline const result_type& get(const mpl::int_<N>&)
1507 mpfi_initializer<result_type>::force_instantiate();
1508 static result_type result;
1509 static bool has_init = false;
1513 mpfi_const_log2(result.backend().data());
1517 static inline result_type get(const mpl::int_<0>&)
1520 mpfi_const_log2(result.backend().data());
1524 template <unsigned Digits10, boost::multiprecision::expression_template_option ExpressionTemplates>
1525 struct constant_euler<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> >
1527 typedef boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> result_type;
1529 static inline result_type const& get(const mpl::int_<N>&)
1531 mpfi_initializer<result_type>::force_instantiate();
1532 static result_type result;
1533 static bool has_init = false;
1537 mpfi_const_euler(result.backend().data());
1541 static inline result_type get(const mpl::int_<0>&)
1544 mpfi_const_euler(result.backend().data());
1548 template <unsigned Digits10, boost::multiprecision::expression_template_option ExpressionTemplates>
1549 struct constant_catalan<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> >
1551 typedef boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> result_type;
1553 static inline result_type const& get(const mpl::int_<N>&)
1555 mpfi_initializer<result_type>::force_instantiate();
1556 static result_type result;
1557 static bool has_init = false;
1561 mpfi_const_catalan(result.backend().data());
1565 static inline result_type get(const mpl::int_<0>&)
1568 mpfi_const_catalan(result.backend().data());
1573 }} // namespace constants::detail
1576 } // namespace boost
1581 // numeric_limits [partial] specializations for the types declared in this header:
1583 template <unsigned Digits10, boost::multiprecision::expression_template_option ExpressionTemplates>
1584 class numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> >
1586 typedef boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> number_type;
1589 BOOST_STATIC_CONSTEXPR bool is_specialized = true;
1590 static number_type(min)()
1592 initializer.do_nothing();
1593 static std::pair<bool, number_type> value;
1598 mpfi_div_2exp(value.second.backend().data(), value.second.backend().data(), -mpfr_get_emin());
1600 return value.second;
1602 static number_type(max)()
1604 initializer.do_nothing();
1605 static std::pair<bool, number_type> value;
1610 mpfi_mul_2exp(value.second.backend().data(), value.second.backend().data(), mpfr_get_emax());
1612 return value.second;
1614 BOOST_STATIC_CONSTEXPR number_type lowest()
1618 BOOST_STATIC_CONSTEXPR int digits = static_cast<int>((Digits10 * 1000L) / 301L + ((Digits10 * 1000L) % 301 ? 2 : 1));
1619 BOOST_STATIC_CONSTEXPR int digits10 = Digits10;
1620 // Is this really correct???
1621 BOOST_STATIC_CONSTEXPR int max_digits10 = Digits10 + 3;
1622 BOOST_STATIC_CONSTEXPR bool is_signed = true;
1623 BOOST_STATIC_CONSTEXPR bool is_integer = false;
1624 BOOST_STATIC_CONSTEXPR bool is_exact = false;
1625 BOOST_STATIC_CONSTEXPR int radix = 2;
1626 static number_type epsilon()
1628 initializer.do_nothing();
1629 static std::pair<bool, number_type> value;
1634 mpfi_div_2exp(value.second.backend().data(), value.second.backend().data(), std::numeric_limits<number_type>::digits - 1);
1636 return value.second;
1638 // What value should this be????
1639 static number_type round_error()
1641 // returns epsilon/2
1642 initializer.do_nothing();
1643 static std::pair<bool, number_type> value;
1648 mpfi_div_2exp(value.second.backend().data(), value.second.backend().data(), 1);
1650 return value.second;
1652 BOOST_STATIC_CONSTEXPR long min_exponent = MPFR_EMIN_DEFAULT;
1653 BOOST_STATIC_CONSTEXPR long min_exponent10 = (MPFR_EMIN_DEFAULT / 1000) * 301L;
1654 BOOST_STATIC_CONSTEXPR long max_exponent = MPFR_EMAX_DEFAULT;
1655 BOOST_STATIC_CONSTEXPR long max_exponent10 = (MPFR_EMAX_DEFAULT / 1000) * 301L;
1656 BOOST_STATIC_CONSTEXPR bool has_infinity = true;
1657 BOOST_STATIC_CONSTEXPR bool has_quiet_NaN = true;
1658 BOOST_STATIC_CONSTEXPR bool has_signaling_NaN = false;
1659 BOOST_STATIC_CONSTEXPR float_denorm_style has_denorm = denorm_absent;
1660 BOOST_STATIC_CONSTEXPR bool has_denorm_loss = false;
1661 static number_type infinity()
1663 initializer.do_nothing();
1664 static std::pair<bool, number_type> value;
1667 boost::multiprecision::mpfr_float_backend<Digits10> t;
1668 mpfr_set_inf(t.data(), 1);
1670 mpfi_set_fr(value.second.backend().data(), t.data());
1672 return value.second;
1674 static number_type quiet_NaN()
1676 initializer.do_nothing();
1677 static std::pair<bool, number_type> value;
1680 boost::multiprecision::mpfr_float_backend<Digits10> t;
1681 mpfr_set_nan(t.data());
1683 mpfi_set_fr(value.second.backend().data(), t.data());
1685 return value.second;
1687 BOOST_STATIC_CONSTEXPR number_type signaling_NaN()
1689 return number_type(0);
1691 BOOST_STATIC_CONSTEXPR number_type denorm_min() { return number_type(0); }
1692 BOOST_STATIC_CONSTEXPR bool is_iec559 = false;
1693 BOOST_STATIC_CONSTEXPR bool is_bounded = true;
1694 BOOST_STATIC_CONSTEXPR bool is_modulo = false;
1695 BOOST_STATIC_CONSTEXPR bool traps = true;
1696 BOOST_STATIC_CONSTEXPR bool tinyness_before = false;
1697 BOOST_STATIC_CONSTEXPR float_round_style round_style = round_to_nearest;
1700 struct data_initializer
1704 std::numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<digits10> > >::epsilon();
1705 std::numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<digits10> > >::round_error();
1706 (std::numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<digits10> > >::min)();
1707 (std::numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<digits10> > >::max)();
1708 std::numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<digits10> > >::infinity();
1709 std::numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<digits10> > >::quiet_NaN();
1711 void do_nothing() const {}
1713 static const data_initializer initializer;
1716 template <unsigned Digits10, boost::multiprecision::expression_template_option ExpressionTemplates>
1717 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;
1719 #ifndef BOOST_NO_INCLASS_MEMBER_INITIALIZATION
1721 template <unsigned Digits10, boost::multiprecision::expression_template_option ExpressionTemplates>
1722 BOOST_CONSTEXPR_OR_CONST int numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> >::digits;
1723 template <unsigned Digits10, boost::multiprecision::expression_template_option ExpressionTemplates>
1724 BOOST_CONSTEXPR_OR_CONST int numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> >::digits10;
1725 template <unsigned Digits10, boost::multiprecision::expression_template_option ExpressionTemplates>
1726 BOOST_CONSTEXPR_OR_CONST int numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> >::max_digits10;
1727 template <unsigned Digits10, boost::multiprecision::expression_template_option ExpressionTemplates>
1728 BOOST_CONSTEXPR_OR_CONST bool numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> >::is_signed;
1729 template <unsigned Digits10, boost::multiprecision::expression_template_option ExpressionTemplates>
1730 BOOST_CONSTEXPR_OR_CONST bool numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> >::is_integer;
1731 template <unsigned Digits10, boost::multiprecision::expression_template_option ExpressionTemplates>
1732 BOOST_CONSTEXPR_OR_CONST bool numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> >::is_exact;
1733 template <unsigned Digits10, boost::multiprecision::expression_template_option ExpressionTemplates>
1734 BOOST_CONSTEXPR_OR_CONST int numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> >::radix;
1735 template <unsigned Digits10, boost::multiprecision::expression_template_option ExpressionTemplates>
1736 BOOST_CONSTEXPR_OR_CONST long numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> >::min_exponent;
1737 template <unsigned Digits10, boost::multiprecision::expression_template_option ExpressionTemplates>
1738 BOOST_CONSTEXPR_OR_CONST long numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> >::min_exponent10;
1739 template <unsigned Digits10, boost::multiprecision::expression_template_option ExpressionTemplates>
1740 BOOST_CONSTEXPR_OR_CONST long numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> >::max_exponent;
1741 template <unsigned Digits10, boost::multiprecision::expression_template_option ExpressionTemplates>
1742 BOOST_CONSTEXPR_OR_CONST long numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> >::max_exponent10;
1743 template <unsigned Digits10, boost::multiprecision::expression_template_option ExpressionTemplates>
1744 BOOST_CONSTEXPR_OR_CONST bool numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> >::has_infinity;
1745 template <unsigned Digits10, boost::multiprecision::expression_template_option ExpressionTemplates>
1746 BOOST_CONSTEXPR_OR_CONST bool numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> >::has_quiet_NaN;
1747 template <unsigned Digits10, boost::multiprecision::expression_template_option ExpressionTemplates>
1748 BOOST_CONSTEXPR_OR_CONST bool numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> >::has_signaling_NaN;
1749 template <unsigned Digits10, boost::multiprecision::expression_template_option ExpressionTemplates>
1750 BOOST_CONSTEXPR_OR_CONST float_denorm_style numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> >::has_denorm;
1751 template <unsigned Digits10, boost::multiprecision::expression_template_option ExpressionTemplates>
1752 BOOST_CONSTEXPR_OR_CONST bool numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> >::has_denorm_loss;
1753 template <unsigned Digits10, boost::multiprecision::expression_template_option ExpressionTemplates>
1754 BOOST_CONSTEXPR_OR_CONST bool numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> >::is_iec559;
1755 template <unsigned Digits10, boost::multiprecision::expression_template_option ExpressionTemplates>
1756 BOOST_CONSTEXPR_OR_CONST bool numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> >::is_bounded;
1757 template <unsigned Digits10, boost::multiprecision::expression_template_option ExpressionTemplates>
1758 BOOST_CONSTEXPR_OR_CONST bool numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> >::is_modulo;
1759 template <unsigned Digits10, boost::multiprecision::expression_template_option ExpressionTemplates>
1760 BOOST_CONSTEXPR_OR_CONST bool numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> >::traps;
1761 template <unsigned Digits10, boost::multiprecision::expression_template_option ExpressionTemplates>
1762 BOOST_CONSTEXPR_OR_CONST bool numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> >::tinyness_before;
1763 template <unsigned Digits10, boost::multiprecision::expression_template_option ExpressionTemplates>
1764 BOOST_CONSTEXPR_OR_CONST float_round_style numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<Digits10>, ExpressionTemplates> >::round_style;
1768 template <boost::multiprecision::expression_template_option ExpressionTemplates>
1769 class numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<0>, ExpressionTemplates> >
1771 typedef boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<0>, ExpressionTemplates> number_type;
1774 BOOST_STATIC_CONSTEXPR bool is_specialized = false;
1775 static number_type(min)() { return number_type(0); }
1776 static number_type(max)() { return number_type(0); }
1777 static number_type lowest() { return number_type(0); }
1778 BOOST_STATIC_CONSTEXPR int digits = 0;
1779 BOOST_STATIC_CONSTEXPR int digits10 = 0;
1780 BOOST_STATIC_CONSTEXPR int max_digits10 = 0;
1781 BOOST_STATIC_CONSTEXPR bool is_signed = false;
1782 BOOST_STATIC_CONSTEXPR bool is_integer = false;
1783 BOOST_STATIC_CONSTEXPR bool is_exact = false;
1784 BOOST_STATIC_CONSTEXPR int radix = 0;
1785 static number_type epsilon() { return number_type(0); }
1786 static number_type round_error() { return number_type(0); }
1787 BOOST_STATIC_CONSTEXPR int min_exponent = 0;
1788 BOOST_STATIC_CONSTEXPR int min_exponent10 = 0;
1789 BOOST_STATIC_CONSTEXPR int max_exponent = 0;
1790 BOOST_STATIC_CONSTEXPR int max_exponent10 = 0;
1791 BOOST_STATIC_CONSTEXPR bool has_infinity = false;
1792 BOOST_STATIC_CONSTEXPR bool has_quiet_NaN = false;
1793 BOOST_STATIC_CONSTEXPR bool has_signaling_NaN = false;
1794 BOOST_STATIC_CONSTEXPR float_denorm_style has_denorm = denorm_absent;
1795 BOOST_STATIC_CONSTEXPR bool has_denorm_loss = false;
1796 static number_type infinity() { return number_type(0); }
1797 static number_type quiet_NaN() { return number_type(0); }
1798 static number_type signaling_NaN() { return number_type(0); }
1799 static number_type denorm_min() { return number_type(0); }
1800 BOOST_STATIC_CONSTEXPR bool is_iec559 = false;
1801 BOOST_STATIC_CONSTEXPR bool is_bounded = false;
1802 BOOST_STATIC_CONSTEXPR bool is_modulo = false;
1803 BOOST_STATIC_CONSTEXPR bool traps = false;
1804 BOOST_STATIC_CONSTEXPR bool tinyness_before = false;
1805 BOOST_STATIC_CONSTEXPR float_round_style round_style = round_toward_zero;
1808 #ifndef BOOST_NO_INCLASS_MEMBER_INITIALIZATION
1810 template <boost::multiprecision::expression_template_option ExpressionTemplates>
1811 BOOST_CONSTEXPR_OR_CONST int numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<0>, ExpressionTemplates> >::digits;
1812 template <boost::multiprecision::expression_template_option ExpressionTemplates>
1813 BOOST_CONSTEXPR_OR_CONST int numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<0>, ExpressionTemplates> >::digits10;
1814 template <boost::multiprecision::expression_template_option ExpressionTemplates>
1815 BOOST_CONSTEXPR_OR_CONST int numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<0>, ExpressionTemplates> >::max_digits10;
1816 template <boost::multiprecision::expression_template_option ExpressionTemplates>
1817 BOOST_CONSTEXPR_OR_CONST bool numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<0>, ExpressionTemplates> >::is_signed;
1818 template <boost::multiprecision::expression_template_option ExpressionTemplates>
1819 BOOST_CONSTEXPR_OR_CONST bool numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<0>, ExpressionTemplates> >::is_integer;
1820 template <boost::multiprecision::expression_template_option ExpressionTemplates>
1821 BOOST_CONSTEXPR_OR_CONST bool numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<0>, ExpressionTemplates> >::is_exact;
1822 template <boost::multiprecision::expression_template_option ExpressionTemplates>
1823 BOOST_CONSTEXPR_OR_CONST int numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<0>, ExpressionTemplates> >::radix;
1824 template <boost::multiprecision::expression_template_option ExpressionTemplates>
1825 BOOST_CONSTEXPR_OR_CONST int numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<0>, ExpressionTemplates> >::min_exponent;
1826 template <boost::multiprecision::expression_template_option ExpressionTemplates>
1827 BOOST_CONSTEXPR_OR_CONST int numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<0>, ExpressionTemplates> >::min_exponent10;
1828 template <boost::multiprecision::expression_template_option ExpressionTemplates>
1829 BOOST_CONSTEXPR_OR_CONST int numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<0>, ExpressionTemplates> >::max_exponent;
1830 template <boost::multiprecision::expression_template_option ExpressionTemplates>
1831 BOOST_CONSTEXPR_OR_CONST int numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<0>, ExpressionTemplates> >::max_exponent10;
1832 template <boost::multiprecision::expression_template_option ExpressionTemplates>
1833 BOOST_CONSTEXPR_OR_CONST bool numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<0>, ExpressionTemplates> >::has_infinity;
1834 template <boost::multiprecision::expression_template_option ExpressionTemplates>
1835 BOOST_CONSTEXPR_OR_CONST bool numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<0>, ExpressionTemplates> >::has_quiet_NaN;
1836 template <boost::multiprecision::expression_template_option ExpressionTemplates>
1837 BOOST_CONSTEXPR_OR_CONST bool numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<0>, ExpressionTemplates> >::has_signaling_NaN;
1838 template <boost::multiprecision::expression_template_option ExpressionTemplates>
1839 BOOST_CONSTEXPR_OR_CONST float_denorm_style numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<0>, ExpressionTemplates> >::has_denorm;
1840 template <boost::multiprecision::expression_template_option ExpressionTemplates>
1841 BOOST_CONSTEXPR_OR_CONST bool numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<0>, ExpressionTemplates> >::has_denorm_loss;
1842 template <boost::multiprecision::expression_template_option ExpressionTemplates>
1843 BOOST_CONSTEXPR_OR_CONST bool numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<0>, ExpressionTemplates> >::is_iec559;
1844 template <boost::multiprecision::expression_template_option ExpressionTemplates>
1845 BOOST_CONSTEXPR_OR_CONST bool numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<0>, ExpressionTemplates> >::is_bounded;
1846 template <boost::multiprecision::expression_template_option ExpressionTemplates>
1847 BOOST_CONSTEXPR_OR_CONST bool numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<0>, ExpressionTemplates> >::is_modulo;
1848 template <boost::multiprecision::expression_template_option ExpressionTemplates>
1849 BOOST_CONSTEXPR_OR_CONST bool numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<0>, ExpressionTemplates> >::traps;
1850 template <boost::multiprecision::expression_template_option ExpressionTemplates>
1851 BOOST_CONSTEXPR_OR_CONST bool numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<0>, ExpressionTemplates> >::tinyness_before;
1852 template <boost::multiprecision::expression_template_option ExpressionTemplates>
1853 BOOST_CONSTEXPR_OR_CONST float_round_style numeric_limits<boost::multiprecision::number<boost::multiprecision::mpfi_float_backend<0>, ExpressionTemplates> >::round_style;