Imported Upstream version 1.72.0
[platform/upstream/boost.git] / libs / multiprecision / test / test_arithmetic.hpp
index 5c44eff..58699fd 100644 (file)
 #include "test.hpp"
 
 template <class T>
-struct is_boost_rational : public boost::mpl::false_{};
+struct is_boost_rational : public boost::mpl::false_
+{};
 template <class T>
-struct is_checked_cpp_int : public boost::mpl::false_ {};
+struct is_checked_cpp_int : public boost::mpl::false_
+{};
 
 #ifdef BOOST_MSVC
 // warning C4127: conditional expression is constant
-#pragma warning(disable:4127)
+#pragma warning(disable : 4127)
 #endif
 
 template <class Target, class Source>
@@ -33,7 +35,7 @@ Target checked_lexical_cast(const Source& val)
       return boost::lexical_cast<Target>(val);
 #ifndef BOOST_NO_EXCEPTIONS
    }
-   catch(...)
+   catch (...)
    {
       std::cerr << "Error in lexical cast\nSource type = " << typeid(Source).name() << " \"" << val << "\"\n";
       std::cerr << "Target type = " << typeid(Target).name() << std::endl;
@@ -42,26 +44,27 @@ Target checked_lexical_cast(const Source& val)
 #endif
 }
 
+bool isfloat(float) { return true; }
+bool isfloat(double) { return true; }
+bool isfloat(long double) { return true; }
+template <class T>
+bool isfloat(T) { return false; }
 
-bool isfloat(float){ return true; }
-bool isfloat(double){ return true; }
-bool isfloat(long double){ return true; }
-template <class T> bool isfloat(T){ return false; }
-
-namespace detail{
+namespace detail {
 
-template<class tag, class Arg1, class Arg2, class Arg3, class Arg4>
+template <class tag, class Arg1, class Arg2, class Arg3, class Arg4>
 typename boost::multiprecision::detail::expression<tag, Arg1, Arg2, Arg3, Arg4>::result_type
-   abs(boost::multiprecision::detail::expression<tag, Arg1, Arg2, Arg3, Arg4> const& v)
+abs(boost::multiprecision::detail::expression<tag, Arg1, Arg2, Arg3, Arg4> const& v)
 {
    typedef typename boost::multiprecision::detail::expression<tag, Arg1, Arg2, Arg3, Arg4>::result_type result_type;
    return v < 0 ? result_type(-v) : result_type(v);
 }
 
-}
+} // namespace detail
 
 template <class T>
-struct is_twos_complement_integer : public boost::mpl::true_ {};
+struct is_twos_complement_integer : public boost::mpl::true_
+{};
 
 template <class T>
 struct related_type
@@ -79,8 +82,8 @@ int normalize_compare_result(int r)
 }
 
 template <class Real, class Val>
