1 ///////////////////////////////////////////////////////////////
2 // Copyright 2012 John Maddock. Distributed under the Boost
3 // Software License, Version 1.0. (See accompanying file
4 // LICENSE_1_0.txt or copy at https://www.boost.org/LICENSE_1_0.txt
6 #ifndef BOOST_MATH_FLOAT_BACKEND_HPP
7 #define BOOST_MATH_FLOAT_BACKEND_HPP
12 #include <boost/cstdint.hpp>
13 #include <boost/lexical_cast.hpp>
14 #include <boost/math/concepts/real_concept.hpp>
15 #include <boost/multiprecision/number.hpp>
16 #include <boost/integer/common_factor_rt.hpp>
17 #include <boost/type_traits/common_type.hpp>
18 #include <boost/container_hash/hash.hpp>
21 namespace multiprecision {
26 #pragma warning(disable : 4389 4244 4018 4244 4127)
29 template <class Arithmetic>
30 struct arithmetic_backend
32 typedef mpl::list<short, int, long, long long> signed_types;
33 typedef mpl::list<unsigned short, unsigned, unsigned long, unsigned long long> unsigned_types;
34 typedef mpl::list<float, double, long double> float_types;
35 typedef int exponent_type;
37 BOOST_MP_CXX14_CONSTEXPR arithmetic_backend() : m_value(0) {}
38 BOOST_MP_CXX14_CONSTEXPR arithmetic_backend(const arithmetic_backend& o) : m_value(o.m_value) {}
40 BOOST_MP_CXX14_CONSTEXPR arithmetic_backend(const A& o, const typename enable_if<is_arithmetic<A> >::type* = 0) : m_value(o) {}
42 BOOST_MP_CXX14_CONSTEXPR arithmetic_backend(const arithmetic_backend<A>& o) : m_value(o.data()) {}
43 BOOST_MP_CXX14_CONSTEXPR arithmetic_backend& operator=(const arithmetic_backend& o)
49 BOOST_MP_CXX14_CONSTEXPR typename enable_if<is_arithmetic<A>, arithmetic_backend&>::type operator=(A i)
55 BOOST_MP_CXX14_CONSTEXPR arithmetic_backend& operator=(const arithmetic_backend<A>& i)
60 arithmetic_backend& operator=(const char* s)
62 #ifndef BOOST_NO_EXCEPTIONS
66 m_value = boost::lexical_cast<Arithmetic>(s);
67 #ifndef BOOST_NO_EXCEPTIONS
69 catch (const bad_lexical_cast&)
71 throw std::runtime_error(std::string("Unable to interpret the string provided: \"") + s + std::string("\" as a compatible number type."));
76 BOOST_MP_CXX14_CONSTEXPR void swap(arithmetic_backend& o)
78 std::swap(m_value, o.m_value);
80 std::string str(std::streamsize digits, std::ios_base::fmtflags f) const
84 ss << std::setprecision(digits ? digits : std::numeric_limits<Arithmetic>::digits10 + 4) << m_value;
87 BOOST_MP_CXX14_CONSTEXPR void do_negate(const mpl::true_&)
89 m_value = 1 + ~m_value;
91 BOOST_MP_CXX14_CONSTEXPR void do_negate(const mpl::false_&)
95 BOOST_MP_CXX14_CONSTEXPR void negate()
97 do_negate(mpl::bool_<is_unsigned<Arithmetic>::value>());
99 BOOST_MP_CXX14_CONSTEXPR int compare(const arithmetic_backend& o) const
101 return m_value > o.m_value ? 1 : (m_value < o.m_value ? -1 : 0);
104 BOOST_MP_CXX14_CONSTEXPR typename enable_if<is_arithmetic<A>, int>::type compare(A i) const
106 return m_value > static_cast<Arithmetic>(i) ? 1 : (m_value < static_cast<Arithmetic>(i) ? -1 : 0);
108 BOOST_MP_CXX14_CONSTEXPR Arithmetic& data() { return m_value; }
109 BOOST_MP_CXX14_CONSTEXPR const Arithmetic& data() const { return m_value; }
115 template <class R, class Arithmetic>
116 inline BOOST_MP_CXX14_CONSTEXPR typename enable_if_c<boost::is_integral<R>::value>::type eval_convert_to(R* result, const arithmetic_backend<Arithmetic>& backend)
118 typedef typename boost::common_type<R, Arithmetic>::type c_type;
119 BOOST_CONSTEXPR const c_type max = static_cast<c_type>((std::numeric_limits<R>::max)());
120 BOOST_CONSTEXPR const c_type min = static_cast<c_type>((std::numeric_limits<R>::min)());
121 c_type ct = static_cast<c_type>(backend.data());
122 if ((backend.data() < 0) && !std::numeric_limits<R>::is_signed)
123 BOOST_THROW_EXCEPTION(std::range_error("Attempt to convert negative number to unsigned type."));
125 *result = boost::is_signed<R>::value ? (std::numeric_limits<R>::max)() : backend.data();
126 else if (std::numeric_limits<Arithmetic>::is_signed && (ct < min))
127 *result = (std::numeric_limits<R>::min)();
129 *result = backend.data();
132 template <class R, class Arithmetic>
133 inline BOOST_MP_CXX14_CONSTEXPR typename disable_if_c<boost::is_integral<R>::value>::type eval_convert_to(R* result, const arithmetic_backend<Arithmetic>& backend)
135 *result = backend.data();
138 template <class Arithmetic>
139 inline BOOST_MP_CXX14_CONSTEXPR bool eval_eq(const arithmetic_backend<Arithmetic>& a, const arithmetic_backend<Arithmetic>& b)
141 return a.data() == b.data();
143 template <class Arithmetic, class A2>
144 inline BOOST_MP_CXX14_CONSTEXPR typename enable_if<is_arithmetic<A2>, bool>::type eval_eq(const arithmetic_backend<Arithmetic>& a, const A2& b)
146 return a.data() == static_cast<Arithmetic>(b);
148 template <class Arithmetic>
149 inline BOOST_MP_CXX14_CONSTEXPR bool eval_lt(const arithmetic_backend<Arithmetic>& a, const arithmetic_backend<Arithmetic>& b)
151 return a.data() < b.data();
153 template <class Arithmetic, class A2>
154 inline BOOST_MP_CXX14_CONSTEXPR typename enable_if<is_arithmetic<A2>, bool>::type eval_lt(const arithmetic_backend<Arithmetic>& a, const A2& b)
156 return a.data() < static_cast<Arithmetic>(b);
158 template <class Arithmetic>
159 inline BOOST_MP_CXX14_CONSTEXPR bool eval_gt(const arithmetic_backend<Arithmetic>& a, const arithmetic_backend<Arithmetic>& b)
161 return a.data() > b.data();
163 template <class Arithmetic, class A2>
164 inline BOOST_MP_CXX14_CONSTEXPR typename enable_if<is_arithmetic<A2>, bool>::type eval_gt(const arithmetic_backend<Arithmetic>& a, const A2& b)
166 return a.data() > static_cast<Arithmetic>(b);
169 template <class Arithmetic>
170 inline BOOST_MP_CXX14_CONSTEXPR void eval_add(arithmetic_backend<Arithmetic>& result, const arithmetic_backend<Arithmetic>& o)
172 result.data() += o.data();
174 template <class Arithmetic>
175 inline BOOST_MP_CXX14_CONSTEXPR void eval_subtract(arithmetic_backend<Arithmetic>& result, const arithmetic_backend<Arithmetic>& o)
177 result.data() -= o.data();
179 template <class Arithmetic>
180 inline BOOST_MP_CXX14_CONSTEXPR void eval_multiply(arithmetic_backend<Arithmetic>& result, const arithmetic_backend<Arithmetic>& o)
182 result.data() *= o.data();
184 template <class Arithmetic>
185 inline BOOST_MP_CXX14_CONSTEXPR typename enable_if_c<std::numeric_limits<Arithmetic>::has_infinity>::type eval_divide(arithmetic_backend<Arithmetic>& result, const arithmetic_backend<Arithmetic>& o)
187 result.data() /= o.data();
189 template <class Arithmetic>
190 inline BOOST_MP_CXX14_CONSTEXPR typename disable_if_c<std::numeric_limits<Arithmetic>::has_infinity>::type eval_divide(arithmetic_backend<Arithmetic>& result, const arithmetic_backend<Arithmetic>& o)
193 BOOST_THROW_EXCEPTION(std::overflow_error("Divide by zero"));
194 result.data() /= o.data();
197 template <class Arithmetic, class A2>
198 inline BOOST_MP_CXX14_CONSTEXPR typename enable_if<is_arithmetic<A2> >::type eval_add(arithmetic_backend<Arithmetic>& result, const A2& o)
202 template <class Arithmetic, class A2>
203 inline BOOST_MP_CXX14_CONSTEXPR typename enable_if<is_arithmetic<A2> >::type eval_subtract(arithmetic_backend<Arithmetic>& result, const A2& o)
207 template <class Arithmetic, class A2>
208 inline BOOST_MP_CXX14_CONSTEXPR typename enable_if<is_arithmetic<A2> >::type eval_multiply(arithmetic_backend<Arithmetic>& result, const A2& o)
212 template <class Arithmetic, class A2>
213 inline BOOST_MP_CXX14_CONSTEXPR typename enable_if_c<(is_arithmetic<A2>::value && !std::numeric_limits<Arithmetic>::has_infinity)>::type
214 eval_divide(arithmetic_backend<Arithmetic>& result, const A2& o)
217 BOOST_THROW_EXCEPTION(std::overflow_error("Divide by zero"));
220 template <class Arithmetic, class A2>
221 inline BOOST_MP_CXX14_CONSTEXPR typename enable_if_c<(is_arithmetic<A2>::value && std::numeric_limits<Arithmetic>::has_infinity)>::type
222 eval_divide(arithmetic_backend<Arithmetic>& result, const A2& o)
227 template <class Arithmetic>
228 inline BOOST_MP_CXX14_CONSTEXPR void eval_add(arithmetic_backend<Arithmetic>& result, const arithmetic_backend<Arithmetic>& a, const arithmetic_backend<Arithmetic>& b)
230 result.data() = a.data() + b.data();
232 template <class Arithmetic>
233 inline BOOST_MP_CXX14_CONSTEXPR void eval_subtract(arithmetic_backend<Arithmetic>& result, const arithmetic_backend<Arithmetic>& a, const arithmetic_backend<Arithmetic>& b)
235 result.data() = a.data() - b.data();
237 template <class Arithmetic>
238 inline BOOST_MP_CXX14_CONSTEXPR void eval_multiply(arithmetic_backend<Arithmetic>& result, const arithmetic_backend<Arithmetic>& a, const arithmetic_backend<Arithmetic>& b)
240 result.data() = a.data() * b.data();
242 template <class Arithmetic>
243 inline BOOST_MP_CXX14_CONSTEXPR typename enable_if_c<std::numeric_limits<Arithmetic>::has_infinity>::type eval_divide(arithmetic_backend<Arithmetic>& result, const arithmetic_backend<Arithmetic>& a, const arithmetic_backend<Arithmetic>& b)
245 result.data() = a.data() / b.data();
247 template <class Arithmetic>
248 inline BOOST_MP_CXX14_CONSTEXPR typename disable_if_c<std::numeric_limits<Arithmetic>::has_infinity>::type eval_divide(arithmetic_backend<Arithmetic>& result, const arithmetic_backend<Arithmetic>& a, const arithmetic_backend<Arithmetic>& b)
251 BOOST_THROW_EXCEPTION(std::overflow_error("Divide by zero"));
252 result.data() = a.data() / b.data();
255 template <class Arithmetic, class A2>
256 inline BOOST_MP_CXX14_CONSTEXPR typename enable_if<is_arithmetic<A2> >::type eval_add(arithmetic_backend<Arithmetic>& result, const arithmetic_backend<Arithmetic>& a, const A2& b)
258 result.data() = a.data() + b;
260 template <class Arithmetic, class A2>
261 inline BOOST_MP_CXX14_CONSTEXPR typename enable_if<is_arithmetic<A2> >::type eval_subtract(arithmetic_backend<Arithmetic>& result, const arithmetic_backend<Arithmetic>& a, const A2& b)
263 result.data() = a.data() - b;
265 template <class Arithmetic, class A2>
266 inline BOOST_MP_CXX14_CONSTEXPR typename enable_if<is_arithmetic<A2> >::type eval_multiply(arithmetic_backend<Arithmetic>& result, const arithmetic_backend<Arithmetic>& a, const A2& b)
268 result.data() = a.data() * b;
270 template <class Arithmetic, class A2>
271 inline BOOST_MP_CXX14_CONSTEXPR typename enable_if_c<(is_arithmetic<A2>::value && !std::numeric_limits<Arithmetic>::has_infinity)>::type
272 eval_divide(arithmetic_backend<Arithmetic>& result, const arithmetic_backend<Arithmetic>& a, const A2& b)
275 BOOST_THROW_EXCEPTION(std::overflow_error("Divide by zero"));
276 result.data() = a.data() / b;
278 template <class Arithmetic, class A2>
279 inline BOOST_MP_CXX14_CONSTEXPR typename enable_if_c<(is_arithmetic<A2>::value && std::numeric_limits<Arithmetic>::has_infinity)>::type
280 eval_divide(arithmetic_backend<Arithmetic>& result, const arithmetic_backend<Arithmetic>& a, const A2& b)
282 result.data() = a.data() / b;
285 template <class Arithmetic>
286 inline BOOST_MP_CXX14_CONSTEXPR bool eval_is_zero(const arithmetic_backend<Arithmetic>& val)
288 return val.data() == 0;
291 template <class Arithmetic>
292 inline BOOST_MP_CXX14_CONSTEXPR typename enable_if_c<
293 (!std::numeric_limits<Arithmetic>::is_specialized || std::numeric_limits<Arithmetic>::is_signed), int>::type
294 eval_get_sign(const arithmetic_backend<Arithmetic>& val)
296 return val.data() == 0 ? 0 : val.data() < 0 ? -1 : 1;
298 template <class Arithmetic>
299 inline BOOST_MP_CXX14_CONSTEXPR typename disable_if_c<
300 (std::numeric_limits<Arithmetic>::is_specialized || std::numeric_limits<Arithmetic>::is_signed), int>::type
301 eval_get_sign(const arithmetic_backend<Arithmetic>& val)
303 return val.data() == 0 ? 0 : 1;
307 inline BOOST_MP_CXX14_CONSTEXPR typename enable_if<is_unsigned<T>, T>::type abs(T v) { return v; }
309 template <class Arithmetic>
310 inline BOOST_MP_CXX14_CONSTEXPR void eval_abs(arithmetic_backend<Arithmetic>& result, const arithmetic_backend<Arithmetic>& o)
312 using boost::multiprecision::backends::abs;
314 result.data() = abs(o.data());
317 template <class Arithmetic>
318 inline BOOST_MP_CXX14_CONSTEXPR void eval_fabs(arithmetic_backend<Arithmetic>& result, const arithmetic_backend<Arithmetic>& o)
320 result.data() = std::abs(o.data());
323 template <class Arithmetic>
324 inline BOOST_MP_CXX14_CONSTEXPR void eval_floor(arithmetic_backend<Arithmetic>& result, const arithmetic_backend<Arithmetic>& o)
327 result.data() = floor(o.data());
330 template <class Arithmetic>
331 inline BOOST_MP_CXX14_CONSTEXPR void eval_ceil(arithmetic_backend<Arithmetic>& result, const arithmetic_backend<Arithmetic>& o)
334 result.data() = ceil(o.data());
337 template <class Arithmetic>
338 inline BOOST_MP_CXX14_CONSTEXPR void eval_sqrt(arithmetic_backend<Arithmetic>& result, const arithmetic_backend<Arithmetic>& o)
341 result.data() = sqrt(o.data());
344 template <class Arithmetic>
345 inline BOOST_MP_CXX14_CONSTEXPR int eval_fpclassify(const arithmetic_backend<Arithmetic>& o)
347 return (boost::math::fpclassify)(o.data());
350 template <class Arithmetic>
351 inline BOOST_MP_CXX14_CONSTEXPR void eval_trunc(arithmetic_backend<Arithmetic>& result, const arithmetic_backend<Arithmetic>& o)
354 result.data() = trunc(o.data());
357 template <class Arithmetic>
358 inline BOOST_MP_CXX14_CONSTEXPR void eval_round(arithmetic_backend<Arithmetic>& result, const arithmetic_backend<Arithmetic>& o)
361 result.data() = round(o.data());
364 template <class Arithmetic>
365 inline BOOST_MP_CXX14_CONSTEXPR void eval_frexp(arithmetic_backend<Arithmetic>& result, const arithmetic_backend<Arithmetic>& a, int* v)
368 result.data() = frexp(a.data(), v);
371 template <class Arithmetic>
372 inline BOOST_MP_CXX14_CONSTEXPR void eval_ldexp(arithmetic_backend<Arithmetic>& result, const arithmetic_backend<Arithmetic>& a, int v)
375 result.data() = ldexp(a.data(), v);
378 template <class Arithmetic>
379 inline BOOST_MP_CXX14_CONSTEXPR void eval_exp(arithmetic_backend<Arithmetic>& result, const arithmetic_backend<Arithmetic>& o)
382 result.data() = exp(o.data());
385 template <class Arithmetic>
386 inline BOOST_MP_CXX14_CONSTEXPR void eval_log(arithmetic_backend<Arithmetic>& result, const arithmetic_backend<Arithmetic>& o)
389 result.data() = log(o.data());
392 template <class Arithmetic>
393 inline BOOST_MP_CXX14_CONSTEXPR void eval_log10(arithmetic_backend<Arithmetic>& result, const arithmetic_backend<Arithmetic>& o)
396 result.data() = log10(o.data());
399 template <class Arithmetic>
400 inline BOOST_MP_CXX14_CONSTEXPR void eval_sin(arithmetic_backend<Arithmetic>& result, const arithmetic_backend<Arithmetic>& o)
403 result.data() = sin(o.data());
406 template <class Arithmetic>
407 inline BOOST_MP_CXX14_CONSTEXPR void eval_cos(arithmetic_backend<Arithmetic>& result, const arithmetic_backend<Arithmetic>& o)
410 result.data() = cos(o.data());
413 template <class Arithmetic>
414 inline BOOST_MP_CXX14_CONSTEXPR void eval_tan(arithmetic_backend<Arithmetic>& result, const arithmetic_backend<Arithmetic>& o)
417 result.data() = tan(o.data());
420 template <class Arithmetic>
421 inline BOOST_MP_CXX14_CONSTEXPR void eval_acos(arithmetic_backend<Arithmetic>& result, const arithmetic_backend<Arithmetic>& o)
424 result.data() = acos(o.data());
427 template <class Arithmetic>
428 inline BOOST_MP_CXX14_CONSTEXPR void eval_asin(arithmetic_backend<Arithmetic>& result, const arithmetic_backend<Arithmetic>& o)
431 result.data() = asin(o.data());
434 template <class Arithmetic>
435 inline BOOST_MP_CXX14_CONSTEXPR void eval_atan(arithmetic_backend<Arithmetic>& result, const arithmetic_backend<Arithmetic>& o)
438 result.data() = atan(o.data());
441 template <class Arithmetic>
442 inline BOOST_MP_CXX14_CONSTEXPR void eval_sinh(arithmetic_backend<Arithmetic>& result, const arithmetic_backend<Arithmetic>& o)
445 result.data() = sinh(o.data());
448 template <class Arithmetic>
449 inline BOOST_MP_CXX14_CONSTEXPR void eval_cosh(arithmetic_backend<Arithmetic>& result, const arithmetic_backend<Arithmetic>& o)
452 result.data() = cosh(o.data());
455 template <class Arithmetic>
456 inline BOOST_MP_CXX14_CONSTEXPR void eval_tanh(arithmetic_backend<Arithmetic>& result, const arithmetic_backend<Arithmetic>& o)
459 result.data() = tanh(o.data());
462 template <class Arithmetic>
463 inline BOOST_MP_CXX14_CONSTEXPR void eval_fmod(arithmetic_backend<Arithmetic>& result, const arithmetic_backend<Arithmetic>& a, const arithmetic_backend<Arithmetic>& b)
466 result.data() = fmod(a.data(), b.data());
469 template <class Arithmetic>
470 inline BOOST_MP_CXX14_CONSTEXPR void eval_pow(arithmetic_backend<Arithmetic>& result, const arithmetic_backend<Arithmetic>& a, const arithmetic_backend<Arithmetic>& b)
473 result.data() = pow(a.data(), b.data());
476 template <class Arithmetic>
477 inline BOOST_MP_CXX14_CONSTEXPR void eval_atan2(arithmetic_backend<Arithmetic>& result, const arithmetic_backend<Arithmetic>& a, const arithmetic_backend<Arithmetic>& b)
480 result.data() = atan2(a.data(), b.data());
483 template <class Arithmetic, class I>
484 inline BOOST_MP_CXX14_CONSTEXPR void eval_left_shift(arithmetic_backend<Arithmetic>& result, I val)
486 result.data() <<= val;
489 template <class Arithmetic, class I>
490 inline BOOST_MP_CXX14_CONSTEXPR void eval_right_shift(arithmetic_backend<Arithmetic>& result, I val)
492 result.data() >>= val;
495 template <class Arithmetic>
496 inline BOOST_MP_CXX14_CONSTEXPR void eval_modulus(arithmetic_backend<Arithmetic>& result, const arithmetic_backend<Arithmetic>& a)
498 result.data() %= a.data();
501 template <class Arithmetic>
502 inline BOOST_MP_CXX14_CONSTEXPR void eval_bitwise_and(arithmetic_backend<Arithmetic>& result, const arithmetic_backend<Arithmetic>& a)
504 result.data() &= a.data();
507 template <class Arithmetic>
508 inline BOOST_MP_CXX14_CONSTEXPR void eval_bitwise_or(arithmetic_backend<Arithmetic>& result, const arithmetic_backend<Arithmetic>& a)
510 result.data() |= a.data();
513 template <class Arithmetic>
514 inline BOOST_MP_CXX14_CONSTEXPR void eval_bitwise_xor(arithmetic_backend<Arithmetic>& result, const arithmetic_backend<Arithmetic>& a)
516 result.data() ^= a.data();
519 template <class Arithmetic>
520 inline BOOST_MP_CXX14_CONSTEXPR void eval_complement(arithmetic_backend<Arithmetic>& result, const arithmetic_backend<Arithmetic>& a)
522 result.data() = ~a.data();
525 template <class Arithmetic>
526 inline BOOST_MP_CXX14_CONSTEXPR void eval_gcd(arithmetic_backend<Arithmetic>& result, const arithmetic_backend<Arithmetic>& a, const arithmetic_backend<Arithmetic>& b)
528 result.data() = boost::integer::gcd(a.data(), b.data());
531 template <class Arithmetic>
532 inline BOOST_MP_CXX14_CONSTEXPR void eval_lcm(arithmetic_backend<Arithmetic>& result, const arithmetic_backend<Arithmetic>& a, const arithmetic_backend<Arithmetic>& b)
534 result.data() = boost::integer::lcm(a.data(), b.data());
537 template <class Arithmetic>
538 inline BOOST_MP_CXX14_CONSTEXPR std::size_t hash_value(const arithmetic_backend<Arithmetic>& a)
540 boost::hash<Arithmetic> hasher;
541 return hasher(a.data());
548 } // namespace backends
550 using boost::multiprecision::backends::arithmetic_backend;
552 template <class Arithmetic>
553 struct number_category<arithmetic_backend<Arithmetic> > : public mpl::int_<is_integral<Arithmetic>::value ? number_kind_integer : number_kind_floating_point>
558 template <class Backend>
559 struct double_precision_type;
561 template <class Arithmetic, boost::multiprecision::expression_template_option ET>
562 struct double_precision_type<number<arithmetic_backend<Arithmetic>, ET> >
564 typedef number<arithmetic_backend<typename double_precision_type<Arithmetic>::type>, ET> type;
567 struct double_precision_type<arithmetic_backend<boost::int32_t> >
569 typedef arithmetic_backend<boost::int64_t> type;
572 } // namespace detail
574 }} // namespace boost::multiprecision
575 #if !(defined(__SGI_STL_PORT) || defined(BOOST_NO_LIMITS_COMPILE_TIME_CONSTANTS))
577 // We shouldn't need these to get code to compile, however for the sake of
578 // "level playing field" performance comparisons they avoid the very slow
579 // lexical_cast's that would otherwise take place. Definition has to be guarded
580 // by the inverse of pp-logic in real_concept.hpp which defines these as a workaround
581 // for STLPort plus some other old/broken standartd libraries.
583 namespace boost { namespace math { namespace tools {
586 inline unsigned int real_cast<unsigned int, concepts::real_concept>(concepts::real_concept r)
588 return static_cast<unsigned int>(r.value());
592 inline int real_cast<int, concepts::real_concept>(concepts::real_concept r)
594 return static_cast<int>(r.value());
598 inline long real_cast<long, concepts::real_concept>(concepts::real_concept r)
600 return static_cast<long>(r.value());
603 // Converts from T to narrower floating-point types, float, double & long double.
606 inline float real_cast<float, concepts::real_concept>(concepts::real_concept r)
608 return static_cast<float>(r.value());
611 inline double real_cast<double, concepts::real_concept>(concepts::real_concept r)
613 return static_cast<double>(r.value());
616 inline long double real_cast<long double, concepts::real_concept>(concepts::real_concept r)
621 }}} // namespace boost::math::tools
626 template <class Arithmetic, boost::multiprecision::expression_template_option ExpressionTemplates>
627 class numeric_limits<boost::multiprecision::number<boost::multiprecision::arithmetic_backend<Arithmetic>, ExpressionTemplates> > : public std::numeric_limits<Arithmetic>
629 typedef std::numeric_limits<Arithmetic> base_type;
630 typedef boost::multiprecision::number<boost::multiprecision::arithmetic_backend<Arithmetic>, ExpressionTemplates> number_type;
633 BOOST_STATIC_CONSTEXPR number_type(min)() BOOST_NOEXCEPT { return (base_type::min)(); }
634 BOOST_STATIC_CONSTEXPR number_type(max)() BOOST_NOEXCEPT { return (base_type::max)(); }
635 BOOST_STATIC_CONSTEXPR number_type lowest() BOOST_NOEXCEPT { return -(max)(); }
636 BOOST_STATIC_CONSTEXPR number_type epsilon() BOOST_NOEXCEPT { return base_type::epsilon(); }
637 BOOST_STATIC_CONSTEXPR number_type round_error() BOOST_NOEXCEPT { return epsilon() / 2; }
638 BOOST_STATIC_CONSTEXPR number_type infinity() BOOST_NOEXCEPT { return base_type::infinity(); }
639 BOOST_STATIC_CONSTEXPR number_type quiet_NaN() BOOST_NOEXCEPT { return base_type::quiet_NaN(); }
640 BOOST_STATIC_CONSTEXPR number_type signaling_NaN() BOOST_NOEXCEPT { return base_type::signaling_NaN(); }
641 BOOST_STATIC_CONSTEXPR number_type denorm_min() BOOST_NOEXCEPT { return base_type::denorm_min(); }
645 class numeric_limits<boost::math::concepts::real_concept> : public std::numeric_limits<long double>
647 typedef std::numeric_limits<long double> base_type;
648 typedef boost::math::concepts::real_concept number_type;
651 static const number_type(min)() BOOST_NOEXCEPT { return (base_type::min)(); }
652 static const number_type(max)() BOOST_NOEXCEPT { return (base_type::max)(); }
653 static const number_type lowest() BOOST_NOEXCEPT { return -(max)(); }
654 static const number_type epsilon() BOOST_NOEXCEPT { return base_type::epsilon(); }
655 static const number_type round_error() BOOST_NOEXCEPT { return epsilon() / 2; }
656 static const number_type infinity() BOOST_NOEXCEPT { return base_type::infinity(); }
657 static const number_type quiet_NaN() BOOST_NOEXCEPT { return base_type::quiet_NaN(); }
658 static const number_type signaling_NaN() BOOST_NOEXCEPT { return base_type::signaling_NaN(); }
659 static const number_type denorm_min() BOOST_NOEXCEPT { return base_type::denorm_min(); }
664 #include <boost/multiprecision/detail/integer_ops.hpp>