1///////////////////////////////////////////////////////////////////////////////
2// Copyright 2013 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)
5
6#ifndef BOOST_MP_UBLAS_INTEROP_HPP
7#define BOOST_MP_UBLAS_INTEROP_HPP
8
9namespace boost { namespace numeric { namespace ublas {
10
11template <class V>
12class sparse_vector_element;
13
14template <class V, class Backend, multiprecision::expression_template_option ExpressionTemplates>
15inline bool operator==(const sparse_vector_element<V>& a, const ::boost::multiprecision::number<Backend, ExpressionTemplates>& b)
16{
17 using ref_type = typename sparse_vector_element<V>::const_reference;
18 return static_cast<ref_type>(a) == b;
19}
20
21template <class X, class Y>
22struct promote_traits;
23
24template <class Backend1, boost::multiprecision::expression_template_option ExpressionTemplates1, class Backend2, boost::multiprecision::expression_template_option ExpressionTemplates2>
25struct promote_traits<boost::multiprecision::number<Backend1, ExpressionTemplates1>, boost::multiprecision::number<Backend2, ExpressionTemplates2> >
26{
27 using number1_t = boost::multiprecision::number<Backend1, ExpressionTemplates1>;
28 using number2_t = boost::multiprecision::number<Backend2, ExpressionTemplates2>;
29 using promote_type = typename std::conditional<
30 std::is_convertible<number1_t, number2_t>::value && !std::is_convertible<number2_t, number1_t>::value,
31 number2_t, number1_t>::type;
32};
33
34template <class Backend1, boost::multiprecision::expression_template_option ExpressionTemplates1, class Arithmetic>
35struct promote_traits<boost::multiprecision::number<Backend1, ExpressionTemplates1>, Arithmetic>
36{
37 using promote_type = boost::multiprecision::number<Backend1, ExpressionTemplates1>;
38};
39
40template <class Arithmetic, class Backend1, boost::multiprecision::expression_template_option ExpressionTemplates1>
41struct promote_traits<Arithmetic, boost::multiprecision::number<Backend1, ExpressionTemplates1> >
42{
43 using promote_type = boost::multiprecision::number<Backend1, ExpressionTemplates1>;
44};
45
46template <class Backend1, boost::multiprecision::expression_template_option ExpressionTemplates1, class tag, class Arg1, class Arg2, class Arg3, class Arg4>
47struct promote_traits<boost::multiprecision::number<Backend1, ExpressionTemplates1>, boost::multiprecision::detail::expression<tag, Arg1, Arg2, Arg3, Arg4> >
48{
49 using number1_t = boost::multiprecision::number<Backend1, ExpressionTemplates1> ;
50 using expression_type = boost::multiprecision::detail::expression<tag, Arg1, Arg2, Arg3, Arg4>;
51 using number2_t = typename expression_type::result_type ;
52 using promote_type = typename promote_traits<number1_t, number2_t>::promote_type ;
53};
54
55template <class tag, class Arg1, class Arg2, class Arg3, class Arg4, class Backend1, boost::multiprecision::expression_template_option ExpressionTemplates1>
56struct promote_traits<boost::multiprecision::detail::expression<tag, Arg1, Arg2, Arg3, Arg4>, boost::multiprecision::number<Backend1, ExpressionTemplates1> >
57{
58 using number1_t = boost::multiprecision::number<Backend1, ExpressionTemplates1> ;
59 using expression_type = boost::multiprecision::detail::expression<tag, Arg1, Arg2, Arg3, Arg4>;
60 using number2_t = typename expression_type::result_type ;
61 using promote_type = typename promote_traits<number1_t, number2_t>::promote_type ;
62};
63
64template <class tag, class Arg1, class Arg2, class Arg3, class Arg4, class tagb, class Arg1b, class Arg2b, class Arg3b, class Arg4b>
65struct promote_traits<boost::multiprecision::detail::expression<tag, Arg1, Arg2, Arg3, Arg4>, boost::multiprecision::detail::expression<tagb, Arg1b, Arg2b, Arg3b, Arg4b> >
66{
67 using expression1_t = boost::multiprecision::detail::expression<tag, Arg1, Arg2, Arg3, Arg4> ;
68 using number1_t = typename expression1_t::result_type ;
69 using expression2_t = boost::multiprecision::detail::expression<tagb, Arg1b, Arg2b, Arg3b, Arg4b>;
70 using number2_t = typename expression2_t::result_type ;
71};
72
73}}} // namespace boost::numeric::ublas
74
75#endif
76

source code of boost/libs/multiprecision/include/boost/multiprecision/detail/ublas_interop.hpp