-typename boost::disable_if_c<boost::multiprecision::number_category<Real>::value == boost::multiprecision::number_kind_complex>::type 
-   test_comparisons(Val a, Val b, const boost::mpl::true_)
+typename boost::disable_if_c<boost::multiprecision::number_category<Real>::value == boost::multiprecision::number_kind_complex>::type
+test_comparisons(Val a, Val b, const boost::mpl::true_)
 {
    Real r1(a);
    Real r2(b);
@@ -109,40 +112,40 @@ typename boost::disable_if_c<boost::multiprecision::number_category<Real>::value
    BOOST_CHECK_EQUAL(a >= r2, a >= b);
    BOOST_CHECK_EQUAL(a > r2, a > b);
 
-   BOOST_CHECK_EQUAL(r1*z == r2, a == b);
-   BOOST_CHECK_EQUAL(r1*z != r2, a != b);
-   BOOST_CHECK_EQUAL(r1*z <= r2, a <= b);
-   BOOST_CHECK_EQUAL(r1*z < r2, a < b);
-   BOOST_CHECK_EQUAL(r1*z >= r2, a >= b);
-   BOOST_CHECK_EQUAL(r1*z > r2, a > b);
-
-   BOOST_CHECK_EQUAL(r1 == r2*z, a == b);
-   BOOST_CHECK_EQUAL(r1 != r2*z, a != b);
-   BOOST_CHECK_EQUAL(r1 <= r2*z, a <= b);
-   BOOST_CHECK_EQUAL(r1 < r2*z, a < b);
-   BOOST_CHECK_EQUAL(r1 >= r2*z, a >= b);
-   BOOST_CHECK_EQUAL(r1 > r2*z, a > b);
-
-   BOOST_CHECK_EQUAL(r1*z == r2*z, a == b);
-   BOOST_CHECK_EQUAL(r1*z != r2*z, a != b);
-   BOOST_CHECK_EQUAL(r1*z <= r2*z, a <= b);
-   BOOST_CHECK_EQUAL(r1*z < r2*z, a < b);
-   BOOST_CHECK_EQUAL(r1*z >= r2*z, a >= b);
-   BOOST_CHECK_EQUAL(r1*z > r2*z, a > b);
-
-   BOOST_CHECK_EQUAL(r1*z == b, a == b);
-   BOOST_CHECK_EQUAL(r1*z != b, a != b);
-   BOOST_CHECK_EQUAL(r1*z <= b, a <= b);
-   BOOST_CHECK_EQUAL(r1*z < b, a < b);
-   BOOST_CHECK_EQUAL(r1*z >= b, a >= b);
-   BOOST_CHECK_EQUAL(r1*z > b, a > b);
-
-   BOOST_CHECK_EQUAL(a == r2*z, a == b);
-   BOOST_CHECK_EQUAL(a != r2*z, a != b);
-   BOOST_CHECK_EQUAL(a <= r2*z, a <= b);
-   BOOST_CHECK_EQUAL(a < r2*z, a < b);
-   BOOST_CHECK_EQUAL(a >= r2*z, a >= b);
-   BOOST_CHECK_EQUAL(a > r2*z, a > b);
+   BOOST_CHECK_EQUAL(r1 * z == r2, a == b);
+   BOOST_CHECK_EQUAL(r1 * z != r2, a != b);
+   BOOST_CHECK_EQUAL(r1 * z <= r2, a <= b);
+   BOOST_CHECK_EQUAL(r1 * z < r2, a < b);
+   BOOST_CHECK_EQUAL(r1 * z >= r2, a >= b);
+   BOOST_CHECK_EQUAL(r1 * z > r2, a > b);
+
+   BOOST_CHECK_EQUAL(r1 == r2 * z, a == b);
+   BOOST_CHECK_EQUAL(r1 != r2 * z, a != b);
+   BOOST_CHECK_EQUAL(r1 <= r2 * z, a <= b);
+   BOOST_CHECK_EQUAL(r1 < r2 * z, a < b);
+   BOOST_CHECK_EQUAL(r1 >= r2 * z, a >= b);
+   BOOST_CHECK_EQUAL(r1 > r2 * z, a > b);
+
+   BOOST_CHECK_EQUAL(r1 * z == r2 * z, a == b);
+   BOOST_CHECK_EQUAL(r1 * z != r2 * z, a != b);
+   BOOST_CHECK_EQUAL(r1 * z <= r2 * z, a <= b);
+   BOOST_CHECK_EQUAL(r1 * z < r2 * z, a < b);
+   BOOST_CHECK_EQUAL(r1 * z >= r2 * z, a >= b);
+   BOOST_CHECK_EQUAL(r1 * z > r2 * z, a > b);
+
+   BOOST_CHECK_EQUAL(r1 * z == b, a == b);
+   BOOST_CHECK_EQUAL(r1 * z != b, a != b);
+   BOOST_CHECK_EQUAL(r1 * z <= b, a <= b);
+   BOOST_CHECK_EQUAL(r1 * z < b, a < b);
+   BOOST_CHECK_EQUAL(r1 * z >= b, a >= b);
+   BOOST_CHECK_EQUAL(r1 * z > b, a > b);
+
+   BOOST_CHECK_EQUAL(a == r2 * z, a == b);
+   BOOST_CHECK_EQUAL(a != r2 * z, a != b);
+   BOOST_CHECK_EQUAL(a <= r2 * z, a <= b);
+   BOOST_CHECK_EQUAL(a < r2 * z, a < b);
+   BOOST_CHECK_EQUAL(a >= r2 * z, a >= b);
+   BOOST_CHECK_EQUAL(a > r2 * z, a > b);
 
    BOOST_CHECK_EQUAL(normalize_compare_result(r1.compare(r2)), cr);
    BOOST_CHECK_EQUAL(normalize_compare_result(r2.compare(r1)), -cr);
@@ -152,7 +155,7 @@ typename boost::disable_if_c<boost::multiprecision::number_category<Real>::value
 
 template <class Real, class Val>
 typename boost::enable_if_c<boost::multiprecision::number_category<Real>::value == boost::multiprecision::number_kind_complex>::type
-  test_comparisons(Val a, Val b, const boost::mpl::true_)
+test_comparisons(Val a, Val b, const boost::mpl::true_)
 {
    Real r1(a);
    Real r2(b);
@@ -170,17 +173,17 @@ typename boost::enable_if_c<boost::multiprecision::number_category<Real>::value
    BOOST_CHECK_EQUAL(a == r2, a == b);
    BOOST_CHECK_EQUAL(a != r2, a != b);
 
-   BOOST_CHECK_EQUAL(r1*z == r2, a == b);
-   BOOST_CHECK_EQUAL(r1*z != r2, a != b);
+   BOOST_CHECK_EQUAL(r1 * z == r2, a == b);
+   BOOST_CHECK_EQUAL(r1 * z != r2, a != b);
 
    BOOST_CHECK_EQUAL(r1 == r2 * z, a == b);
    BOOST_CHECK_EQUAL(r1 != r2 * z, a != b);
 
-   BOOST_CHECK_EQUAL(r1*z == r2 * z, a == b);
-   BOOST_CHECK_EQUAL(r1*z != r2 * z, a != b);
+   BOOST_CHECK_EQUAL(r1 * z == r2 * z, a == b);
+   BOOST_CHECK_EQUAL(r1 * z != r2 * z, a != b);
 
-   BOOST_CHECK_EQUAL(r1*z == b, a == b);
-   BOOST_CHECK_EQUAL(r1*z != b, a != b);
+   BOOST_CHECK_EQUAL(r1 * z == b, a == b);
+   BOOST_CHECK_EQUAL(r1 * z != b, a != b);
 
    BOOST_CHECK_EQUAL(a == r2 * z, a == b);
    BOOST_CHECK_EQUAL(a != r2 * z, a != b);
@@ -207,7 +210,7 @@ void test_conditional(Real v, Exp e)
    //
    // Verify that Exp is usable in Boolean contexts, and has the same value as v:
    //
-   if(e)
+   if (e)
    {
       BOOST_CHECK(v);
    }
@@ -215,7 +218,7 @@ void test_conditional(Real v, Exp e)
    {
       BOOST_CHECK(!v);
    }
-   if(!e)
+   if (!e)
    {
       BOOST_CHECK(!v);
    }
@@ -228,20 +231,20 @@ void test_conditional(Real v, Exp e)
 template <class Real>
 void test_complement(Real a, Real b, Real c, const boost::mpl::true_&)
 {
-   int i = 1020304;
-   int j = 56789123;
+   int i         = 1020304;
+   int j         = 56789123;
    int sign_mask = ~0;
-   if(std::numeric_limits<Real>::is_signed)
+   if (std::numeric_limits<Real>::is_signed)
    {
-      BOOST_CHECK_EQUAL(~a ,  (~i & sign_mask));
+      BOOST_CHECK_EQUAL(~a, (~i & sign_mask));
       c = a & ~b;
-      BOOST_CHECK_EQUAL(c ,  (i & (~j & sign_mask)));
+      BOOST_CHECK_EQUAL(c, (i & (~j & sign_mask)));
       c = ~(a | b);
-      BOOST_CHECK_EQUAL(c ,  (~(i | j) & sign_mask));
+      BOOST_CHECK_EQUAL(c, (~(i | j) & sign_mask));
    }
    else
    {
-      BOOST_CHECK_EQUAL((~a & a) ,  0);
+      BOOST_CHECK_EQUAL((~a & a), 0);
    }
 }
 
@@ -251,18 +254,18 @@ void test_complement(Real, Real, Real, const boost::mpl::false_&)
 }
 
 template <class Real, class T>
-void test_integer_ops(const T&){}
+void test_integer_ops(const T&) {}
 
 template <class Real>
 void test_rational(const boost::mpl::true_&)
 {
    Real a(2);
    a /= 3;
-   BOOST_CHECK_EQUAL(numerator(a) ,  2);
-   BOOST_CHECK_EQUAL(denominator(a) ,  3);
+   BOOST_CHECK_EQUAL(numerator(a), 2);
+   BOOST_CHECK_EQUAL(denominator(a), 3);
    Real b(4);
    b /= 6;
-   BOOST_CHECK_EQUAL(a ,  b);
+   BOOST_CHECK_EQUAL(a, b);
 
    //
    // Check IO code:
@@ -278,11 +281,11 @@ void test_rational(const boost::mpl::false_&)
 {
    Real a(2);
    a /= 3;
-   BOOST_CHECK_EQUAL(numerator(a) ,  2);
-   BOOST_CHECK_EQUAL(denominator(a) ,  3);
+   BOOST_CHECK_EQUAL(numerator(a), 2);
+   BOOST_CHECK_EQUAL(denominator(a), 3);
    Real b(4);
    b /= 6;
-   BOOST_CHECK_EQUAL(a ,  b);
+   BOOST_CHECK_EQUAL(a, b);
 
 #ifndef BOOST_NO_EXCEPTIONS
    BOOST_CHECK_THROW(Real(a / 0), std::overflow_error);
@@ -311,103 +314,103 @@ void test_signed_integer_ops(const boost::mpl::true_&)
    Real a(20);
    Real b(7);
    Real c(5);
-   BOOST_CHECK_EQUAL(-a % c ,  0);
-   BOOST_CHECK_EQUAL(-a % b ,  -20 % 7);
-   BOOST_CHECK_EQUAL(-a % -b ,  -20 % -7);
-   BOOST_CHECK_EQUAL(a % -b ,  20 % -7);
-   BOOST_CHECK_EQUAL(-a % 7 ,  -20 % 7);
-   BOOST_CHECK_EQUAL(-a % -7 ,  -20 % -7);
-   BOOST_CHECK_EQUAL(a % -7 ,  20 % -7);
-   BOOST_CHECK_EQUAL(-a % 7u ,  -20 % 7);
-   BOOST_CHECK_EQUAL(-a % a ,  0);
-   BOOST_CHECK_EQUAL(-a % 5 ,  0);
-   BOOST_CHECK_EQUAL(-a % -5 ,  0);
-   BOOST_CHECK_EQUAL(a % -5 ,  0);
+   BOOST_CHECK_EQUAL(-a % c, 0);
+   BOOST_CHECK_EQUAL(-a % b, -20 % 7);
+   BOOST_CHECK_EQUAL(-a % -b, -20 % -7);
+   BOOST_CHECK_EQUAL(a % -b, 20 % -7);
+   BOOST_CHECK_EQUAL(-a % 7, -20 % 7);
+   BOOST_CHECK_EQUAL(-a % -7, -20 % -7);
+   BOOST_CHECK_EQUAL(a % -7, 20 % -7);
+   BOOST_CHECK_EQUAL(-a % 7u, -20 % 7);
+   BOOST_CHECK_EQUAL(-a % a, 0);
+   BOOST_CHECK_EQUAL(-a % 5, 0);
+   BOOST_CHECK_EQUAL(-a % -5, 0);
+   BOOST_CHECK_EQUAL(a % -5, 0);
 
    b = -b;
-   BOOST_CHECK_EQUAL(a % b ,  20 % -7);
+   BOOST_CHECK_EQUAL(a % b, 20 % -7);
    a = -a;
-   BOOST_CHECK_EQUAL(a % b ,  -20 % -7);
-   BOOST_CHECK_EQUAL(a % -7 ,  -20 % -7);
+   BOOST_CHECK_EQUAL(a % b, -20 % -7);
+   BOOST_CHECK_EQUAL(a % -7, -20 % -7);
    b = 7;
-   BOOST_CHECK_EQUAL(a % b ,  -20 % 7);
-   BOOST_CHECK_EQUAL(a % 7 ,  -20 % 7);
-   BOOST_CHECK_EQUAL(a % 7u ,  -20 % 7);
+   BOOST_CHECK_EQUAL(a % b, -20 % 7);
+   BOOST_CHECK_EQUAL(a % 7, -20 % 7);
+   BOOST_CHECK_EQUAL(a % 7u, -20 % 7);
 
    a = 20;
    a %= b;
-   BOOST_CHECK_EQUAL(a ,  20 % 7);
+   BOOST_CHECK_EQUAL(a, 20 % 7);
    a = -20;
    a %= b;
-   BOOST_CHECK_EQUAL(a ,  -20 % 7);
+   BOOST_CHECK_EQUAL(a, -20 % 7);
    a = 20;
    a %= -b;
-   BOOST_CHECK_EQUAL(a ,  20 % -7);
+   BOOST_CHECK_EQUAL(a, 20 % -7);
    a = -20;
    a %= -b;
-   BOOST_CHECK_EQUAL(a ,  -20 % -7);
+   BOOST_CHECK_EQUAL(a, -20 % -7);
    a = 5;
    a %= b - a;
-   BOOST_CHECK_EQUAL(a ,  5 % (7-5));
+   BOOST_CHECK_EQUAL(a, 5 % (7 - 5));
    a = -20;
    a %= 7;
-   BOOST_CHECK_EQUAL(a ,  -20 % 7);
+   BOOST_CHECK_EQUAL(a, -20 % 7);
    a = 20;
    a %= -7;
-   BOOST_CHECK_EQUAL(a ,  20 % -7);
+   BOOST_CHECK_EQUAL(a, 20 % -7);
    a = -20;
    a %= -7;
-   BOOST_CHECK_EQUAL(a ,  -20 % -7);
+   BOOST_CHECK_EQUAL(a, -20 % -7);
 #ifndef BOOST_NO_LONG_LONG
    a = -20;
    a %= 7uLL;
-   BOOST_CHECK_EQUAL(a ,  -20 % 7);
+   BOOST_CHECK_EQUAL(a, -20 % 7);
    a = 20;
    a %= -7LL;
-   BOOST_CHECK_EQUAL(a ,  20 % -7);
+   BOOST_CHECK_EQUAL(a, 20 % -7);
    a = -20;
    a %= -7LL;
-   BOOST_CHECK_EQUAL(a ,  -20 % -7);
+   BOOST_CHECK_EQUAL(a, -20 % -7);
 #endif
    a = 400;
    b = 45;
-   BOOST_CHECK_EQUAL(gcd(a, -45) ,  boost::integer::gcd(400, 45));
-   BOOST_CHECK_EQUAL(lcm(a, -45) ,  boost::integer::lcm(400, 45));
-   BOOST_CHECK_EQUAL(gcd(-400, b) ,  boost::integer::gcd(400, 45));
-   BOOST_CHECK_EQUAL(lcm(-400, b) ,  boost::integer::lcm(400, 45));
+   BOOST_CHECK_EQUAL(gcd(a, -45), boost::integer::gcd(400, 45));
+   BOOST_CHECK_EQUAL(lcm(a, -45), boost::integer::lcm(400, 45));
+   BOOST_CHECK_EQUAL(gcd(-400, b), boost::integer::gcd(400, 45));
+   BOOST_CHECK_EQUAL(lcm(-400, b), boost::integer::lcm(400, 45));
    a = -20;
-   BOOST_CHECK_EQUAL(abs(a) ,  20);
-   BOOST_CHECK_EQUAL(abs(-a) ,  20);
-   BOOST_CHECK_EQUAL(abs(+a) ,  20);
+   BOOST_CHECK_EQUAL(abs(a), 20);
+   BOOST_CHECK_EQUAL(abs(-a), 20);
+   BOOST_CHECK_EQUAL(abs(+a), 20);
    a = 20;
-   BOOST_CHECK_EQUAL(abs(a) ,  20);
-   BOOST_CHECK_EQUAL(abs(-a) ,  20);
-   BOOST_CHECK_EQUAL(abs(+a) ,  20);
+   BOOST_CHECK_EQUAL(abs(a), 20);
+   BOOST_CHECK_EQUAL(abs(-a), 20);
+   BOOST_CHECK_EQUAL(abs(+a), 20);
    a = -400;
    b = 45;
-   BOOST_CHECK_EQUAL(gcd(a, b) ,  boost::integer::gcd(-400, 45));
-   BOOST_CHECK_EQUAL(lcm(a, b) ,  boost::integer::lcm(-400, 45));
-   BOOST_CHECK_EQUAL(gcd(a, 45) ,  boost::integer::gcd(-400, 45));
-   BOOST_CHECK_EQUAL(lcm(a, 45) ,  boost::integer::lcm(-400, 45));
-   BOOST_CHECK_EQUAL(gcd(-400, b) ,  boost::integer::gcd(-400, 45));
-   BOOST_CHECK_EQUAL(lcm(-400, b) ,  boost::integer::lcm(-400, 45));
+   BOOST_CHECK_EQUAL(gcd(a, b), boost::integer::gcd(-400, 45));
+   BOOST_CHECK_EQUAL(lcm(a, b), boost::integer::lcm(-400, 45));
+   BOOST_CHECK_EQUAL(gcd(a, 45), boost::integer::gcd(-400, 45));
+   BOOST_CHECK_EQUAL(lcm(a, 45), boost::integer::lcm(-400, 45));
+   BOOST_CHECK_EQUAL(gcd(-400, b), boost::integer::gcd(-400, 45));
+   BOOST_CHECK_EQUAL(lcm(-400, b), boost::integer::lcm(-400, 45));
    Real r;
    divide_qr(a, b, c, r);
-   BOOST_CHECK_EQUAL(c ,  a / b);
-   BOOST_CHECK_EQUAL(r ,  a % b);
-   BOOST_CHECK_EQUAL(integer_modulus(a, 57) ,  abs(a % 57));
+   BOOST_CHECK_EQUAL(c, a / b);
+   BOOST_CHECK_EQUAL(r, a % b);
+   BOOST_CHECK_EQUAL(integer_modulus(a, 57), abs(a % 57));
    b = -57;
    divide_qr(a, b, c, r);
-   BOOST_CHECK_EQUAL(c ,  a / b);
-   BOOST_CHECK_EQUAL(r ,  a % b);
-   BOOST_CHECK_EQUAL(integer_modulus(a, -57) ,  abs(a % -57));
+   BOOST_CHECK_EQUAL(c, a / b);
+   BOOST_CHECK_EQUAL(r, a % b);
+   BOOST_CHECK_EQUAL(integer_modulus(a, -57), abs(a % -57));
    a = 458;
    divide_qr(a, b, c, r);
-   BOOST_CHECK_EQUAL(c ,  a / b);
-   BOOST_CHECK_EQUAL(r ,  a % b);
-   BOOST_CHECK_EQUAL(integer_modulus(a, -57) ,  abs(a % -57));
+   BOOST_CHECK_EQUAL(c, a / b);
+   BOOST_CHECK_EQUAL(r, a % b);
+   BOOST_CHECK_EQUAL(integer_modulus(a, -57), abs(a % -57));
 #ifndef TEST_CHECKED_INT
-   if(is_checked_cpp_int<Real>::value)
+   if (is_checked_cpp_int<Real>::value)
    {
       a = -1;
 #ifndef BOOST_NO_EXCEPTIONS
@@ -452,13 +455,13 @@ void test_integer_overflow()
    if (std::numeric_limits<Real>::digits > std::numeric_limits<Int>::digits)
    {
       Real m((std::numeric_limits<Int>::max)());
-      Int r;
+      Int  r;
       ++m;
       if (is_checked_cpp_int<Real>::value)
       {
          BOOST_CHECK_THROW(m.template convert_to<Int>(), std::overflow_error);
       }
-      else if(boost::is_signed<Int>::value)
+      else if (boost::is_signed<Int>::value)
       {
          r = m.template convert_to<Int>();
          BOOST_CHECK_EQUAL(r, (std::numeric_limits<Int>::max)());
@@ -533,12 +536,12 @@ void test_integer_overflow()
 template <class Real, class Int>
 void test_integer_round_trip()
 {
-   if(std::numeric_limits<Real>::digits >= std::numeric_limits<Int>::digits)
+   if (std::numeric_limits<Real>::digits >= std::numeric_limits<Int>::digits)
    {
       Real m((std::numeric_limits<Int>::max)());
-      Int r = m.template convert_to<Int>();
+      Int  r = m.template convert_to<Int>();
       BOOST_CHECK_EQUAL(m, r);
-      if(std::numeric_limits<Real>::is_signed && (std::numeric_limits<Real>::digits > std::numeric_limits<Int>::digits))
+      if (std::numeric_limits<Real>::is_signed && (std::numeric_limits<Real>::digits > std::numeric_limits<Int>::digits))
       {
          m = (std::numeric_limits<Int>::min)();
          r = m.template convert_to<Int>();
@@ -556,48 +559,48 @@ void test_integer_ops(const boost::mpl::int_<boost::multiprecision::number_kind_
    Real a(20);
    Real b(7);
    Real c(5);
-   BOOST_CHECK_EQUAL(a % b ,  20 % 7);
-   BOOST_CHECK_EQUAL(a % 7 ,  20 % 7);
-   BOOST_CHECK_EQUAL(a % 7u ,  20 % 7);
-   BOOST_CHECK_EQUAL(a % a ,  0);
-   BOOST_CHECK_EQUAL(a % c ,  0);
-   BOOST_CHECK_EQUAL(a % 5 ,  0);
+   BOOST_CHECK_EQUAL(a % b, 20 % 7);
+   BOOST_CHECK_EQUAL(a % 7, 20 % 7);
+   BOOST_CHECK_EQUAL(a % 7u, 20 % 7);
+   BOOST_CHECK_EQUAL(a % a, 0);
+   BOOST_CHECK_EQUAL(a % c, 0);
+   BOOST_CHECK_EQUAL(a % 5, 0);
    a = a % (b + 0);
-   BOOST_CHECK_EQUAL(a ,  20 % 7);
+   BOOST_CHECK_EQUAL(a, 20 % 7);
    a = 20;
    c = (a + 2) % (a - 1);
-   BOOST_CHECK_EQUAL(c ,  22 % 19);
+   BOOST_CHECK_EQUAL(c, 22 % 19);
    c = 5;
    a = b % (a - 15);
-   BOOST_CHECK_EQUAL(a ,  7 % 5);
+   BOOST_CHECK_EQUAL(a, 7 % 5);
    a = 20;
 
    a = 20;
    a %= 7;
-   BOOST_CHECK_EQUAL(a ,  20 % 7);
+   BOOST_CHECK_EQUAL(a, 20 % 7);
 #ifndef BOOST_NO_LONG_LONG
    a = 20;
    a %= 7uLL;
-   BOOST_CHECK_EQUAL(a ,  20 % 7);
+   BOOST_CHECK_EQUAL(a, 20 % 7);
 #endif
    a = 20;
    ++a;
-   BOOST_CHECK_EQUAL(a ,  21);
+   BOOST_CHECK_EQUAL(a, 21);
    --a;
-   BOOST_CHECK_EQUAL(a ,  20);
-   BOOST_CHECK_EQUAL(a++ ,  20);
-   BOOST_CHECK_EQUAL(a ,  21);
-   BOOST_CHECK_EQUAL(a-- ,  21);
-   BOOST_CHECK_EQUAL(a ,  20);
+   BOOST_CHECK_EQUAL(a, 20);
+   BOOST_CHECK_EQUAL(a++, 20);
+   BOOST_CHECK_EQUAL(a, 21);
+   BOOST_CHECK_EQUAL(a--, 21);
+   BOOST_CHECK_EQUAL(a, 20);
    a = 2000;
    a <<= 20;
-   BOOST_CHECK_EQUAL(a ,  2000L << 20);
+   BOOST_CHECK_EQUAL(a, 2000L << 20);
    a >>= 20;
-   BOOST_CHECK_EQUAL(a ,  2000);
+   BOOST_CHECK_EQUAL(a, 2000);
    a <<= 20u;
-   BOOST_CHECK_EQUAL(a ,  2000L << 20);
+   BOOST_CHECK_EQUAL(a, 2000L << 20);
    a >>= 20u;
-   BOOST_CHECK_EQUAL(a ,  2000);
+   BOOST_CHECK_EQUAL(a, 2000);
 #ifndef BOOST_NO_EXCEPTIONS
    BOOST_CHECK_THROW(a <<= -20, std::out_of_range);
    BOOST_CHECK_THROW(a >>= -20, std::out_of_range);
@@ -605,7 +608,7 @@ void test_integer_ops(const boost::mpl::int_<boost::multiprecision::number_kind_
    BOOST_CHECK_THROW(Real(a >> -20), std::out_of_range);
 #endif
 #ifndef BOOST_NO_LONG_LONG
-   if(sizeof(long long) > sizeof(std::size_t))
+   if (sizeof(long long) > sizeof(std::size_t))
    {
       // extreme values should trigger an exception:
 #ifndef BOOST_NO_EXCEPTIONS
@@ -619,10 +622,10 @@ void test_integer_ops(const boost::mpl::int_<boost::multiprecision::number_kind_
       // Unless they fit within range:
       a = 2000L;
       a <<= 20uLL;
-      BOOST_CHECK_EQUAL(a,  (2000L << 20));
+      BOOST_CHECK_EQUAL(a, (2000L << 20));
       a = 2000;
       a <<= 20LL;
-      BOOST_CHECK_EQUAL(a,  (2000L << 20));
+      BOOST_CHECK_EQUAL(a, (2000L << 20));
 
 #ifndef BOOST_NO_EXCEPTIONS
       BOOST_CHECK_THROW(Real(a >> (1uLL << (sizeof(long long) * CHAR_BIT - 2))), std::out_of_range);
@@ -634,52 +637,52 @@ void test_integer_ops(const boost::mpl::int_<boost::multiprecision::number_kind_
 #endif
       // Unless they fit within range:
       a = 2000L;
-      BOOST_CHECK_EQUAL(Real(a << 20uLL) ,  (2000L << 20));
+      BOOST_CHECK_EQUAL(Real(a << 20uLL), (2000L << 20));
       a = 2000;
-      BOOST_CHECK_EQUAL(Real(a << 20LL)  ,  (2000L << 20));
+      BOOST_CHECK_EQUAL(Real(a << 20LL), (2000L << 20));
    }
 #endif
    a = 20;
    b = a << 20;
-   BOOST_CHECK_EQUAL(b ,  (20 << 20));
+   BOOST_CHECK_EQUAL(b, (20 << 20));
    b = a >> 2;
-   BOOST_CHECK_EQUAL(b ,  (20 >> 2));
+   BOOST_CHECK_EQUAL(b, (20 >> 2));
    b = (a + 2) << 10;
-   BOOST_CHECK_EQUAL(b ,  (22 << 10));
+   BOOST_CHECK_EQUAL(b, (22 << 10));
    b = (a + 3) >> 3;
-   BOOST_CHECK_EQUAL(b ,  (23 >> 3));
+   BOOST_CHECK_EQUAL(b, (23 >> 3));
    //
    // Bit fiddling:
    //
    int i = 1020304;
    int j = 56789123;
    int k = 4523187;
-   a = i;
-   b = j;
-   c = a;
+   a     = i;
+   b     = j;
+   c     = a;
    c &= b;
-   BOOST_CHECK_EQUAL(c ,  (i & j));
+   BOOST_CHECK_EQUAL(c, (i & j));
    c = a;
    c &= j;
-   BOOST_CHECK_EQUAL(c ,  (i & j));
+   BOOST_CHECK_EQUAL(c, (i & j));
    c = a;
    c &= a + b;
-   BOOST_CHECK_EQUAL(c ,  (i & (i + j)));
-   BOOST_CHECK_EQUAL((a & b) ,  (i & j));
+   BOOST_CHECK_EQUAL(c, (i & (i + j)));
+   BOOST_CHECK_EQUAL((a & b), (i & j));
    c = k;
    a = a & (b + k);
-   BOOST_CHECK_EQUAL(a ,  (i & (j + k)));
+   BOOST_CHECK_EQUAL(a, (i & (j + k)));
    a = i;
    a = (b + k) & a;
-   BOOST_CHECK_EQUAL(a ,  (i & (j + k)));
+   BOOST_CHECK_EQUAL(a, (i & (j + k)));
    a = i;
    c = a & b & k;
-   BOOST_CHECK_EQUAL(c ,  (i&j&k));
+   BOOST_CHECK_EQUAL(c, (i & j & k));
    c = a;
-   c &= (c+b);
-   BOOST_CHECK_EQUAL(c ,  (i & (i+j)));
+   c &= (c + b);
+   BOOST_CHECK_EQUAL(c, (i & (i + j)));
    c = a & (b | 1);
-   BOOST_CHECK_EQUAL(c ,  (i & (j | 1)));
+   BOOST_CHECK_EQUAL(c, (i & (j | 1)));
 
    test_complement<Real>(a, b, c, typename is_twos_complement_integer<Real>::type());
 
@@ -687,55 +690,55 @@ void test_integer_ops(const boost::mpl::int_<boost::multiprecision::number_kind_
    b = j;
    c = a;
    c |= b;
-   BOOST_CHECK_EQUAL(c ,  (i | j));
+   BOOST_CHECK_EQUAL(c, (i | j));
    c = a;
    c |= j;
-   BOOST_CHECK_EQUAL(c ,  (i | j));
+   BOOST_CHECK_EQUAL(c, (i | j));
    c = a;
    c |= a + b;
-   BOOST_CHECK_EQUAL(c ,  (i | (i + j)));
-   BOOST_CHECK_EQUAL((a | b) ,  (i | j));
+   BOOST_CHECK_EQUAL(c, (i | (i + j)));
+   BOOST_CHECK_EQUAL((a | b), (i | j));
    c = k;
    a = a | (b + k);
-   BOOST_CHECK_EQUAL(a ,  (i | (j + k)));
+   BOOST_CHECK_EQUAL(a, (i | (j + k)));
    a = i;
    a = (b + k) | a;
-   BOOST_CHECK_EQUAL(a ,  (i | (j + k)));
+   BOOST_CHECK_EQUAL(a, (i | (j + k)));
    a = i;
    c = a | b | k;
-   BOOST_CHECK_EQUAL(c ,  (i|j|k));
+   BOOST_CHECK_EQUAL(c, (i | j | k));
    c = a;
    c |= (c + b);
-   BOOST_CHECK_EQUAL(c ,  (i | (i+j)));
+   BOOST_CHECK_EQUAL(c, (i | (i + j)));
    c = a | (b | 1);
-   BOOST_CHECK_EQUAL(c ,  (i | (j | 1)));
+   BOOST_CHECK_EQUAL(c, (i | (j | 1)));
 
    a = i;
    b = j;
    c = a;
    c ^= b;
-   BOOST_CHECK_EQUAL(c ,  (i ^ j));
+   BOOST_CHECK_EQUAL(c, (i ^ j));
    c = a;
    c ^= j;
-   BOOST_CHECK_EQUAL(c ,  (i ^ j));
+   BOOST_CHECK_EQUAL(c, (i ^ j));
    c = a;
    c ^= a + b;
-   BOOST_CHECK_EQUAL(c ,  (i ^ (i + j)));
-   BOOST_CHECK_EQUAL((a ^ b) ,  (i ^ j));
+   BOOST_CHECK_EQUAL(c, (i ^ (i + j)));
+   BOOST_CHECK_EQUAL((a ^ b), (i ^ j));
    c = k;
    a = a ^ (b + k);
-   BOOST_CHECK_EQUAL(a ,  (i ^ (j + k)));
+   BOOST_CHECK_EQUAL(a, (i ^ (j + k)));
    a = i;
    a = (b + k) ^ a;
-   BOOST_CHECK_EQUAL(a ,  (i ^ (j + k)));
+   BOOST_CHECK_EQUAL(a, (i ^ (j + k)));
    a = i;
    c = a ^ b ^ k;
-   BOOST_CHECK_EQUAL(c ,  (i^j^k));
+   BOOST_CHECK_EQUAL(c, (i ^ j ^ k));
    c = a;
    c ^= (c + b);
-   BOOST_CHECK_EQUAL(c ,  (i ^ (i+j)));
+   BOOST_CHECK_EQUAL(c, (i ^ (i + j)));
    c = a ^ (b | 1);
-   BOOST_CHECK_EQUAL(c ,  (i ^ (j | 1)));
+   BOOST_CHECK_EQUAL(c, (i ^ (j | 1)));
 
    a = i;
    b = j;
@@ -745,16 +748,16 @@ void test_integer_ops(const boost::mpl::int_<boost::multiprecision::number_kind_
    //
    a = 400;
    b = 45;
-   BOOST_CHECK_EQUAL(gcd(a, b) ,  boost::integer::gcd(400, 45));
-   BOOST_CHECK_EQUAL(lcm(a, b) ,  boost::integer::lcm(400, 45));
-   BOOST_CHECK_EQUAL(gcd(a, 45) ,  boost::integer::gcd(400, 45));
-   BOOST_CHECK_EQUAL(lcm(a, 45) ,  boost::integer::lcm(400, 45));
-   BOOST_CHECK_EQUAL(gcd(a, 45u) ,  boost::integer::gcd(400, 45));
-   BOOST_CHECK_EQUAL(lcm(a, 45u) ,  boost::integer::lcm(400, 45));
-   BOOST_CHECK_EQUAL(gcd(400, b) ,  boost::integer::gcd(400, 45));
-   BOOST_CHECK_EQUAL(lcm(400, b) ,  boost::integer::lcm(400, 45));
-   BOOST_CHECK_EQUAL(gcd(400u, b) ,  boost::integer::gcd(400, 45));
-   BOOST_CHECK_EQUAL(lcm(400u, b) ,  boost::integer::lcm(400, 45));
+   BOOST_CHECK_EQUAL(gcd(a, b), boost::integer::gcd(400, 45));
+   BOOST_CHECK_EQUAL(lcm(a, b), boost::integer::lcm(400, 45));
+   BOOST_CHECK_EQUAL(gcd(a, 45), boost::integer::gcd(400, 45));
+   BOOST_CHECK_EQUAL(lcm(a, 45), boost::integer::lcm(400, 45));
+   BOOST_CHECK_EQUAL(gcd(a, 45u), boost::integer::gcd(400, 45));
+   BOOST_CHECK_EQUAL(lcm(a, 45u), boost::integer::lcm(400, 45));
+   BOOST_CHECK_EQUAL(gcd(400, b), boost::integer::gcd(400, 45));
+   BOOST_CHECK_EQUAL(lcm(400, b), boost::integer::lcm(400, 45));
+   BOOST_CHECK_EQUAL(gcd(400u, b), boost::integer::gcd(400, 45));
+   BOOST_CHECK_EQUAL(lcm(400u, b), boost::integer::lcm(400, 45));
 
    //
    // Conditionals involving 2 arg functions:
@@ -763,63 +766,63 @@ void test_integer_ops(const boost::mpl::int_<boost::multiprecision::number_kind_
 
    Real r;
    divide_qr(a, b, c, r);
-   BOOST_CHECK_EQUAL(c ,  a / b);
-   BOOST_CHECK_EQUAL(r ,  a % b);
+   BOOST_CHECK_EQUAL(c, a / b);
+   BOOST_CHECK_EQUAL(r, a % b);
    divide_qr(a + 0, b, c, r);
-   BOOST_CHECK_EQUAL(c ,  a / b);
-   BOOST_CHECK_EQUAL(r ,  a % b);
-   divide_qr(a, b+0, c, r);
-   BOOST_CHECK_EQUAL(c ,  a / b);
-   BOOST_CHECK_EQUAL(r ,  a % b);
-   divide_qr(a+0, b+0, c, r);
-   BOOST_CHECK_EQUAL(c ,  a / b);
-   BOOST_CHECK_EQUAL(r ,  a % b);
-   BOOST_CHECK_EQUAL(integer_modulus(a, 57) ,  a % 57);
-   for(i = 0; i < 20; ++i)
+   BOOST_CHECK_EQUAL(c, a / b);
+   BOOST_CHECK_EQUAL(r, a % b);
+   divide_qr(a, b + 0, c, r);
+   BOOST_CHECK_EQUAL(c, a / b);
+   BOOST_CHECK_EQUAL(r, a % b);
+   divide_qr(a + 0, b + 0, c, r);
+   BOOST_CHECK_EQUAL(c, a / b);
+   BOOST_CHECK_EQUAL(r, a % b);
+   BOOST_CHECK_EQUAL(integer_modulus(a, 57), a % 57);
+   for (i = 0; i < 20; ++i)
    {
-      if(std::numeric_limits<Real>::is_specialized && (!std::numeric_limits<Real>::is_bounded || ((int)i * 17 < std::numeric_limits<Real>::digits)))
+      if (std::numeric_limits<Real>::is_specialized && (!std::numeric_limits<Real>::is_bounded || ((int)i * 17 < std::numeric_limits<Real>::digits)))
       {
-         BOOST_CHECK_EQUAL(lsb(Real(1) << (i * 17)),  static_cast<unsigned>(i * 17));
+         BOOST_CHECK_EQUAL(lsb(Real(1) << (i * 17)), static_cast<unsigned>(i * 17));
          BOOST_CHECK_EQUAL(msb(Real(1) << (i * 17)), static_cast<unsigned>(i * 17));
          BOOST_CHECK(bit_test(Real(1) << (i * 17), i * 17));
          BOOST_CHECK(!bit_test(Real(1) << (i * 17), i * 17 + 1));
-         if(i)
+         if (i)
          {
             BOOST_CHECK(!bit_test(Real(1) << (i * 17), i * 17 - 1));
          }
          Real zero(0);
          BOOST_CHECK(bit_test(bit_set(zero, i * 17), i * 17));
          zero = 0;
-         BOOST_CHECK_EQUAL(bit_flip(zero, i*17) ,  Real(1) << i * 17);
+         BOOST_CHECK_EQUAL(bit_flip(zero, i * 17), Real(1) << i * 17);
          zero = Real(1) << i * 17;
-         BOOST_CHECK_EQUAL(bit_flip(zero, i * 17) ,  0);
+         BOOST_CHECK_EQUAL(bit_flip(zero, i * 17), 0);
          zero = Real(1) << i * 17;
-         BOOST_CHECK_EQUAL(bit_unset(zero, i * 17) ,  0);
+         BOOST_CHECK_EQUAL(bit_unset(zero, i * 17), 0);
       }
    }
    //
    // pow, powm:
    //
-   BOOST_CHECK_EQUAL(pow(Real(3), 4u) ,  81);
-   BOOST_CHECK_EQUAL(pow(Real(3) + Real(0), 4u) ,  81);
-   BOOST_CHECK_EQUAL(powm(Real(3), Real(4), Real(13)) ,  81 % 13);
-   BOOST_CHECK_EQUAL(powm(Real(3), Real(4), 13) ,  81 % 13);
-   BOOST_CHECK_EQUAL(powm(Real(3), Real(4), Real(13) + 0) ,  81 % 13);
-   BOOST_CHECK_EQUAL(powm(Real(3), Real(4) + 0, Real(13)) ,  81 % 13);
-   BOOST_CHECK_EQUAL(powm(Real(3), Real(4) + 0, 13) ,  81 % 13);
-   BOOST_CHECK_EQUAL(powm(Real(3), Real(4) + 0, Real(13) + 0) ,  81 % 13);
-   BOOST_CHECK_EQUAL(powm(Real(3), 4 + 0, Real(13)) ,  81 % 13);
-   BOOST_CHECK_EQUAL(powm(Real(3), 4 + 0, 13) ,  81 % 13);
-   BOOST_CHECK_EQUAL(powm(Real(3), 4 + 0, Real(13) + 0) ,  81 % 13);
-   BOOST_CHECK_EQUAL(powm(Real(3) + 0, Real(4), Real(13)) ,  81 % 13);
-   BOOST_CHECK_EQUAL(powm(Real(3) + 0, Real(4), 13) ,  81 % 13);
-   BOOST_CHECK_EQUAL(powm(Real(3) + 0, Real(4), Real(13) + 0) ,  81 % 13);
-   BOOST_CHECK_EQUAL(powm(Real(3) + 0, Real(4) + 0, Real(13)) ,  81 % 13);
-   BOOST_CHECK_EQUAL(powm(Real(3) + 0, Real(4) + 0, 13) ,  81 % 13);
-   BOOST_CHECK_EQUAL(powm(Real(3) + 0, Real(4) + 0, Real(13) + 0) ,  81 % 13);
-   BOOST_CHECK_EQUAL(powm(Real(3) + 0, 4 + 0, Real(13)) ,  81 % 13);
-   BOOST_CHECK_EQUAL(powm(Real(3) + 0, 4 + 0, 13) ,  81 % 13);
-   BOOST_CHECK_EQUAL(powm(Real(3) + 0, 4 + 0, Real(13) + 0) ,  81 % 13);
+   BOOST_CHECK_EQUAL(pow(Real(3), 4u), 81);
+   BOOST_CHECK_EQUAL(pow(Real(3) + Real(0), 4u), 81);
+   BOOST_CHECK_EQUAL(powm(Real(3), Real(4), Real(13)), 81 % 13);
+   BOOST_CHECK_EQUAL(powm(Real(3), Real(4), 13), 81 % 13);
+   BOOST_CHECK_EQUAL(powm(Real(3), Real(4), Real(13) + 0), 81 % 13);
+   BOOST_CHECK_EQUAL(powm(Real(3), Real(4) + 0, Real(13)), 81 % 13);
+   BOOST_CHECK_EQUAL(powm(Real(3), Real(4) + 0, 13), 81 % 13);
+   BOOST_CHECK_EQUAL(powm(Real(3), Real(4) + 0, Real(13) + 0), 81 % 13);
+   BOOST_CHECK_EQUAL(powm(Real(3), 4 + 0, Real(13)), 81 % 13);
+   BOOST_CHECK_EQUAL(powm(Real(3), 4 + 0, 13), 81 % 13);
+   BOOST_CHECK_EQUAL(powm(Real(3), 4 + 0, Real(13) + 0), 81 % 13);
+   BOOST_CHECK_EQUAL(powm(Real(3) + 0, Real(4), Real(13)), 81 % 13);
+   BOOST_CHECK_EQUAL(powm(Real(3) + 0, Real(4), 13), 81 % 13);
+   BOOST_CHECK_EQUAL(powm(Real(3) + 0, Real(4), Real(13) + 0), 81 % 13);
+   BOOST_CHECK_EQUAL(powm(Real(3) + 0, Real(4) + 0, Real(13)), 81 % 13);
+   BOOST_CHECK_EQUAL(powm(Real(3) + 0, Real(4) + 0, 13), 81 % 13);
+   BOOST_CHECK_EQUAL(powm(Real(3) + 0, Real(4) + 0, Real(13) + 0), 81 % 13);
+   BOOST_CHECK_EQUAL(powm(Real(3) + 0, 4 + 0, Real(13)), 81 % 13);
+   BOOST_CHECK_EQUAL(powm(Real(3) + 0, 4 + 0, 13), 81 % 13);
+   BOOST_CHECK_EQUAL(powm(Real(3) + 0, 4 + 0, Real(13) + 0), 81 % 13);
    //
    // Conditionals involving 3 arg functions:
    //
@@ -839,32 +842,31 @@ void test_integer_ops(const boost::mpl::int_<boost::multiprecision::number_kind_
    a = 20;
    b = 7;
    c = 20 % b;
-   BOOST_CHECK_EQUAL(c ,  (20 % 7));
+   BOOST_CHECK_EQUAL(c, (20 % 7));
    c = 20 % (b + 0);
-   BOOST_CHECK_EQUAL(c ,  (20 % 7));
+   BOOST_CHECK_EQUAL(c, (20 % 7));
    c = a & 10;
-   BOOST_CHECK_EQUAL(c ,  (20 & 10));
+   BOOST_CHECK_EQUAL(c, (20 & 10));
    c = 10 & a;
-   BOOST_CHECK_EQUAL(c ,  (20 & 10));
+   BOOST_CHECK_EQUAL(c, (20 & 10));
    c = (a + 0) & (b + 0);
-   BOOST_CHECK_EQUAL(c ,  (20 & 7));
+   BOOST_CHECK_EQUAL(c, (20 & 7));
    c = 10 & (a + 0);
-   BOOST_CHECK_EQUAL(c ,  (20 & 10));
+   BOOST_CHECK_EQUAL(c, (20 & 10));
    c = 10 | a;
-   BOOST_CHECK_EQUAL(c ,  (20 | 10));
+   BOOST_CHECK_EQUAL(c, (20 | 10));
    c = (a + 0) | (b + 0);
    BOOST_CHECK(c == (20 | 7))
    c = 20 | (b + 0);
-   BOOST_CHECK_EQUAL(c ,  (20 | 7));
+   BOOST_CHECK_EQUAL(c, (20 | 7));
    c = a ^ 7;
-   BOOST_CHECK_EQUAL(c ,  (20 ^ 7));
+   BOOST_CHECK_EQUAL(c, (20 ^ 7));
    c = 20 ^ b;
-   BOOST_CHECK_EQUAL(c ,  (20 ^ 7));
+   BOOST_CHECK_EQUAL(c, (20 ^ 7));
    c = (a + 0) ^ (b + 0);
-   BOOST_CHECK_EQUAL(c ,  (20 ^ 7));
+   BOOST_CHECK_EQUAL(c, (20 ^ 7));
    c = 20 ^ (b + 0);
-   BOOST_CHECK_EQUAL(c ,  (20 ^ 7));
-
+   BOOST_CHECK_EQUAL(c, (20 ^ 7));
 
    //
    // Round tripping of built in integers:
@@ -882,12 +884,12 @@ void test_integer_ops(const boost::mpl::int_<boost::multiprecision::number_kind_
 }
 
 template <class Real, class T>
-void test_float_funcs(const T&){}
+void test_float_funcs(const T&) {}
 
 template <class Real>
 void test_float_funcs(const boost::mpl::true_&)
 {
-   if(boost::multiprecision::is_interval_number<Real>::value)
+   if (boost::multiprecision::is_interval_number<Real>::value)
       return;
    //
    // Test variable reuse in function calls, see https://svn.boost.org/trac/boost/ticket/8326
@@ -927,8 +929,8 @@ void test_float_funcs(const boost::mpl::true_&)
    BOOST_CHECK_EQUAL(a, 0.5);
 
    Real tol = std::numeric_limits<Real>::epsilon() * 3;
-   a = 4;
-   a = sqrt(a);
+   a        = 4;
+   a        = sqrt(a);
    BOOST_CHECK_CLOSE_FRACTION(a, 2, tol);
    a = 3;
    a = exp(a);
@@ -1022,7 +1024,7 @@ void test_float_funcs(const boost::mpl::true_&)
    BOOST_CHECK_EQUAL(b > 0, a > 0);
    BOOST_CHECK_EQUAL(c > 0, a > 0);
 
-   if(std::numeric_limits<Real>::has_infinity)
+   if (std::numeric_limits<Real>::has_infinity)
    {
       a = std::numeric_limits<Real>::infinity();
       b = modf(a, &c);
@@ -1033,7 +1035,7 @@ void test_float_funcs(const boost::mpl::true_&)
       BOOST_CHECK_EQUAL(a, c);
       BOOST_CHECK_EQUAL(b, 0);
    }
-   if(std::numeric_limits<Real>::has_quiet_NaN)
+   if (std::numeric_limits<Real>::has_quiet_NaN)
    {
       a = std::numeric_limits<Real>::quiet_NaN();
       b = modf(a, &c);
@@ -1122,67 +1124,67 @@ void compare_NaNs(const T& a, const U& b)
 }
 
 template <class Real, class T>
-void test_float_ops(const T&){}
+void test_float_ops(const T&) {}
 
 template <class Real>
 void test_float_ops(const boost::mpl::int_<boost::multiprecision::number_kind_floating_point>&)
 {
-   BOOST_CHECK_EQUAL(abs(Real(2)) ,  2);
-   BOOST_CHECK_EQUAL(abs(Real(-2)) ,  2);
-   BOOST_CHECK_EQUAL(fabs(Real(2)) ,  2);
-   BOOST_CHECK_EQUAL(fabs(Real(-2)) ,  2);
-   BOOST_CHECK_EQUAL(floor(Real(5) / 2) ,  2);
-   BOOST_CHECK_EQUAL(ceil(Real(5) / 2) ,  3);
-   BOOST_CHECK_EQUAL(floor(Real(-5) / 2) ,  -3);
-   BOOST_CHECK_EQUAL(ceil(Real(-5) / 2) ,  -2);
-   BOOST_CHECK_EQUAL(trunc(Real(5) / 2) ,  2);
-   BOOST_CHECK_EQUAL(trunc(Real(-5) / 2) ,  -2);
+   BOOST_CHECK_EQUAL(abs(Real(2)), 2);
+   BOOST_CHECK_EQUAL(abs(Real(-2)), 2);
+   BOOST_CHECK_EQUAL(fabs(Real(2)), 2);
+   BOOST_CHECK_EQUAL(fabs(Real(-2)), 2);
+   BOOST_CHECK_EQUAL(floor(Real(5) / 2), 2);
+   BOOST_CHECK_EQUAL(ceil(Real(5) / 2), 3);
+   BOOST_CHECK_EQUAL(floor(Real(-5) / 2), -3);
+   BOOST_CHECK_EQUAL(ceil(Real(-5) / 2), -2);
+   BOOST_CHECK_EQUAL(trunc(Real(5) / 2), 2);
+   BOOST_CHECK_EQUAL(trunc(Real(-5) / 2), -2);
    //
    // ldexp and frexp, these pretty much have to be implemented by each backend:
    //
    typedef typename Real::backend_type::exponent_type e_type;
-   BOOST_CHECK_EQUAL(ldexp(Real(2), 5) ,  64);
-   BOOST_CHECK_EQUAL(ldexp(Real(2), -5) ,  Real(2) / 32);
-   Real v(512);
+   BOOST_CHECK_EQUAL(ldexp(Real(2), 5), 64);
+   BOOST_CHECK_EQUAL(ldexp(Real(2), -5), Real(2) / 32);
+   Real   v(512);
    e_type exponent;
-   Real r = frexp(v, &exponent);
-   BOOST_CHECK_EQUAL(r ,  0.5);
-   BOOST_CHECK_EQUAL(exponent ,  10);
-   BOOST_CHECK_EQUAL(v ,  512);
+   Real   r = frexp(v, &exponent);
+   BOOST_CHECK_EQUAL(r, 0.5);
+   BOOST_CHECK_EQUAL(exponent, 10);
+   BOOST_CHECK_EQUAL(v, 512);
    v = 1 / v;
    r = frexp(v, &exponent);
-   BOOST_CHECK_EQUAL(r ,  0.5);
-   BOOST_CHECK_EQUAL(exponent ,  -8);
-   BOOST_CHECK_EQUAL(ldexp(Real(2), e_type(5)) ,  64);
-   BOOST_CHECK_EQUAL(ldexp(Real(2), e_type(-5)) ,  Real(2) / 32);
+   BOOST_CHECK_EQUAL(r, 0.5);
+   BOOST_CHECK_EQUAL(exponent, -8);
+   BOOST_CHECK_EQUAL(ldexp(Real(2), e_type(5)), 64);
+   BOOST_CHECK_EQUAL(ldexp(Real(2), e_type(-5)), Real(2) / 32);
    v = 512;
    e_type exp2;
    r = frexp(v, &exp2);
-   BOOST_CHECK_EQUAL(r ,  0.5);
-   BOOST_CHECK_EQUAL(exp2 ,  10);
-   BOOST_CHECK_EQUAL(v ,  512);
+   BOOST_CHECK_EQUAL(r, 0.5);
+   BOOST_CHECK_EQUAL(exp2, 10);
+   BOOST_CHECK_EQUAL(v, 512);
    v = 1 / v;
    r = frexp(v, &exp2);
-   BOOST_CHECK_EQUAL(r ,  0.5);
-   BOOST_CHECK_EQUAL(exp2 ,  -8);
+   BOOST_CHECK_EQUAL(r, 0.5);
+   BOOST_CHECK_EQUAL(exp2, -8);
    //
    // scalbn and logb, these are the same as ldexp and frexp unless the radix is
    // something other than 2:
    //
-   if(std::numeric_limits<Real>::is_specialized && std::numeric_limits<Real>::radix)
+   if (std::numeric_limits<Real>::is_specialized && std::numeric_limits<Real>::radix)
    {
       BOOST_CHECK_EQUAL(scalbn(Real(2), 5), 2 * pow(double(std::numeric_limits<Real>::radix), 5));
       BOOST_CHECK_EQUAL(scalbn(Real(2), -5), Real(2) / pow(double(std::numeric_limits<Real>::radix), 5));
-      v = 512;
+      v        = 512;
       exponent = ilogb(v);
-      r = scalbn(v, -exponent);
+      r        = scalbn(v, -exponent);
       BOOST_CHECK(r >= 1);
       BOOST_CHECK(r < std::numeric_limits<Real>::radix);
       BOOST_CHECK_EQUAL(exponent, logb(v));
       BOOST_CHECK_EQUAL(v, scalbn(r, exponent));
-      v = 1 / v;
+      v        = 1 / v;
       exponent = ilogb(v);
-      r = scalbn(v, -exponent);
+      r        = scalbn(v, -exponent);
       BOOST_CHECK(r >= 1);
       BOOST_CHECK(r < std::numeric_limits<Real>::radix);
       BOOST_CHECK_EQUAL(exponent, logb(v));
@@ -1193,43 +1195,43 @@ void test_float_ops(const boost::mpl::int_<boost::multiprecision::number_kind_fl
    //
    v = 3.25;
    r = pow(v, 0);
-   BOOST_CHECK_EQUAL(r ,  1);
+   BOOST_CHECK_EQUAL(r, 1);
    r = pow(v, 1);
-   BOOST_CHECK_EQUAL(r ,  3.25);
+   BOOST_CHECK_EQUAL(r, 3.25);
    r = pow(v, 2);
-   BOOST_CHECK_EQUAL(r ,  boost::math::pow<2>(3.25));
+   BOOST_CHECK_EQUAL(r, boost::math::pow<2>(3.25));
    r = pow(v, 3);
-   BOOST_CHECK_EQUAL(r ,  boost::math::pow<3>(3.25));
+   BOOST_CHECK_EQUAL(r, boost::math::pow<3>(3.25));
    r = pow(v, 4);
-   BOOST_CHECK_EQUAL(r ,  boost::math::pow<4>(3.25));
+   BOOST_CHECK_EQUAL(r, boost::math::pow<4>(3.25));
    r = pow(v, 5);
-   BOOST_CHECK_EQUAL(r ,  boost::math::pow<5>(3.25));
+   BOOST_CHECK_EQUAL(r, boost::math::pow<5>(3.25));
    r = pow(v, 6);
-   BOOST_CHECK_EQUAL(r ,  boost::math::pow<6>(3.25));
+   BOOST_CHECK_EQUAL(r, boost::math::pow<6>(3.25));
    r = pow(v, 25);
-   BOOST_CHECK_EQUAL(r ,  boost::math::pow<25>(Real(3.25)));
+   BOOST_CHECK_EQUAL(r, boost::math::pow<25>(Real(3.25)));
 
 #ifndef BOOST_NO_EXCEPTIONS
    //
    // Things that are expected errors:
    //
    BOOST_CHECK_THROW(Real("3.14L"), std::runtime_error);
-   if(std::numeric_limits<Real>::is_specialized)
+   if (std::numeric_limits<Real>::is_specialized)
    {
-      if(std::numeric_limits<Real>::has_infinity)
+      if (std::numeric_limits<Real>::has_infinity)
       {
          BOOST_CHECK((boost::math::isinf)(Real(20) / 0u));
       }
       else
       {
-         BOOST_CHECK_THROW(Real(Real(20) / 0u), std::overflow_error);
+         BOOST_CHECK_THROW(r = Real(Real(20) / 0u), std::overflow_error);
       }
    }
 #endif
    //
    // Comparisons of NaN's should always fail:
    //
-   if(std::numeric_limits<Real>::has_quiet_NaN)
+   if (std::numeric_limits<Real>::has_quiet_NaN)
    {
       r = v = std::numeric_limits<Real>::quiet_NaN();
       compare_NaNs(r, v);
@@ -1247,7 +1249,7 @@ void test_float_ops(const boost::mpl::int_<boost::multiprecision::number_kind_fl
       //
       compare_NaNs(v, 0.5);
       compare_NaNs(0.5, v);
-      if(std::numeric_limits<double>::has_quiet_NaN)
+      if (std::numeric_limits<double>::has_quiet_NaN)
       {
          compare_NaNs(r, std::numeric_limits<double>::quiet_NaN());
          compare_NaNs(std::numeric_limits<double>::quiet_NaN(), r);
@@ -1257,7 +1259,7 @@ void test_float_ops(const boost::mpl::int_<boost::multiprecision::number_kind_fl
    //
    // Operations involving NaN's as one argument:
    //
-   if(std::numeric_limits<Real>::has_quiet_NaN)
+   if (std::numeric_limits<Real>::has_quiet_NaN)
    {
       v = 20.25;
       r = std::numeric_limits<Real>::quiet_NaN();
@@ -1289,7 +1291,7 @@ void test_float_ops(const boost::mpl::int_<boost::multiprecision::number_kind_fl
    //
    // Operations involving infinities as one argument:
    //
-   if(std::numeric_limits<Real>::has_infinity)
+   if (std::numeric_limits<Real>::has_infinity)
    {
       v = 20.25;
       r = std::numeric_limits<Real>::infinity();
@@ -1324,21 +1326,21 @@ void test_float_ops(const boost::mpl::int_<boost::multiprecision::number_kind_fl
    //
    // Operations that should produce NaN as a result:
    //
-   if(std::numeric_limits<Real>::has_quiet_NaN)
+   if (std::numeric_limits<Real>::has_quiet_NaN)
    {
-      v = r = 0;
+      v = r  = 0;
       Real t = v / r;
       BOOST_CHECK((boost::math::isnan)(t));
       v /= r;
       BOOST_CHECK((boost::math::isnan)(v));
       t = v / 0;
       BOOST_CHECK((boost::math::isnan)(v));
-      if(std::numeric_limits<Real>::has_infinity)
+      if (std::numeric_limits<Real>::has_infinity)
       {
          v = 0;
          r = std::numeric_limits<Real>::infinity();
          t = v * r;
-         if(!boost::multiprecision::is_interval_number<Real>::value)
+         if (!boost::multiprecision::is_interval_number<Real>::value)
          {
             BOOST_CHECK((boost::math::isnan)(t));
             t = r * 0;
@@ -1357,20 +1359,18 @@ template <class T>
 struct lexical_cast_target_type
 {
    typedef typename boost::mpl::if_<
-      boost::is_signed<T>,
-      boost::intmax_t,
-      typename boost::mpl::if_<
-         boost::is_unsigned<T>,
-         boost::uintmax_t,
-         T
-      >::type
-   >::type type;
+       boost::is_signed<T>,
+       boost::intmax_t,
+       typename boost::mpl::if_<
+           boost::is_unsigned<T>,
+           boost::uintmax_t,
+           T>::type>::type type;
 };
 
 template <class Real, class Num>
 void test_negative_mixed_minmax(boost::mpl::true_ const&)
 {
-   if(!std::numeric_limits<Real>::is_bounded || (std::numeric_limits<Real>::digits >= std::numeric_limits<Num>::digits))
+   if (!std::numeric_limits<Real>::is_bounded || (std::numeric_limits<Real>::digits >= std::numeric_limits<Num>::digits))
    {
       Real mx1((std::numeric_limits<Num>::max)() - 1);
       ++mx1;
@@ -1381,7 +1381,7 @@ void test_negative_mixed_minmax(boost::mpl::true_ const&)
       mx2 = (std::numeric_limits<Num>::max)();
       BOOST_CHECK_EQUAL(mx1, mx2);
 
-      if(!std::numeric_limits<Real>::is_bounded || (std::numeric_limits<Real>::digits > std::numeric_limits<Num>::digits))
+      if (!std::numeric_limits<Real>::is_bounded || (std::numeric_limits<Real>::digits > std::numeric_limits<Num>::digits))
       {
          Real mx3((std::numeric_limits<Num>::min)() + 1);
          --mx3;
@@ -1408,17 +1408,18 @@ void test_negative_mixed_numeric_limits(boost::mpl::true_ const&)
 #else
    Num tol = 0;
 #endif
-   static const int left_shift = std::numeric_limits<Num>::digits - 1;
-   Num n1 = -static_cast<Num>(1uLL << ((left_shift < 63) && (left_shift > 0) ? left_shift : 10));
-   Num n2 = -1;
-   Num n3 = 0;
-   Num n4 = -20;
-   std::ios_base::fmtflags f = boost::is_floating_point<Num>::value ? std::ios_base::scientific : std::ios_base::fmtflags(0);
-   int digits_to_print = boost::is_floating_point<Num>::value && std::numeric_limits<Num>::is_specialized
-   ? std::numeric_limits<Num>::digits10 + 5 : 0;
-   if(std::numeric_limits<target_type>::digits <= std::numeric_limits<Real>::digits)
+   static const int        left_shift      = std::numeric_limits<Num>::digits - 1;
+   Num                     n1              = -static_cast<Num>(1uLL << ((left_shift < 63) && (left_shift > 0) ? left_shift : 10));
+   Num                     n2              = -1;
+   Num                     n3              = 0;
+   Num                     n4              = -20;
+   std::ios_base::fmtflags f               = boost::is_floating_point<Num>::value ? std::ios_base::scientific : std::ios_base::fmtflags(0);
+   int                     digits_to_print = boost::is_floating_point<Num>::value && std::numeric_limits<Num>::is_specialized
+                             ? std::numeric_limits<Num>::digits10 + 5
+                             : 0;
+   if (std::numeric_limits<target_type>::digits <= std::numeric_limits<Real>::digits)
    {
-   BOOST_CHECK_CLOSE(n1, checked_lexical_cast<target_type>(Real(n1).str(digits_to_print, f)), tol);
+      BOOST_CHECK_CLOSE(n1, checked_lexical_cast<target_type>(Real(n1).str(digits_to_print, f)), tol);
    }
    BOOST_CHECK_CLOSE(n2, checked_lexical_cast<target_type>(Real(n2).str(digits_to_print, f)), 0);
    BOOST_CHECK_CLOSE(n3, checked_lexical_cast<target_type>(Real(n3).str(digits_to_print, f)), 0);
@@ -1426,28 +1427,26 @@ void test_negative_mixed_numeric_limits(boost::mpl::true_ const&)
 }
 
 template <class Real, class Num>
-void test_negative_mixed_numeric_limits(boost::mpl::false_ const&){}
+void test_negative_mixed_numeric_limits(boost::mpl::false_ const&) {}
 
 template <class Real, class Num>
 void test_negative_mixed(boost::mpl::true_ const&)
 {
    typedef typename boost::mpl::if_<
-         boost::is_convertible<Num, Real>,
-         typename boost::mpl::if_c<boost::is_integral<Num>::value && (sizeof(Num) < sizeof(int)), int, Num>::type,
-         Real
-      >::type cast_type;
+       boost::is_convertible<Num, Real>,
+       typename boost::mpl::if_c<boost::is_integral<Num>::value && (sizeof(Num) < sizeof(int)), int, Num>::type,
+       Real>::type cast_type;
    typedef typename boost::mpl::if_<
-         boost::is_convertible<Num, Real>,
-         Num,
-         Real
-      >::type simple_cast_type;
-   std::cout << "Testing mixed arithmetic with type: " << typeid(Real).name()  << " and " << typeid(Num).name() << std::endl;
+       boost::is_convertible<Num, Real>,
+       Num,
+       Real>::type simple_cast_type;
+   std::cout << "Testing mixed arithmetic with type: " << typeid(Real).name() << " and " << typeid(Num).name() << std::endl;
    static const int left_shift = std::numeric_limits<Num>::digits - 1;
-   Num n1 = -static_cast<Num>(1uLL << ((left_shift < 63) && (left_shift > 0) ? left_shift : 10));
-   Num n2 = -1;
-   Num n3 = 0;
-   Num n4 = -20;
-   Num n5 = -8;
+   Num              n1         = -static_cast<Num>(1uLL << ((left_shift < 63) && (left_shift > 0) ? left_shift : 10));
+   Num              n2         = -1;
+   Num              n3         = 0;
+   Num              n4         = -20;
+   Num              n5         = -8;
 
    test_comparisons<Real>(n1, n2, boost::is_convertible<Num, Real>());
    test_comparisons<Real>(n1, n3, boost::is_convertible<Num, Real>());
@@ -1457,26 +1456,26 @@ void test_negative_mixed(boost::mpl::true_ const&)
    test_comparisons<Real>(n3, n3, boost::is_convertible<Num, Real>());
 
    // Default construct:
-   BOOST_CHECK_EQUAL(Real(n1) ,  static_cast<cast_type>(n1));
-   BOOST_CHECK_EQUAL(Real(n2) ,  static_cast<cast_type>(n2));
-   BOOST_CHECK_EQUAL(Real(n3) ,  static_cast<cast_type>(n3));
-   BOOST_CHECK_EQUAL(Real(n4) ,  static_cast<cast_type>(n4));
-   BOOST_CHECK_EQUAL(static_cast<cast_type>(n1) ,  Real(n1));
-   BOOST_CHECK_EQUAL(static_cast<cast_type>(n2) ,  Real(n2));
-   BOOST_CHECK_EQUAL(static_cast<cast_type>(n3) ,  Real(n3));
-   BOOST_CHECK_EQUAL(static_cast<cast_type>(n4) ,  Real(n4));
-   BOOST_CHECK_EQUAL(Real(n1).template convert_to<Num>() ,  n1);
-   BOOST_CHECK_EQUAL(Real(n2).template convert_to<Num>() ,  n2);
-   BOOST_CHECK_EQUAL(Real(n3).template convert_to<Num>() ,  n3);
-   BOOST_CHECK_EQUAL(Real(n4).template convert_to<Num>() ,  n4);
+   BOOST_CHECK_EQUAL(Real(n1), static_cast<cast_type>(n1));
+   BOOST_CHECK_EQUAL(Real(n2), static_cast<cast_type>(n2));
+   BOOST_CHECK_EQUAL(Real(n3), static_cast<cast_type>(n3));
+   BOOST_CHECK_EQUAL(Real(n4), static_cast<cast_type>(n4));
+   BOOST_CHECK_EQUAL(static_cast<cast_type>(n1), Real(n1));
+   BOOST_CHECK_EQUAL(static_cast<cast_type>(n2), Real(n2));
+   BOOST_CHECK_EQUAL(static_cast<cast_type>(n3), Real(n3));
+   BOOST_CHECK_EQUAL(static_cast<cast_type>(n4), Real(n4));
+   BOOST_CHECK_EQUAL(Real(n1).template convert_to<Num>(), n1);
+   BOOST_CHECK_EQUAL(Real(n2).template convert_to<Num>(), n2);
+   BOOST_CHECK_EQUAL(Real(n3).template convert_to<Num>(), n3);
+   BOOST_CHECK_EQUAL(Real(n4).template convert_to<Num>(), n4);
 #ifndef BOOST_MP_NO_CXX11_EXPLICIT_CONVERSION_OPERATORS
-   BOOST_CHECK_EQUAL(static_cast<Num>(Real(n1)) ,  n1);
-   BOOST_CHECK_EQUAL(static_cast<Num>(Real(n2)) ,  n2);
-   BOOST_CHECK_EQUAL(static_cast<Num>(Real(n3)) ,  n3);
-   BOOST_CHECK_EQUAL(static_cast<Num>(Real(n4)) ,  n4);
+   BOOST_CHECK_EQUAL(static_cast<Num>(Real(n1)), n1);
+   BOOST_CHECK_EQUAL(static_cast<Num>(Real(n2)), n2);
+   BOOST_CHECK_EQUAL(static_cast<Num>(Real(n3)), n3);
+   BOOST_CHECK_EQUAL(static_cast<Num>(Real(n4)), n4);
 #endif
    // Conversions when source is an expression template:
-   BOOST_CHECK_EQUAL((Real(n1) + 0).template convert_to<Num>() ,  n1);
+   BOOST_CHECK_EQUAL((Real(n1) + 0).template convert_to<Num>(), n1);
    BOOST_CHECK_EQUAL((Real(n2) + 0).template convert_to<Num>(), n2);
    BOOST_CHECK_EQUAL((Real(n3) + 0).template convert_to<Num>(), n3);
    BOOST_CHECK_EQUAL((Real(n4) + 0).template convert_to<Num>(), n4);
@@ -1491,218 +1490,218 @@ void test_negative_mixed(boost::mpl::true_ const&)
    Real r(0);
    BOOST_CHECK(r != static_cast<cast_type>(n1));
    r = static_cast<simple_cast_type>(n1);
-   BOOST_CHECK_EQUAL(r ,  static_cast<cast_type>(n1));
+   BOOST_CHECK_EQUAL(r, static_cast<cast_type>(n1));
    r = static_cast<simple_cast_type>(n2);
-   BOOST_CHECK_EQUAL(r ,  static_cast<cast_type>(n2));
+   BOOST_CHECK_EQUAL(r, static_cast<cast_type>(n2));
    r = static_cast<simple_cast_type>(n3);
-   BOOST_CHECK_EQUAL(r ,  static_cast<cast_type>(n3));
+   BOOST_CHECK_EQUAL(r, static_cast<cast_type>(n3));
    r = static_cast<simple_cast_type>(n4);
-   BOOST_CHECK_EQUAL(r ,  static_cast<cast_type>(n4));
+   BOOST_CHECK_EQUAL(r, static_cast<cast_type>(n4));
    // Addition:
    r = static_cast<simple_cast_type>(n2);
-   BOOST_CHECK_EQUAL(r + static_cast<simple_cast_type>(n4) ,  static_cast<cast_type>(n2 + n4));
-   BOOST_CHECK_EQUAL(Real(r + static_cast<simple_cast_type>(n4)) ,  static_cast<cast_type>(n2 + n4));
+   BOOST_CHECK_EQUAL(r + static_cast<simple_cast_type>(n4), static_cast<cast_type>(n2 + n4));
+   BOOST_CHECK_EQUAL(Real(r + static_cast<simple_cast_type>(n4)), static_cast<cast_type>(n2 + n4));
    r += static_cast<simple_cast_type>(n4);
-   BOOST_CHECK_EQUAL(r ,  static_cast<cast_type>(n2 + n4));
+   BOOST_CHECK_EQUAL(r, static_cast<cast_type>(n2 + n4));
    // subtraction:
    r = static_cast<simple_cast_type>(n4);
-   BOOST_CHECK_EQUAL(r - static_cast<simple_cast_type>(n5) ,  static_cast<cast_type>(n4 - n5));
-   BOOST_CHECK_EQUAL(Real(r - static_cast<simple_cast_type>(n5)) ,  static_cast<cast_type>(n4 - n5));
+   BOOST_CHECK_EQUAL(r - static_cast<simple_cast_type>(n5), static_cast<cast_type>(n4 - n5));
+   BOOST_CHECK_EQUAL(Real(r - static_cast<simple_cast_type>(n5)), static_cast<cast_type>(n4 - n5));
    r -= static_cast<simple_cast_type>(n5);
-   BOOST_CHECK_EQUAL(r ,  static_cast<cast_type>(n4 - n5));
+   BOOST_CHECK_EQUAL(r, static_cast<cast_type>(n4 - n5));
    // Multiplication:
    r = static_cast<simple_cast_type>(n2);
-   BOOST_CHECK_EQUAL(r * static_cast<simple_cast_type>(n4) ,  static_cast<cast_type>(n2 * n4));
-   BOOST_CHECK_EQUAL(Real(r * static_cast<simple_cast_type>(n4)) ,  static_cast<cast_type>(n2 * n4));
+   BOOST_CHECK_EQUAL(r * static_cast<simple_cast_type>(n4), static_cast<cast_type>(n2 * n4));
+   BOOST_CHECK_EQUAL(Real(r * static_cast<simple_cast_type>(n4)), static_cast<cast_type>(n2 * n4));
    r *= static_cast<simple_cast_type>(n4);
-   BOOST_CHECK_EQUAL(r ,  static_cast<cast_type>(n2 * n4));
+   BOOST_CHECK_EQUAL(r, static_cast<cast_type>(n2 * n4));
    // Division:
    r = static_cast<simple_cast_type>(n1);
-   BOOST_CHECK_EQUAL(r / static_cast<simple_cast_type>(n5) ,  static_cast<cast_type>(n1 / n5));
-   BOOST_CHECK_EQUAL(Real(r / static_cast<simple_cast_type>(n5)) ,  static_cast<cast_type>(n1 / n5));
+   BOOST_CHECK_EQUAL(r / static_cast<simple_cast_type>(n5), static_cast<cast_type>(n1 / n5));
+   BOOST_CHECK_EQUAL(Real(r / static_cast<simple_cast_type>(n5)), static_cast<cast_type>(n1 / n5));
    r /= static_cast<simple_cast_type>(n5);
-   BOOST_CHECK_EQUAL(r ,  static_cast<cast_type>(n1 / n5));
+   BOOST_CHECK_EQUAL(r, static_cast<cast_type>(n1 / n5));
    //
    // Extra cases for full coverage:
    //
    r = Real(n4) + static_cast<simple_cast_type>(n5);
-   BOOST_CHECK_EQUAL(r ,  static_cast<cast_type>(n4 + n5));
+   BOOST_CHECK_EQUAL(r, static_cast<cast_type>(n4 + n5));
    r = static_cast<simple_cast_type>(n4) + Real(n5);
-   BOOST_CHECK_EQUAL(r ,  static_cast<cast_type>(n4 + n5));
+   BOOST_CHECK_EQUAL(r, static_cast<cast_type>(n4 + n5));
    r = Real(n4) - static_cast<simple_cast_type>(n5);
-   BOOST_CHECK_EQUAL(r ,  static_cast<cast_type>(n4 - n5));
+   BOOST_CHECK_EQUAL(r, static_cast<cast_type>(n4 - n5));
    r = static_cast<simple_cast_type>(n4) - Real(n5);
-   BOOST_CHECK_EQUAL(r ,  static_cast<cast_type>(n4 - n5));
+   BOOST_CHECK_EQUAL(r, static_cast<cast_type>(n4 - n5));
    r = static_cast<simple_cast_type>(n4) * Real(n5);
-   BOOST_CHECK_EQUAL(r ,  static_cast<cast_type>(n4 * n5));
+   BOOST_CHECK_EQUAL(r, static_cast<cast_type>(n4 * n5));
    r = static_cast<cast_type>(Num(4) * n4) / Real(4);
-   BOOST_CHECK_EQUAL(r ,  static_cast<cast_type>(n4));
+   BOOST_CHECK_EQUAL(r, static_cast<cast_type>(n4));
 
    Real a, b, c;
    a = 20;
    b = 30;
    c = -a + b;
-   BOOST_CHECK_EQUAL(c ,  10);
+   BOOST_CHECK_EQUAL(c, 10);
    c = b + -a;
-   BOOST_CHECK_EQUAL(c ,  10);
+   BOOST_CHECK_EQUAL(c, 10);
    n4 = 30;
-   c = -a + static_cast<cast_type>(n4);
-   BOOST_CHECK_EQUAL(c ,  10);
+   c  = -a + static_cast<cast_type>(n4);
+   BOOST_CHECK_EQUAL(c, 10);
    c = static_cast<cast_type>(n4) + -a;
-   BOOST_CHECK_EQUAL(c ,  10);
+   BOOST_CHECK_EQUAL(c, 10);
    c = -a + -b;
-   BOOST_CHECK_EQUAL(c ,  -50);
+   BOOST_CHECK_EQUAL(c, -50);
    n4 = 4;
-   c = -(a + b) + static_cast<cast_type>(n4);
-   BOOST_CHECK_EQUAL(c ,  -50+4);
+   c  = -(a + b) + static_cast<cast_type>(n4);
+   BOOST_CHECK_EQUAL(c, -50 + 4);
    n4 = 50;
+   c  = (a + b) - static_cast<cast_type>(n4);
+   BOOST_CHECK_EQUAL(c, 0);
    c = (a + b) - static_cast<cast_type>(n4);
-   BOOST_CHECK_EQUAL(c ,  0);
-   c = (a + b) - static_cast<cast_type>(n4);
-   BOOST_CHECK_EQUAL(c ,  0);
+   BOOST_CHECK_EQUAL(c, 0);
    c = a - -(b + static_cast<cast_type>(n4));
-   BOOST_CHECK_EQUAL(c ,  20 - -(30 + 50));
+   BOOST_CHECK_EQUAL(c, 20 - -(30 + 50));
    c = -(b + static_cast<cast_type>(n4)) - a;
-   BOOST_CHECK_EQUAL(c ,  -(30 + 50) - 20);
+   BOOST_CHECK_EQUAL(c, -(30 + 50) - 20);
    c = a - -b;
-   BOOST_CHECK_EQUAL(c ,  50);
+   BOOST_CHECK_EQUAL(c, 50);
    c = -a - b;
-   BOOST_CHECK_EQUAL(c ,  -50);
+   BOOST_CHECK_EQUAL(c, -50);
    c = -a - static_cast<cast_type>(n4);
-   BOOST_CHECK_EQUAL(c ,  -20 - 50);
+   BOOST_CHECK_EQUAL(c, -20 - 50);
    c = static_cast<cast_type>(n4) - -a;
-   BOOST_CHECK_EQUAL(c ,  50 + 20);
+   BOOST_CHECK_EQUAL(c, 50 + 20);
    c = -(a + b) - Real(n4);
-   BOOST_CHECK_EQUAL(c ,  -(20 + 30) - 50);
+   BOOST_CHECK_EQUAL(c, -(20 + 30) - 50);
    c = static_cast<cast_type>(n4) - (a + b);
-   BOOST_CHECK_EQUAL(c ,  0);
+   BOOST_CHECK_EQUAL(c, 0);
    c = (a + b) * static_cast<cast_type>(n4);
-   BOOST_CHECK_EQUAL(c ,  50 * 50);
+   BOOST_CHECK_EQUAL(c, 50 * 50);
    c = static_cast<cast_type>(n4) * (a + b);
-   BOOST_CHECK_EQUAL(c ,  50 * 50);
+   BOOST_CHECK_EQUAL(c, 50 * 50);
    c = a * -(b + static_cast<cast_type>(n4));
-   BOOST_CHECK_EQUAL(c ,  20 * -(30 + 50));
+   BOOST_CHECK_EQUAL(c, 20 * -(30 + 50));
    c = -(b + static_cast<cast_type>(n4)) * a;
-   BOOST_CHECK_EQUAL(c ,  20 * -(30 + 50));
+   BOOST_CHECK_EQUAL(c, 20 * -(30 + 50));
    c = a * -b;
-   BOOST_CHECK_EQUAL(c ,  20 * -30);
+   BOOST_CHECK_EQUAL(c, 20 * -30);
    c = -a * b;
-   BOOST_CHECK_EQUAL(c ,  20 * -30);
+   BOOST_CHECK_EQUAL(c, 20 * -30);
    c = -a * static_cast<cast_type>(n4);
-   BOOST_CHECK_EQUAL(c ,  -20 * 50);
+   BOOST_CHECK_EQUAL(c, -20 * 50);
    c = static_cast<cast_type>(n4) * -a;
-   BOOST_CHECK_EQUAL(c ,  -20 * 50);
+   BOOST_CHECK_EQUAL(c, -20 * 50);
    c = -(a + b) + a;
    BOOST_CHECK(-50 + 20);
    c = static_cast<cast_type>(n4) - (a + b);
-   BOOST_CHECK_EQUAL(c ,  0);
+   BOOST_CHECK_EQUAL(c, 0);
    Real d = 10;
-   c = (a + b) / d;
-   BOOST_CHECK_EQUAL(c ,  5);
+   c      = (a + b) / d;
+   BOOST_CHECK_EQUAL(c, 5);
    c = (a + b) / (d + 0);
-   BOOST_CHECK_EQUAL(c ,  5);
+   BOOST_CHECK_EQUAL(c, 5);
    c = (a + b) / static_cast<cast_type>(n4);
-   BOOST_CHECK_EQUAL(c ,  1);
+   BOOST_CHECK_EQUAL(c, 1);
    c = static_cast<cast_type>(n4) / (a + b);
-   BOOST_CHECK_EQUAL(c ,  1);
+   BOOST_CHECK_EQUAL(c, 1);
    d = 50;
    c = d / -(a + b);
-   BOOST_CHECK_EQUAL(c ,  -1);
+   BOOST_CHECK_EQUAL(c, -1);
    c = -(a + b) / d;
-   BOOST_CHECK_EQUAL(c ,  -1);
+   BOOST_CHECK_EQUAL(c, -1);
    d = 2;
    c = a / -d;
-   BOOST_CHECK_EQUAL(c ,  20 / -2);
+   BOOST_CHECK_EQUAL(c, 20 / -2);
    c = -a / d;
-   BOOST_CHECK_EQUAL(c ,  20 / -2);
+   BOOST_CHECK_EQUAL(c, 20 / -2);
    d = 50;
    c = -d / static_cast<cast_type>(n4);
-   BOOST_CHECK_EQUAL(c ,  -1);
+   BOOST_CHECK_EQUAL(c, -1);
    c = static_cast<cast_type>(n4) / -d;
-   BOOST_CHECK_EQUAL(c ,  -1);
+   BOOST_CHECK_EQUAL(c, -1);
    c = static_cast<cast_type>(n4) + a;
-   BOOST_CHECK_EQUAL(c ,  70);
+   BOOST_CHECK_EQUAL(c, 70);
    c = static_cast<cast_type>(n4) - a;
-   BOOST_CHECK_EQUAL(c ,  30);
+   BOOST_CHECK_EQUAL(c, 30);
    c = static_cast<cast_type>(n4) * a;
-   BOOST_CHECK_EQUAL(c ,  50 * 20);
+   BOOST_CHECK_EQUAL(c, 50 * 20);
 
    n1 = -2;
    n2 = -3;
    n3 = -4;
-   a = static_cast<cast_type>(n1);
-   b = static_cast<cast_type>(n2);
-   c = static_cast<cast_type>(n3);
-   d = a + b * c;
-   BOOST_CHECK_EQUAL(d ,  -2 + -3 * -4);
+   a  = static_cast<cast_type>(n1);
+   b  = static_cast<cast_type>(n2);
+   c  = static_cast<cast_type>(n3);
+   d  = a + b * c;
+   BOOST_CHECK_EQUAL(d, -2 + -3 * -4);
    d = static_cast<cast_type>(n1) + b * c;
-   BOOST_CHECK_EQUAL(d ,  -2 + -3 * -4);
+   BOOST_CHECK_EQUAL(d, -2 + -3 * -4);
    d = a + static_cast<cast_type>(n2) * c;
-   BOOST_CHECK_EQUAL(d ,  -2 + -3 * -4);
+   BOOST_CHECK_EQUAL(d, -2 + -3 * -4);
    d = a + b * static_cast<cast_type>(n3);
-   BOOST_CHECK_EQUAL(d ,  -2 + -3 * -4);
+   BOOST_CHECK_EQUAL(d, -2 + -3 * -4);
    d = static_cast<cast_type>(n1) + static_cast<cast_type>(n2) * c;
-   BOOST_CHECK_EQUAL(d ,  -2 + -3 * -4);
+   BOOST_CHECK_EQUAL(d, -2 + -3 * -4);
    d = static_cast<cast_type>(n1) + b * static_cast<cast_type>(n3);
-   BOOST_CHECK_EQUAL(d ,  -2 + -3 * -4);
+   BOOST_CHECK_EQUAL(d, -2 + -3 * -4);
    a += static_cast<cast_type>(n2) * c;
-   BOOST_CHECK_EQUAL(a ,  -2 + -3 * -4);
+   BOOST_CHECK_EQUAL(a, -2 + -3 * -4);
    a = static_cast<cast_type>(n1);
    a += b * static_cast<cast_type>(n3);
-   BOOST_CHECK_EQUAL(a ,  -2 + -3 * -4);
+   BOOST_CHECK_EQUAL(a, -2 + -3 * -4);
    a = static_cast<cast_type>(n1);
 
    d = b * c + a;
-   BOOST_CHECK_EQUAL(d ,  -2 + -3 * -4);
+   BOOST_CHECK_EQUAL(d, -2 + -3 * -4);
    d = b * c + static_cast<cast_type>(n1);
-   BOOST_CHECK_EQUAL(d ,  -2 + -3 * -4);
+   BOOST_CHECK_EQUAL(d, -2 + -3 * -4);
    d = static_cast<cast_type>(n2) * c + a;
-   BOOST_CHECK_EQUAL(d ,  -2 + -3 * -4);
+   BOOST_CHECK_EQUAL(d, -2 + -3 * -4);
    d = b * static_cast<cast_type>(n3) + a;
-   BOOST_CHECK_EQUAL(d ,  -2 + -3 * -4);
+   BOOST_CHECK_EQUAL(d, -2 + -3 * -4);
    d = static_cast<cast_type>(n2) * c + static_cast<cast_type>(n1);
-   BOOST_CHECK_EQUAL(d ,  -2 + -3 * -4);
+   BOOST_CHECK_EQUAL(d, -2 + -3 * -4);
    d = b * static_cast<cast_type>(n3) + static_cast<cast_type>(n1);
-   BOOST_CHECK_EQUAL(d ,  -2 + -3 * -4);
+   BOOST_CHECK_EQUAL(d, -2 + -3 * -4);
 
    a = -20;
    d = a - b * c;
-   BOOST_CHECK_EQUAL(d ,  -20 - -3 * -4);
+   BOOST_CHECK_EQUAL(d, -20 - -3 * -4);
    n1 = -20;
-   d = static_cast<cast_type>(n1) - b * c;
-   BOOST_CHECK_EQUAL(d ,  -20 - -3 * -4);
+   d  = static_cast<cast_type>(n1) - b * c;
+   BOOST_CHECK_EQUAL(d, -20 - -3 * -4);
    d = a - static_cast<cast_type>(n2) * c;
-   BOOST_CHECK_EQUAL(d ,  -20 - -3 * -4);
+   BOOST_CHECK_EQUAL(d, -20 - -3 * -4);
    d = a - b * static_cast<cast_type>(n3);
-   BOOST_CHECK_EQUAL(d ,  -20 - -3 * -4);
+   BOOST_CHECK_EQUAL(d, -20 - -3 * -4);
    d = static_cast<cast_type>(n1) - static_cast<cast_type>(n2) * c;
-   BOOST_CHECK_EQUAL(d ,  -20 - -3 * -4);
+   BOOST_CHECK_EQUAL(d, -20 - -3 * -4);
    d = static_cast<cast_type>(n1) - b * static_cast<cast_type>(n3);
-   BOOST_CHECK_EQUAL(d ,  -20 - -3 * -4);
+   BOOST_CHECK_EQUAL(d, -20 - -3 * -4);
    a -= static_cast<cast_type>(n2) * c;
-   BOOST_CHECK_EQUAL(a ,  -20 - -3 * -4);
+   BOOST_CHECK_EQUAL(a, -20 - -3 * -4);
    a = static_cast<cast_type>(n1);
    a -= b * static_cast<cast_type>(n3);
-   BOOST_CHECK_EQUAL(a ,  -20 - -3 * -4);
+   BOOST_CHECK_EQUAL(a, -20 - -3 * -4);
 
    a = -2;
    d = b * c - a;
-   BOOST_CHECK_EQUAL(d ,  -3 * -4 - -2);
+   BOOST_CHECK_EQUAL(d, -3 * -4 - -2);
    n1 = -2;
-   d = b * c - static_cast<cast_type>(n1);
-   BOOST_CHECK_EQUAL(d ,  -3 * -4 - -2);
+   d  = b * c - static_cast<cast_type>(n1);
+   BOOST_CHECK_EQUAL(d, -3 * -4 - -2);
    d = static_cast<cast_type>(n2) * c - a;
-   BOOST_CHECK_EQUAL(d ,  -3 * -4 - -2);
+   BOOST_CHECK_EQUAL(d, -3 * -4 - -2);
    d = b * static_cast<cast_type>(n3) - a;
-   BOOST_CHECK_EQUAL(d ,  -3 * -4 - -2);
+   BOOST_CHECK_EQUAL(d, -3 * -4 - -2);
    d = static_cast<cast_type>(n2) * c - static_cast<cast_type>(n1);
-   BOOST_CHECK_EQUAL(d ,  -3 * -4 - -2);
+   BOOST_CHECK_EQUAL(d, -3 * -4 - -2);
    d = b * static_cast<cast_type>(n3) - static_cast<cast_type>(n1);
-   BOOST_CHECK_EQUAL(d ,  -3 * -4 - -2);
+   BOOST_CHECK_EQUAL(d, -3 * -4 - -2);
    //
    // Conversion from min and max values:
    //
-   test_negative_mixed_minmax<Real, Num>(boost::mpl::bool_<std::numeric_limits<Real>::is_integer && std::numeric_limits<Num>::is_integer>());
+   test_negative_mixed_minmax<Real, Num>(boost::mpl::bool_ < std::numeric_limits<Real>::is_integer && std::numeric_limits<Num>::is_integer > ());
 }
 
 template <class Real, class Num>
@@ -1763,17 +1762,18 @@ void test_mixed_numeric_limits(const boost::mpl::true_&)
    }
 
    static const int left_shift = std::numeric_limits<Num>::digits - 1;
-   Num n1 = static_cast<Num>(1uLL << ((left_shift < 63) && (left_shift > 0) ? left_shift : 10));
-   Num n2 = 1;
-   Num n3 = 0;
-   Num n4 = 20;
-
-   std::ios_base::fmtflags f = boost::is_floating_point<Num>::value ? std::ios_base::scientific : std::ios_base::fmtflags(0);
-   int digits_to_print = boost::is_floating_point<Num>::value && std::numeric_limits<Num>::is_specialized
-      ? std::numeric_limits<Num>::digits10 + 5 : 0;
-   if(std::numeric_limits<target_type>::digits <= std::numeric_limits<Real>::digits)
+   Num              n1         = static_cast<Num>(1uLL << ((left_shift < 63) && (left_shift > 0) ? left_shift : 10));
+   Num              n2         = 1;
+   Num              n3         = 0;
+   Num              n4         = 20;
+
+   std::ios_base::fmtflags f               = boost::is_floating_point<Num>::value ? std::ios_base::scientific : std::ios_base::fmtflags(0);
+   int                     digits_to_print = boost::is_floating_point<Num>::value && std::numeric_limits<Num>::is_specialized
+                             ? std::numeric_limits<Num>::digits10 + 5
+                             : 0;
+   if (std::numeric_limits<target_type>::digits <= std::numeric_limits<Real>::digits)
    {
-   BOOST_CHECK_CLOSE(n1, checked_lexical_cast<target_type>(Real(n1).str(digits_to_print, f)), tol);
+      BOOST_CHECK_CLOSE(n1, checked_lexical_cast<target_type>(Real(n1).str(digits_to_print, f)), tol);
    }
    BOOST_CHECK_CLOSE(n2, checked_lexical_cast<target_type>(Real(n2).str(digits_to_print, f)), 0);
    BOOST_CHECK_CLOSE(n3, checked_lexical_cast<target_type>(Real(n3).str(digits_to_print, f)), 0);
@@ -1788,26 +1788,24 @@ template <class Real, class Num>
 void test_mixed(const boost::mpl::true_&)
 {
    typedef typename boost::mpl::if_<
-         boost::is_convertible<Num, Real>,
-         typename boost::mpl::if_c<boost::is_integral<Num>::value && (sizeof(Num) < sizeof(int)), int, Num>::type,
-         Real
-      >::type cast_type;
+       boost::is_convertible<Num, Real>,
+       typename boost::mpl::if_c<boost::is_integral<Num>::value && (sizeof(Num) < sizeof(int)), int, Num>::type,
+       Real>::type cast_type;
    typedef typename boost::mpl::if_<
-         boost::is_convertible<Num, Real>,
-         Num,
-         Real
-      >::type simple_cast_type;
+       boost::is_convertible<Num, Real>,
+       Num,
+       Real>::type simple_cast_type;
 
-   if(std::numeric_limits<Real>::is_specialized && std::numeric_limits<Real>::is_bounded && std::numeric_limits<Real>::digits < std::numeric_limits<Num>::digits)
+   if (std::numeric_limits<Real>::is_specialized && std::numeric_limits<Real>::is_bounded && std::numeric_limits<Real>::digits < std::numeric_limits<Num>::digits)
       return;
 
-   std::cout << "Testing mixed arithmetic with type: " << typeid(Real).name()  << " and " << typeid(Num).name() << std::endl;
+   std::cout << "Testing mixed arithmetic with type: " << typeid(Real).name() << " and " << typeid(Num).name() << std::endl;
    static const int left_shift = std::numeric_limits<Num>::digits - 1;
-   Num n1 = static_cast<Num>(1uLL << ((left_shift < 63) && (left_shift > 0) ? left_shift : 10));
-   Num n2 = 1;
-   Num n3 = 0;
-   Num n4 = 20;
-   Num n5 = 8;
+   Num              n1         = static_cast<Num>(1uLL << ((left_shift < 63) && (left_shift > 0) ? left_shift : 10));
+   Num              n2         = 1;
+   Num              n3         = 0;
+   Num              n4         = 20;
+   Num              n5         = 8;
 
    test_comparisons<Real>(n1, n2, boost::is_convertible<Num, Real>());
    test_comparisons<Real>(n1, n3, boost::is_convertible<Num, Real>());
@@ -1817,19 +1815,19 @@ void test_mixed(const boost::mpl::true_&)
    test_comparisons<Real>(n3, n3, boost::is_convertible<Num, Real>());
 
    // Default construct:
-   BOOST_CHECK_EQUAL(Real(n1) ,  static_cast<cast_type>(n1));
-   BOOST_CHECK_EQUAL(Real(n2) ,  static_cast<cast_type>(n2));
-   BOOST_CHECK_EQUAL(Real(n3) ,  static_cast<cast_type>(n3));
-   BOOST_CHECK_EQUAL(Real(n4) ,  static_cast<cast_type>(n4));
-   BOOST_CHECK_EQUAL(Real(n1).template convert_to<Num>() ,  n1);
-   BOOST_CHECK_EQUAL(Real(n2).template convert_to<Num>() ,  n2);
-   BOOST_CHECK_EQUAL(Real(n3).template convert_to<Num>() ,  n3);
-   BOOST_CHECK_EQUAL(Real(n4).template convert_to<Num>() ,  n4);
+   BOOST_CHECK_EQUAL(Real(n1), static_cast<cast_type>(n1));
+   BOOST_CHECK_EQUAL(Real(n2), static_cast<cast_type>(n2));
+   BOOST_CHECK_EQUAL(Real(n3), static_cast<cast_type>(n3));
+   BOOST_CHECK_EQUAL(Real(n4), static_cast<cast_type>(n4));
+   BOOST_CHECK_EQUAL(Real(n1).template convert_to<Num>(), n1);
+   BOOST_CHECK_EQUAL(Real(n2).template convert_to<Num>(), n2);
+   BOOST_CHECK_EQUAL(Real(n3).template convert_to<Num>(), n3);
+   BOOST_CHECK_EQUAL(Real(n4).template convert_to<Num>(), n4);
 #ifndef BOOST_MP_NO_CXX11_EXPLICIT_CONVERSION_OPERATORS
-   BOOST_CHECK_EQUAL(static_cast<Num>(Real(n1)) ,  n1);
-   BOOST_CHECK_EQUAL(static_cast<Num>(Real(n2)) ,  n2);
-   BOOST_CHECK_EQUAL(static_cast<Num>(Real(n3)) ,  n3);
-   BOOST_CHECK_EQUAL(static_cast<Num>(Real(n4)) ,  n4);
+   BOOST_CHECK_EQUAL(static_cast<Num>(Real(n1)), n1);
+   BOOST_CHECK_EQUAL(static_cast<Num>(Real(n2)), n2);
+   BOOST_CHECK_EQUAL(static_cast<Num>(Real(n3)), n3);
+   BOOST_CHECK_EQUAL(static_cast<Num>(Real(n4)), n4);
 #endif
    // Again with expression templates:
    BOOST_CHECK_EQUAL((Real(n1) + 0).template convert_to<Num>(), n1);
@@ -1843,59 +1841,59 @@ void test_mixed(const boost::mpl::true_&)
    BOOST_CHECK_EQUAL(static_cast<Num>(Real(n4) + 0), n4);
 #endif
    BOOST_CHECK_EQUAL(static_cast<cast_type>(n1), Real(n1));
-   BOOST_CHECK_EQUAL(static_cast<cast_type>(n2) ,  Real(n2));
-   BOOST_CHECK_EQUAL(static_cast<cast_type>(n3) ,  Real(n3));
-   BOOST_CHECK_EQUAL(static_cast<cast_type>(n4) ,  Real(n4));
+   BOOST_CHECK_EQUAL(static_cast<cast_type>(n2), Real(n2));
+   BOOST_CHECK_EQUAL(static_cast<cast_type>(n3), Real(n3));
+   BOOST_CHECK_EQUAL(static_cast<cast_type>(n4), Real(n4));
    // Assignment:
    Real r(0);
    BOOST_CHECK(r != static_cast<cast_type>(n1));
    r = static_cast<simple_cast_type>(n1);
-   BOOST_CHECK_EQUAL(r ,  static_cast<cast_type>(n1));
+   BOOST_CHECK_EQUAL(r, static_cast<cast_type>(n1));
    r = static_cast<simple_cast_type>(n2);
-   BOOST_CHECK_EQUAL(r ,  static_cast<cast_type>(n2));
+   BOOST_CHECK_EQUAL(r, static_cast<cast_type>(n2));
    r = static_cast<simple_cast_type>(n3);
-   BOOST_CHECK_EQUAL(r ,  static_cast<cast_type>(n3));
+   BOOST_CHECK_EQUAL(r, static_cast<cast_type>(n3));
    r = static_cast<simple_cast_type>(n4);
-   BOOST_CHECK_EQUAL(r ,  static_cast<cast_type>(n4));
+   BOOST_CHECK_EQUAL(r, static_cast<cast_type>(n4));
    // Addition:
    r = static_cast<simple_cast_type>(n2);
-   BOOST_CHECK_EQUAL(r + static_cast<simple_cast_type>(n4) ,  static_cast<cast_type>(n2 + n4));
-   BOOST_CHECK_EQUAL(Real(r + static_cast<simple_cast_type>(n4)) ,  static_cast<cast_type>(n2 + n4));
+   BOOST_CHECK_EQUAL(r + static_cast<simple_cast_type>(n4), static_cast<cast_type>(n2 + n4));
+   BOOST_CHECK_EQUAL(Real(r + static_cast<simple_cast_type>(n4)), static_cast<cast_type>(n2 + n4));
    r += static_cast<simple_cast_type>(n4);
-   BOOST_CHECK_EQUAL(r ,  static_cast<cast_type>(n2 + n4));
+   BOOST_CHECK_EQUAL(r, static_cast<cast_type>(n2 + n4));
    // subtraction:
    r = static_cast<simple_cast_type>(n4);
-   BOOST_CHECK_EQUAL(r - static_cast<simple_cast_type>(n5) ,  static_cast<cast_type>(n4 - n5));
-   BOOST_CHECK_EQUAL(Real(r - static_cast<simple_cast_type>(n5)) ,  static_cast<cast_type>(n4 - n5));
+   BOOST_CHECK_EQUAL(r - static_cast<simple_cast_type>(n5), static_cast<cast_type>(n4 - n5));
+   BOOST_CHECK_EQUAL(Real(r - static_cast<simple_cast_type>(n5)), static_cast<cast_type>(n4 - n5));
    r -= static_cast<simple_cast_type>(n5);
-   BOOST_CHECK_EQUAL(r ,  static_cast<cast_type>(n4 - n5));
+   BOOST_CHECK_EQUAL(r, static_cast<cast_type>(n4 - n5));
    // Multiplication:
    r = static_cast<simple_cast_type>(n2);
-   BOOST_CHECK_EQUAL(r * static_cast<simple_cast_type>(n4) ,  static_cast<cast_type>(n2 * n4));
-   BOOST_CHECK_EQUAL(Real(r * static_cast<simple_cast_type>(n4)) ,  static_cast<cast_type>(n2 * n4));
+   BOOST_CHECK_EQUAL(r * static_cast<simple_cast_type>(n4), static_cast<cast_type>(n2 * n4));
+   BOOST_CHECK_EQUAL(Real(r * static_cast<simple_cast_type>(n4)), static_cast<cast_type>(n2 * n4));
    r *= static_cast<simple_cast_type>(n4);
-   BOOST_CHECK_EQUAL(r ,  static_cast<cast_type>(n2 * n4));
+   BOOST_CHECK_EQUAL(r, static_cast<cast_type>(n2 * n4));
    // Division:
    r = static_cast<simple_cast_type>(n1);
-   BOOST_CHECK_EQUAL(r / static_cast<simple_cast_type>(n5) ,  static_cast<cast_type>(n1 / n5));
-   BOOST_CHECK_EQUAL(Real(r / static_cast<simple_cast_type>(n5)) ,  static_cast<cast_type>(n1 / n5));
+   BOOST_CHECK_EQUAL(r / static_cast<simple_cast_type>(n5), static_cast<cast_type>(n1 / n5));
+   BOOST_CHECK_EQUAL(Real(r / static_cast<simple_cast_type>(n5)), static_cast<cast_type>(n1 / n5));
    r /= static_cast<simple_cast_type>(n5);
-   BOOST_CHECK_EQUAL(r ,  static_cast<cast_type>(n1 / n5));
+   BOOST_CHECK_EQUAL(r, static_cast<cast_type>(n1 / n5));
    //
    // special cases for full coverage:
    //
    r = static_cast<simple_cast_type>(n5) + Real(n4);
-   BOOST_CHECK_EQUAL(r ,  static_cast<cast_type>(n4 + n5));
+   BOOST_CHECK_EQUAL(r, static_cast<cast_type>(n4 + n5));
    r = static_cast<simple_cast_type>(n4) - Real(n5);
-   BOOST_CHECK_EQUAL(r ,  static_cast<cast_type>(n4 - n5));
+   BOOST_CHECK_EQUAL(r, static_cast<cast_type>(n4 - n5));
    r = static_cast<simple_cast_type>(n4) * Real(n5);
-   BOOST_CHECK_EQUAL(r ,  static_cast<cast_type>(n4 * n5));
+   BOOST_CHECK_EQUAL(r, static_cast<cast_type>(n4 * n5));
    r = static_cast<cast_type>(Num(4) * n4) / Real(4);
-   BOOST_CHECK_EQUAL(r ,  static_cast<cast_type>(n4));
+   BOOST_CHECK_EQUAL(r, static_cast<cast_type>(n4));
 
    typedef boost::mpl::bool_<
-      (!std::numeric_limits<Num>::is_specialized || std::numeric_limits<Num>::is_signed)
-      && (!std::numeric_limits<Real>::is_specialized || std::numeric_limits<Real>::is_signed)> signed_tag;
+       (!std::numeric_limits<Num>::is_specialized || std::numeric_limits<Num>::is_signed) && (!std::numeric_limits<Real>::is_specialized || std::numeric_limits<Real>::is_signed)>
+       signed_tag;
 
    test_negative_mixed<Real, Num>(signed_tag());
 
@@ -1904,71 +1902,71 @@ void test_mixed(const boost::mpl::true_&)
    n3 = 4;
    Real a(n1), b(n2), c(n3), d;
    d = a + b * c;
-   BOOST_CHECK_EQUAL(d ,  2 + 3 * 4);
+   BOOST_CHECK_EQUAL(d, 2 + 3 * 4);
    d = static_cast<cast_type>(n1) + b * c;
-   BOOST_CHECK_EQUAL(d ,  2 + 3 * 4);
+   BOOST_CHECK_EQUAL(d, 2 + 3 * 4);
    d = a + static_cast<cast_type>(n2) * c;
-   BOOST_CHECK_EQUAL(d ,  2 + 3 * 4);
+   BOOST_CHECK_EQUAL(d, 2 + 3 * 4);
    d = a + b * static_cast<cast_type>(n3);
-   BOOST_CHECK_EQUAL(d ,  2 + 3 * 4);
+   BOOST_CHECK_EQUAL(d, 2 + 3 * 4);
    d = static_cast<cast_type>(n1) + static_cast<cast_type>(n2) * c;
-   BOOST_CHECK_EQUAL(d ,  2 + 3 * 4);
+   BOOST_CHECK_EQUAL(d, 2 + 3 * 4);
    d = static_cast<cast_type>(n1) + b * static_cast<cast_type>(n3);
-   BOOST_CHECK_EQUAL(d ,  2 + 3 * 4);
+   BOOST_CHECK_EQUAL(d, 2 + 3 * 4);
    a += static_cast<cast_type>(n2) * c;
-   BOOST_CHECK_EQUAL(a ,  2 + 3 * 4);
+   BOOST_CHECK_EQUAL(a, 2 + 3 * 4);
    a = static_cast<cast_type>(n1);
    a += b * static_cast<cast_type>(n3);
-   BOOST_CHECK_EQUAL(a ,  2 + 3 * 4);
+   BOOST_CHECK_EQUAL(a, 2 + 3 * 4);
    a = static_cast<cast_type>(n1);
 
    d = b * c + a;
-   BOOST_CHECK_EQUAL(d ,  2 + 3 * 4);
+   BOOST_CHECK_EQUAL(d, 2 + 3 * 4);
    d = b * c + static_cast<cast_type>(n1);
-   BOOST_CHECK_EQUAL(d ,  2 + 3 * 4);
+   BOOST_CHECK_EQUAL(d, 2 + 3 * 4);
    d = static_cast<cast_type>(n2) * c + a;
-   BOOST_CHECK_EQUAL(d ,  2 + 3 * 4);
+   BOOST_CHECK_EQUAL(d, 2 + 3 * 4);
    d = b * static_cast<cast_type>(n3) + a;
-   BOOST_CHECK_EQUAL(d ,  2 + 3 * 4);
+   BOOST_CHECK_EQUAL(d, 2 + 3 * 4);
    d = static_cast<cast_type>(n2) * c + static_cast<cast_type>(n1);
-   BOOST_CHECK_EQUAL(d ,  2 + 3 * 4);
+   BOOST_CHECK_EQUAL(d, 2 + 3 * 4);
    d = b * static_cast<cast_type>(n3) + static_cast<cast_type>(n1);
-   BOOST_CHECK_EQUAL(d ,  2 + 3 * 4);
+   BOOST_CHECK_EQUAL(d, 2 + 3 * 4);
 
    a = 20;
    d = a - b * c;
-   BOOST_CHECK_EQUAL(d ,  20 - 3 * 4);
+   BOOST_CHECK_EQUAL(d, 20 - 3 * 4);
    n1 = 20;
-   d = static_cast<cast_type>(n1) - b * c;
-   BOOST_CHECK_EQUAL(d ,  20 - 3 * 4);
+   d  = static_cast<cast_type>(n1) - b * c;
+   BOOST_CHECK_EQUAL(d, 20 - 3 * 4);
    d = a - static_cast<cast_type>(n2) * c;
-   BOOST_CHECK_EQUAL(d ,  20 - 3 * 4);
+   BOOST_CHECK_EQUAL(d, 20 - 3 * 4);
    d = a - b * static_cast<cast_type>(n3);
-   BOOST_CHECK_EQUAL(d ,  20 - 3 * 4);
+   BOOST_CHECK_EQUAL(d, 20 - 3 * 4);
    d = static_cast<cast_type>(n1) - static_cast<cast_type>(n2) * c;
-   BOOST_CHECK_EQUAL(d ,  20 - 3 * 4);
+   BOOST_CHECK_EQUAL(d, 20 - 3 * 4);
    d = static_cast<cast_type>(n1) - b * static_cast<cast_type>(n3);
-   BOOST_CHECK_EQUAL(d ,  20 - 3 * 4);
+   BOOST_CHECK_EQUAL(d, 20 - 3 * 4);
    a -= static_cast<cast_type>(n2) * c;
-   BOOST_CHECK_EQUAL(a ,  20 - 3 * 4);
+   BOOST_CHECK_EQUAL(a, 20 - 3 * 4);
    a = static_cast<cast_type>(n1);
    a -= b * static_cast<cast_type>(n3);
-   BOOST_CHECK_EQUAL(a ,  20 - 3 * 4);
+   BOOST_CHECK_EQUAL(a, 20 - 3 * 4);
 
    a = 2;
    d = b * c - a;
-   BOOST_CHECK_EQUAL(d ,  3 * 4 - 2);
+   BOOST_CHECK_EQUAL(d, 3 * 4 - 2);
    n1 = 2;
-   d = b * c - static_cast<cast_type>(n1);
-   BOOST_CHECK_EQUAL(d ,  3 * 4 - 2);
+   d  = b * c - static_cast<cast_type>(n1);
+   BOOST_CHECK_EQUAL(d, 3 * 4 - 2);
    d = static_cast<cast_type>(n2) * c - a;
-   BOOST_CHECK_EQUAL(d ,  3 * 4 - 2);
+   BOOST_CHECK_EQUAL(d, 3 * 4 - 2);
    d = b * static_cast<cast_type>(n3) - a;
-   BOOST_CHECK_EQUAL(d ,  3 * 4 - a);
+   BOOST_CHECK_EQUAL(d, 3 * 4 - a);
    d = static_cast<cast_type>(n2) * c - static_cast<cast_type>(n1);
-   BOOST_CHECK_EQUAL(d ,  3 * 4 - 2);
+   BOOST_CHECK_EQUAL(d, 3 * 4 - 2);
    d = b * static_cast<cast_type>(n3) - static_cast<cast_type>(n1);
-   BOOST_CHECK_EQUAL(d ,  3 * 4 - 2);
+   BOOST_CHECK_EQUAL(d, 3 * 4 - 2);
 
    test_mixed_numeric_limits<Real, Num>(boost::mpl::bool_<std::numeric_limits<Real>::is_specialized>());
 }
@@ -2020,43 +2018,43 @@ typename boost::enable_if_c<boost::multiprecision::number_category<Real>::value
    }
 #ifndef BOOST_NO_CXX17_HDR_STRING_VIEW
    {
-      std::string sx("40.550"), sy("222");
+      std::string      sx("40.550"), sy("222");
       std::string_view vx(sx.c_str(), 4), vy(sy.c_str(), 1);
-      Real d(vx, vy);
+      Real             d(vx, vy);
       BOOST_CHECK_EQUAL(d.real(), 40.5);
       BOOST_CHECK_EQUAL(d.imag(), 2);
    }
 #endif
    {
       typename Real::value_type x(40.5), y(2);
-      Real d(x, y);
+      Real                      d(x, y);
       BOOST_CHECK_EQUAL(d.real(), 40.5);
       BOOST_CHECK_EQUAL(d.imag(), 2);
    }
 #ifdef TEST_MPC
    {
       typename Real::value_type x(40.5), y(2);
-      Real d(x.backend().data(), y.backend().data());
+      Real                      d(x.backend().data(), y.backend().data());
       BOOST_CHECK_EQUAL(d.real(), 40.5);
       BOOST_CHECK_EQUAL(d.imag(), 2);
    }
 #endif
    {
       typename Real::value_type x(40.5);
-      Real d(x, 2);
+      Real                      d(x, 2);
       BOOST_CHECK_EQUAL(d.real(), 40.5);
       BOOST_CHECK_EQUAL(d.imag(), 2);
    }
    {
       typename Real::value_type x(40.5);
-      Real d(2, x);
+      Real                      d(2, x);
       BOOST_CHECK_EQUAL(d.imag(), 40.5);
       BOOST_CHECK_EQUAL(d.real(), 2);
    }
    {
-      typename Real::value_type x(real(a)*real(b) + imag(a)*imag(b)), y(imag(a)*real(b) - real(a)*imag(b));
-      Real d(real(a)*real(b) + imag(a)*imag(b), imag(a)*real(b) - real(a)*imag(b));
-      Real e(x, y);
+      typename Real::value_type x(real(a) * real(b) + imag(a) * imag(b)), y(imag(a) * real(b) - real(a) * imag(b));
+      Real                      d(real(a) * real(b) + imag(a) * imag(b), imag(a) * real(b) - real(a) * imag(b));
+      Real                      e(x, y);
       BOOST_CHECK_EQUAL(d, e);
    }
    //
@@ -2111,8 +2109,7 @@ typename boost::enable_if_c<boost::multiprecision::number_category<Real>::value
    c.imag(29ULL);
    BOOST_CHECK_EQUAL(imag(c), 29);
 #endif
-   
-   
+
    c.real(2).imag(3);
 
    BOOST_CHECK_EQUAL(real(a), 30);
@@ -2123,7 +2120,7 @@ typename boost::enable_if_c<boost::multiprecision::number_category<Real>::value
    BOOST_CHECK_EQUAL(imag(a + 0), 0);
    BOOST_CHECK_EQUAL(real(c + 0), 2);
    BOOST_CHECK_EQUAL(imag(c + 0), 3);
-   
+
    // string construction:
    a = Real("2");
    BOOST_CHECK_EQUAL(real(a), 2);
@@ -2387,21 +2384,21 @@ typename boost::enable_if_c<boost::multiprecision::number_category<Real>::value
    Real b = 30;
    BOOST_CHECK(a.sign() > 0);
    BOOST_CHECK(!a.is_zero());
-   if(std::numeric_limits<Real>::is_signed)
+   if (std::numeric_limits<Real>::is_signed)
    {
       a = -20;
       BOOST_CHECK(a.sign() < 0);
       BOOST_CHECK(!a.is_zero());
    }
    a = 0;
-   BOOST_CHECK_EQUAL(a.sign() ,  0);
+   BOOST_CHECK_EQUAL(a.sign(), 0);
    BOOST_CHECK(a.is_zero());
 
    a = 20;
    b = 30;
    a.swap(b);
-   BOOST_CHECK_EQUAL(a ,  30);
-   BOOST_CHECK_EQUAL(b ,  20);
+   BOOST_CHECK_EQUAL(a, 30);
+   BOOST_CHECK_EQUAL(b, 20);
    //
    // Test complex number functions which are also overloaded for scalar type:
    //
@@ -2434,38 +2431,38 @@ void test_signed_ops(const boost::mpl::true_&)
    Real c(500);
    Real d(1024);
    Real ac;
-   BOOST_CHECK_EQUAL(-a ,  -8);
+   BOOST_CHECK_EQUAL(-a, -8);
    ac = a;
    ac = ac - b;
-   BOOST_CHECK_EQUAL(ac ,  8 - 64);
+   BOOST_CHECK_EQUAL(ac, 8 - 64);
    ac = a;
    ac -= a + b;
-   BOOST_CHECK_EQUAL(ac ,  -64);
+   BOOST_CHECK_EQUAL(ac, -64);
    ac = a;
    ac -= b - a;
-   BOOST_CHECK_EQUAL(ac ,  16 - 64);
+   BOOST_CHECK_EQUAL(ac, 16 - 64);
    ac = -a;
-   BOOST_CHECK_EQUAL(ac ,  -8);
+   BOOST_CHECK_EQUAL(ac, -8);
    ac = a;
    ac -= -a;
-   BOOST_CHECK_EQUAL(ac ,  16);
+   BOOST_CHECK_EQUAL(ac, 16);
    ac = a;
    ac += -a;
-   BOOST_CHECK_EQUAL(ac ,  0);
+   BOOST_CHECK_EQUAL(ac, 0);
    ac = b;
    ac /= -a;
-   BOOST_CHECK_EQUAL(ac ,  -8);
+   BOOST_CHECK_EQUAL(ac, -8);
    ac = a;
    ac *= -a;
-   BOOST_CHECK_EQUAL(ac ,  -64);
+   BOOST_CHECK_EQUAL(ac, -64);
    ac = a + -b;
-   BOOST_CHECK_EQUAL(ac ,  8 - 64);
+   BOOST_CHECK_EQUAL(ac, 8 - 64);
    ac = -a + b;
-   BOOST_CHECK_EQUAL(ac ,  -8+64);
+   BOOST_CHECK_EQUAL(ac, -8 + 64);
    ac = -a + -b;
-   BOOST_CHECK_EQUAL(ac ,  -72);
-   ac = a + - + -b; // lots of unary operators!!
-   BOOST_CHECK_EQUAL(ac ,  72);
+   BOOST_CHECK_EQUAL(ac, -72);
+   ac = a + -+-b; // lots of unary operators!!
+   BOOST_CHECK_EQUAL(ac, 72);
    test_conditional(Real(-a), -a);
 }
 template <class Real>
@@ -2476,47 +2473,53 @@ void test_signed_ops(const boost::mpl::false_&)
 template <class Real>
 void test_basic_conditionals(Real a, Real b)
 {
-   if(a)
+   if (a)
    {
       BOOST_ERROR("Unexpected non-zero result");
    }
-   if(!a){}
+   if (!a)
+   {
+   }
    else
    {
       BOOST_ERROR("Unexpected zero result");
    }
    b = 2;
-   if(!b)
+   if (!b)
    {
       BOOST_ERROR("Unexpected zero result");
    }
-   if(b){}
+   if (b)
+   {
+   }
    else
    {
       BOOST_ERROR("Unexpected non-zero result");
    }
-   if(a && b)
+   if (a && b)
    {
       BOOST_ERROR("Unexpected zero result");
    }
-   if(!(a || b))
+   if (!(a || b))
    {
       BOOST_ERROR("Unexpected zero result");
    }
-   if(a + b){}
+   if (a + b)
+   {
+   }
    else
    {
       BOOST_ERROR("Unexpected zero result");
    }
-   if(b - 2)
+   if (b - 2)
    {
       BOOST_ERROR("Unexpected non-zero result");
    }
 }
 
 template <class T>
-typename boost::enable_if_c<boost::multiprecision::number_category<T>::value == boost::multiprecision::number_kind_complex>::type 
-   test_relationals(T a, T b)
+typename boost::enable_if_c<boost::multiprecision::number_category<T>::value == boost::multiprecision::number_kind_complex>::type
+test_relationals(T a, T b)
 {
    BOOST_CHECK_EQUAL((a == b), false);
    BOOST_CHECK_EQUAL((a != b), true);
@@ -2620,41 +2623,44 @@ test_relationals(T a, T b)
    // min and max overloads:
    //
 #if !defined(min) && !defined(max)
-   using std::max;
-   using std::min;
+//   using std::max;
+//   using std::min;
+// This works, but still causes complaints from inspect.exe, so use brackets to prevent macrosubstitution,
+// and to explicitly specify type T seems necessary, for reasons unclear.
    a = 2;
    b = 5;
    c = 6;
-   BOOST_CHECK_EQUAL(min(a, b), a);
-   BOOST_CHECK_EQUAL(min(b, a), a);
-   BOOST_CHECK_EQUAL(max(a, b), b);
-   BOOST_CHECK_EQUAL(max(b, a), b);
-   BOOST_CHECK_EQUAL(min(a, b + c), a);
-   BOOST_CHECK_EQUAL(min(b + c, a), a);
-   BOOST_CHECK_EQUAL(min(a, c - b), 1);
-   BOOST_CHECK_EQUAL(min(c - b, a), 1);
-   BOOST_CHECK_EQUAL(max(a, b + c), 11);
-   BOOST_CHECK_EQUAL(max(b + c, a), 11);
-   BOOST_CHECK_EQUAL(max(a, c - b), a);
-   BOOST_CHECK_EQUAL(max(c - b, a), a);
-   BOOST_CHECK_EQUAL(min(a + b, b + c), 7);
-   BOOST_CHECK_EQUAL(min(b + c, a + b), 7);
-   BOOST_CHECK_EQUAL(max(a + b, b + c), 11);
-   BOOST_CHECK_EQUAL(max(b + c, a + b), 11);
-   BOOST_CHECK_EQUAL(min(a + b, c - a), 4);
-   BOOST_CHECK_EQUAL(min(c - a, a + b), 4);
-   BOOST_CHECK_EQUAL(max(a + b, c - a), 7);
-   BOOST_CHECK_EQUAL(max(c - a, a + b), 7);
+   BOOST_CHECK_EQUAL( (std::min<T>)(a, b),  a);
+   BOOST_CHECK_EQUAL( (std::min<T>)(b, a), a);
+   BOOST_CHECK_EQUAL( (std::max<T>)(a, b), b);
+   BOOST_CHECK_EQUAL( (std::max<T>)(b, a), b);
+   BOOST_CHECK_EQUAL( (std::min<T>)(a, b + c), a);
+   BOOST_CHECK_EQUAL( (std::min<T>)(b + c, a), a);
+   BOOST_CHECK_EQUAL( (std::min<T>)(a, c - b), 1);
+   BOOST_CHECK_EQUAL( (std::min<T>)(c - b, a), 1);
+   BOOST_CHECK_EQUAL( (std::max<T>)(a, b + c), 11);
+   BOOST_CHECK_EQUAL( (std::max<T>)(b + c, a), 11);
+   BOOST_CHECK_EQUAL( (std::max<T>)(a, c - b), a);
+   BOOST_CHECK_EQUAL( (std::max<T>)(c - b, a), a);
+   BOOST_CHECK_EQUAL( (std::min<T>)(a + b, b + c), 7);
+   BOOST_CHECK_EQUAL( (std::min<T>)(b + c, a + b), 7);
+   BOOST_CHECK_EQUAL( (std::max<T>)(a + b, b + c), 11);
+   BOOST_CHECK_EQUAL( (std::max<T>)(b + c, a + b), 11);
+   BOOST_CHECK_EQUAL( (std::min<T>)(a + b, c - a), 4);
+   BOOST_CHECK_EQUAL( (std::min<T>)(c - a, a + b), 4);
+   BOOST_CHECK_EQUAL( (std::max<T>)(a + b, c - a), 7);
+   BOOST_CHECK_EQUAL( (std::max<T>)(c - a, a + b), 7);
 
    long l1(2), l2(3), l3;
-   l3 = min(l1, l2) + max(l1, l2) + max<long>(l1, l2) + min<long>(l1, l2);
+   l3 = (std::min)(l1, l2) + (std::max)(l1, l2) + (std::max<long>)(l1, l2) + (std::min<long>)(l1, l2);
    BOOST_CHECK_EQUAL(l3, 10);
 
 #endif
-
-
 }
 
+template <class T>
+const T& self(const T& a) { return a; }
+
 template <class Real>
 void test()
 {
@@ -2677,8 +2683,8 @@ void test()
    test_mixed<Real, double>(tag);
    test_mixed<Real, long double>(tag);
 
-   typedef typename related_type<Real>::type related_type;
-   boost::mpl::bool_<boost::multiprecision::is_number<Real>::value && !boost::is_same<related_type, Real>::value>  tag2;
+   typedef typename related_type<Real>::type                                                                      related_type;
+   boost::mpl::bool_<boost::multiprecision::is_number<Real>::value && !boost::is_same<related_type, Real>::value> tag2;
 
    test_mixed<Real, related_type>(tag2);
 
@@ -2704,145 +2710,145 @@ void test()
    Real b(64);
    Real c(500);
    Real d(1024);
-   BOOST_CHECK_EQUAL(a + b ,  72);
+   BOOST_CHECK_EQUAL(a + b, 72);
    a += b;
-   BOOST_CHECK_EQUAL(a ,  72);
-   BOOST_CHECK_EQUAL(a - b ,  8);
+   BOOST_CHECK_EQUAL(a, 72);
+   BOOST_CHECK_EQUAL(a - b, 8);
    a -= b;
-   BOOST_CHECK_EQUAL(a ,  8);
-   BOOST_CHECK_EQUAL(a * b ,  8*64L);
+   BOOST_CHECK_EQUAL(a, 8);
+   BOOST_CHECK_EQUAL(a * b, 8 * 64L);
    a *= b;
-   BOOST_CHECK_EQUAL(a ,  8*64L);
-   BOOST_CHECK_EQUAL(a / b ,  8);
+   BOOST_CHECK_EQUAL(a, 8 * 64L);
+   BOOST_CHECK_EQUAL(a / b, 8);
    a /= b;
-   BOOST_CHECK_EQUAL(a ,  8);
+   BOOST_CHECK_EQUAL(a, 8);
    Real ac(a);
-   BOOST_CHECK_EQUAL(ac ,  a);
+   BOOST_CHECK_EQUAL(ac, a);
    ac = a * c;
-   BOOST_CHECK_EQUAL(ac ,  8*500L);
-   ac = 8*500L;
+   BOOST_CHECK_EQUAL(ac, 8 * 500L);
+   ac = 8 * 500L;
    ac = ac + b + c;
-   BOOST_CHECK_EQUAL(ac ,  8*500L+64+500);
+   BOOST_CHECK_EQUAL(ac, 8 * 500L + 64 + 500);
    ac = a;
    ac = b + c + ac;
-   BOOST_CHECK_EQUAL(ac ,  8+64+500);
+   BOOST_CHECK_EQUAL(ac, 8 + 64 + 500);
    ac = ac - b + c;
-   BOOST_CHECK_EQUAL(ac ,  8+64+500-64+500);
+   BOOST_CHECK_EQUAL(ac, 8 + 64 + 500 - 64 + 500);
    ac = a;
    ac = b + c - ac;
-   BOOST_CHECK_EQUAL(ac ,  -8+64+500);
+   BOOST_CHECK_EQUAL(ac, -8 + 64 + 500);
    ac = a;
    ac = ac * b;
-   BOOST_CHECK_EQUAL(ac ,  8*64);
+   BOOST_CHECK_EQUAL(ac, 8 * 64);
    ac = a;
    ac *= b * ac;
-   BOOST_CHECK_EQUAL(ac ,  8*8*64);
+   BOOST_CHECK_EQUAL(ac, 8 * 8 * 64);
    ac = b;
    ac = ac / a;
-   BOOST_CHECK_EQUAL(ac ,  64/8);
+   BOOST_CHECK_EQUAL(ac, 64 / 8);
    ac = b;
    ac /= ac / a;
-   BOOST_CHECK_EQUAL(ac ,  64 / (64/8));
+   BOOST_CHECK_EQUAL(ac, 64 / (64 / 8));
    ac = a;
    ac = b + ac * a;
-   BOOST_CHECK_EQUAL(ac ,  64 * 2);
+   BOOST_CHECK_EQUAL(ac, 64 * 2);
    ac = a;
    ac = b - ac * a;
-   BOOST_CHECK_EQUAL(ac ,  0);
+   BOOST_CHECK_EQUAL(ac, 0);
    ac = a;
    ac = b * (ac + a);
-   BOOST_CHECK_EQUAL(ac ,  64 * (16));
+   BOOST_CHECK_EQUAL(ac, 64 * (16));
    ac = a;
    ac = b / (ac * 1);
-   BOOST_CHECK_EQUAL(ac ,  64 / 8);
+   BOOST_CHECK_EQUAL(ac, 64 / 8);
    ac = a;
    ac = ac + b;
-   BOOST_CHECK_EQUAL(ac ,  8 + 64);
+   BOOST_CHECK_EQUAL(ac, 8 + 64);
    ac = a;
    ac = a + ac;
-   BOOST_CHECK_EQUAL(ac ,  16);
+   BOOST_CHECK_EQUAL(ac, 16);
    ac = a;
    ac = a - ac;
-   BOOST_CHECK_EQUAL(ac ,  0);
+   BOOST_CHECK_EQUAL(ac, 0);
    ac = a;
    ac += a + b;
-   BOOST_CHECK_EQUAL(ac ,  80);
+   BOOST_CHECK_EQUAL(ac, 80);
    ac = a;
    ac += b + a;
-   BOOST_CHECK_EQUAL(ac ,  80);
+   BOOST_CHECK_EQUAL(ac, 80);
    ac = +a;
-   BOOST_CHECK_EQUAL(ac ,  8);
+   BOOST_CHECK_EQUAL(ac, 8);
    ac = 8;
    ac = a * ac;
-   BOOST_CHECK_EQUAL(ac ,  8*8);
+   BOOST_CHECK_EQUAL(ac, 8 * 8);
    ac = a;
    ac = a;
    ac += +a;
-   BOOST_CHECK_EQUAL(ac ,  16);
+   BOOST_CHECK_EQUAL(ac, 16);
    ac = a;
    ac += b - a;
-   BOOST_CHECK_EQUAL(ac ,  8 + 64-8);
+   BOOST_CHECK_EQUAL(ac, 8 + 64 - 8);
    ac = a;
-   ac += b*c;
-   BOOST_CHECK_EQUAL(ac ,  8 + 64 * 500);
+   ac += b * c;
+   BOOST_CHECK_EQUAL(ac, 8 + 64 * 500);
    ac = a;
    ac = a;
    ac -= +a;
-   BOOST_CHECK_EQUAL(ac ,  0);
+   BOOST_CHECK_EQUAL(ac, 0);
    ac = a;
-   if(std::numeric_limits<Real>::is_signed || is_twos_complement_integer<Real>::value)
+   if (std::numeric_limits<Real>::is_signed || is_twos_complement_integer<Real>::value)
    {
       ac = a;
       ac -= c - b;
-      BOOST_CHECK_EQUAL(ac ,  8 - (500-64));
+      BOOST_CHECK_EQUAL(ac, 8 - (500 - 64));
       ac = a;
-      ac -= b*c;
-      BOOST_CHECK_EQUAL(ac ,  8 - 500*64);
+      ac -= b * c;
+      BOOST_CHECK_EQUAL(ac, 8 - 500 * 64);
    }
    ac = a;
    ac += ac * b;
-   BOOST_CHECK_EQUAL(ac ,  8 + 8 * 64);
-   if(std::numeric_limits<Real>::is_signed || is_twos_complement_integer<Real>::value)
+   BOOST_CHECK_EQUAL(ac, 8 + 8 * 64);
+   if (std::numeric_limits<Real>::is_signed || is_twos_complement_integer<Real>::value)
    {
       ac = a;
       ac -= ac * b;
-      BOOST_CHECK_EQUAL(ac ,  8 - 8 * 64);
+      BOOST_CHECK_EQUAL(ac, 8 - 8 * 64);
    }
    ac = a * 8;
    ac *= +a;
-   BOOST_CHECK_EQUAL(ac ,  64 * 8);
+   BOOST_CHECK_EQUAL(ac, 64 * 8);
    ac = a;
    ac *= b * c;
-   BOOST_CHECK_EQUAL(ac ,  8 * 64 * 500);
+   BOOST_CHECK_EQUAL(ac, 8 * 64 * 500);
    ac = a;
    ac *= b / a;
-   BOOST_CHECK_EQUAL(ac ,  8 * 64 / 8);
+   BOOST_CHECK_EQUAL(ac, 8 * 64 / 8);
    ac = a;
    ac *= b + c;
-   BOOST_CHECK_EQUAL(ac ,  8 * (64 + 500));
+   BOOST_CHECK_EQUAL(ac, 8 * (64 + 500));
    ac = b;
    ac /= +a;
-   BOOST_CHECK_EQUAL(ac ,  8);
+   BOOST_CHECK_EQUAL(ac, 8);
    ac = b;
    ac /= b / a;
-   BOOST_CHECK_EQUAL(ac ,  64 / (64/8));
+   BOOST_CHECK_EQUAL(ac, 64 / (64 / 8));
    ac = b;
    ac /= a + Real(0);
-   BOOST_CHECK_EQUAL(ac ,  8);
+   BOOST_CHECK_EQUAL(ac, 8);
    //
    // simple tests with immediate values, these calls can be optimised in many backends:
    //
    ac = a + b;
-   BOOST_CHECK_EQUAL(ac ,  72);
+   BOOST_CHECK_EQUAL(ac, 72);
    ac = a + +b;
-   BOOST_CHECK_EQUAL(ac ,  72);
+   BOOST_CHECK_EQUAL(ac, 72);
    ac = +a + b;
-   BOOST_CHECK_EQUAL(ac ,  72);
+   BOOST_CHECK_EQUAL(ac, 72);
    ac = +a + +b;
-   BOOST_CHECK_EQUAL(ac ,  72);
+   BOOST_CHECK_EQUAL(ac, 72);
    ac = a;
    ac = b / ac;
-   BOOST_CHECK_EQUAL(ac ,  b / a);
+   BOOST_CHECK_EQUAL(ac, b / a);
    //
    // Comparisons:
    //
@@ -2862,64 +2868,64 @@ void test()
    b = 2;
    ss << a;
    ss >> c;
-   BOOST_CHECK_EQUAL(a ,  c);
+   BOOST_CHECK_EQUAL(a, c);
    ss.clear();
    ss << a + b;
    ss >> c;
-   BOOST_CHECK_EQUAL(c ,  22);
-   BOOST_CHECK_EQUAL(c ,  a + b);
+   BOOST_CHECK_EQUAL(c, 22);
+   BOOST_CHECK_EQUAL(c, a + b);
    //
    // More cases for complete code coverage:
    //
    a = 20;
    b = 30;
    swap(a, b);
-   BOOST_CHECK_EQUAL(a ,  30);
-   BOOST_CHECK_EQUAL(b ,  20);
+   BOOST_CHECK_EQUAL(a, 30);
+   BOOST_CHECK_EQUAL(b, 20);
    a = 20;
    b = 30;
    std::swap(a, b);
-   BOOST_CHECK_EQUAL(a ,  30);
-   BOOST_CHECK_EQUAL(b ,  20);
+   BOOST_CHECK_EQUAL(a, 30);
+   BOOST_CHECK_EQUAL(b, 20);
    a = 20;
    b = 30;
    a = a + b * 2;
-   BOOST_CHECK_EQUAL(a ,  20 + 30 * 2);
+   BOOST_CHECK_EQUAL(a, 20 + 30 * 2);
    a = 100;
    a = a - b * 2;
-   BOOST_CHECK_EQUAL(a ,  100 - 30 * 2);
+   BOOST_CHECK_EQUAL(a, 100 - 30 * 2);
    a = 20;
    a = a * (b + 2);
-   BOOST_CHECK_EQUAL(a ,  20 * (32));
+   BOOST_CHECK_EQUAL(a, 20 * (32));
    a = 20;
    a = (b + 2) * a;
-   BOOST_CHECK_EQUAL(a ,  20 * (32));
+   BOOST_CHECK_EQUAL(a, 20 * (32));
    a = 90;
    b = 2;
    a = a / (b + 0);
-   BOOST_CHECK_EQUAL(a ,  45);
+   BOOST_CHECK_EQUAL(a, 45);
    a = 20;
    b = 30;
    c = (a * b) + 22;
-   BOOST_CHECK_EQUAL(c ,  20 * 30 + 22);
+   BOOST_CHECK_EQUAL(c, 20 * 30 + 22);
    c = 22 + (a * b);
-   BOOST_CHECK_EQUAL(c ,  20 * 30 + 22);
-   c = 10;
+   BOOST_CHECK_EQUAL(c, 20 * 30 + 22);
+   c  = 10;
    ac = a + b * c;
-   BOOST_CHECK_EQUAL(ac ,  20 + 30 * 10);
+   BOOST_CHECK_EQUAL(ac, 20 + 30 * 10);
    ac = b * c + a;
-   BOOST_CHECK_EQUAL(ac ,  20 + 30 * 10);
+   BOOST_CHECK_EQUAL(ac, 20 + 30 * 10);
    a = a + b * c;
-   BOOST_CHECK_EQUAL(a ,  20 + 30 * 10);
+   BOOST_CHECK_EQUAL(a, 20 + 30 * 10);
    a = 20;
    b = a + b * c;
-   BOOST_CHECK_EQUAL(b ,  20 + 30 * 10);
+   BOOST_CHECK_EQUAL(b, 20 + 30 * 10);
    b = 30;
    c = a + b * c;
-   BOOST_CHECK_EQUAL(c ,  20 + 30 * 10);
+   BOOST_CHECK_EQUAL(c, 20 + 30 * 10);
    c = 10;
    c = a + b / c;
-   BOOST_CHECK_EQUAL(c ,  20 + 30 / 10);
+   BOOST_CHECK_EQUAL(c, 20 + 30 / 10);
 
    //
    // Test conditionals:
@@ -2933,7 +2939,7 @@ void test()
    // Test hashing:
    //
    boost::hash<Real> hasher;
-   std::size_t s = hasher(a);
+   std::size_t       s = hasher(a);
    BOOST_CHECK_NE(s, 0);
 #ifndef BOOST_NO_CXX11_HDR_FUNCTIONAL
    std::hash<Real> hasher2;
@@ -2949,7 +2955,7 @@ void test()
    BOOST_CHECK_EQUAL(m, 20);
    // Move from already moved from object:
    Real m2(static_cast<Real&&>(a));
-   // assign from moved from object 
+   // assign from moved from object
    // (may result in "a" being left in valid state as implementation artifact):
    c = static_cast<Real&&>(a);
    // assignment to moved-from objects:
@@ -2965,18 +2971,18 @@ void test()
    // string and string_view:
    //
    {
-      std::string s("2");
-      Real x(s);
+      std::string s1("2");
+      Real        x(s1);
       BOOST_CHECK_EQUAL(x, 2);
-      s = "3";
-      x.assign(s);
+      s1 = "3";
+      x.assign(s1);
       BOOST_CHECK_EQUAL(x, 3);
 #ifndef BOOST_NO_CXX17_HDR_STRING_VIEW
-      s = "20";
-      std::string_view v(s.c_str(), 1);
-      Real y(v);
+      s1 = "20";
+      std::string_view v(s1.c_str(), 1);
+      Real             y(v);
       BOOST_CHECK_EQUAL(y, 2);
-      std::string_view v2(s.c_str(), 2);
+      std::string_view v2(s1.c_str(), 2);
       y.assign(v2);
       BOOST_CHECK_EQUAL(y, 20);
 #endif
@@ -2986,7 +2992,7 @@ void test()
    // Bug cases, self assignment first:
    //
    a = 20;
-   a = a;
+   a = self(a);
    BOOST_CHECK_EQUAL(a, 20);
 
    a = 2;