//
#ifdef _MSC_VER
-# define _SCL_SECURE_NO_WARNINGS
+#define _SCL_SECURE_NO_WARNINGS
#endif
#include <boost/multiprecision/gmp.hpp>
template <class T>
T generate_random(unsigned bits_wanted)
{
- static boost::random::mt19937 gen;
+ static boost::random::mt19937 gen;
typedef boost::random::mt19937::result_type random_type;
- T max_val;
+ T max_val;
unsigned digits;
- if(std::numeric_limits<T>::is_bounded && (bits_wanted == std::numeric_limits<T>::digits))
+ if (std::numeric_limits<T>::is_bounded && (bits_wanted == std::numeric_limits<T>::digits))
{
max_val = (std::numeric_limits<T>::max)();
- digits = std::numeric_limits<T>::digits;
+ digits = std::numeric_limits<T>::digits;
}
else
{
max_val = T(1) << bits_wanted;
- digits = bits_wanted;
+ digits = bits_wanted;
}
unsigned bits_per_r_val = std::numeric_limits<random_type>::digits - 1;
- while((random_type(1) << bits_per_r_val) > (gen.max)()) --bits_per_r_val;
+ while ((random_type(1) << bits_per_r_val) > (gen.max)())
+ --bits_per_r_val;
unsigned terms_needed = digits / bits_per_r_val + 1;
T val = 0;
- for(unsigned i = 0; i < terms_needed; ++i)
+ for (unsigned i = 0; i < terms_needed; ++i)
{
val *= (gen.max)();
val += gen();
{
using namespace boost::multiprecision;
typedef number<fixed_int<1024, true> > packed_type;
- unsigned last_error_count = 0;
- for(int i = 0; i < 1000; ++i)
+ unsigned last_error_count = 0;
+ for (int i = 0; i < 1000; ++i)
{
mpz_int a = generate_random<mpz_int>(1000);
mpz_int b = generate_random<mpz_int>(512);
BOOST_CHECK_EQUAL(b.str(), b1.str());
BOOST_CHECK_EQUAL(c.str(), c1.str());
BOOST_CHECK_EQUAL(d.str(), d1.str());
- BOOST_CHECK_EQUAL(mpz_int(a+b).str(), packed_type(a1 + b1).str());
- BOOST_CHECK_EQUAL(mpz_int(a-b).str(), packed_type(a1 - b1).str());
- BOOST_CHECK_EQUAL(mpz_int(mpz_int(-a)+b).str(), packed_type(packed_type(-a1) + b1).str());
- BOOST_CHECK_EQUAL(mpz_int(mpz_int(-a)-b).str(), packed_type(packed_type(-a1) - b1).str());
+ BOOST_CHECK_EQUAL(mpz_int(a + b).str(), packed_type(a1 + b1).str());
+ BOOST_CHECK_EQUAL(mpz_int(a - b).str(), packed_type(a1 - b1).str());
+ BOOST_CHECK_EQUAL(mpz_int(mpz_int(-a) + b).str(), packed_type(packed_type(-a1) + b1).str());
+ BOOST_CHECK_EQUAL(mpz_int(mpz_int(-a) - b).str(), packed_type(packed_type(-a1) - b1).str());
BOOST_CHECK_EQUAL(mpz_int(c * d).str(), packed_type(c1 * d1).str());
BOOST_CHECK_EQUAL(mpz_int(c * -d).str(), packed_type(c1 * -d1).str());
BOOST_CHECK_EQUAL(mpz_int(-c * d).str(), packed_type(-c1 * d1).str());
BOOST_CHECK_EQUAL(mpz_int(-a % b).str(), packed_type(-a1 % b1).str());
BOOST_CHECK_EQUAL(mpz_int(a % d).str(), packed_type(a1 % d1).str());
// bitwise ops:
- BOOST_CHECK_EQUAL(mpz_int(a|b).str(), packed_type(a1 | b1).str());
- BOOST_CHECK_EQUAL(mpz_int(a&b).str(), packed_type(a1 & b1).str());
- BOOST_CHECK_EQUAL(mpz_int(a^b).str(), packed_type(a1 ^ b1).str());
+ BOOST_CHECK_EQUAL(mpz_int(a | b).str(), packed_type(a1 | b1).str());
+ BOOST_CHECK_EQUAL(mpz_int(a & b).str(), packed_type(a1 & b1).str());
+ BOOST_CHECK_EQUAL(mpz_int(a ^ b).str(), packed_type(a1 ^ b1).str());
// Now check operations involving integers:
BOOST_CHECK_EQUAL(mpz_int(a + si).str(), packed_type(a1 + si).str());
BOOST_CHECK_EQUAL(mpz_int(a + -si).str(), packed_type(a1 + -si).str());
BOOST_CHECK_EQUAL(mpz_int(a % si).str(), packed_type(a1 % si).str());
BOOST_CHECK_EQUAL(mpz_int(a % -si).str(), packed_type(a1 % -si).str());
BOOST_CHECK_EQUAL(mpz_int(-a % si).str(), packed_type(-a1 % si).str());
- BOOST_CHECK_EQUAL(mpz_int(a|si).str(), packed_type(a1 | si).str());
- BOOST_CHECK_EQUAL(mpz_int(a&si).str(), packed_type(a1 & si).str());
- BOOST_CHECK_EQUAL(mpz_int(a^si).str(), packed_type(a1 ^ si).str());
- BOOST_CHECK_EQUAL(mpz_int(si|a).str(), packed_type(si|a1).str());
- BOOST_CHECK_EQUAL(mpz_int(si&a).str(), packed_type(si&a1).str());
- BOOST_CHECK_EQUAL(mpz_int(si^a).str(), packed_type(si^a1).str());
+ BOOST_CHECK_EQUAL(mpz_int(a | si).str(), packed_type(a1 | si).str());
+ BOOST_CHECK_EQUAL(mpz_int(a & si).str(), packed_type(a1 & si).str());
+ BOOST_CHECK_EQUAL(mpz_int(a ^ si).str(), packed_type(a1 ^ si).str());
+ BOOST_CHECK_EQUAL(mpz_int(si | a).str(), packed_type(si | a1).str());
+ BOOST_CHECK_EQUAL(mpz_int(si & a).str(), packed_type(si & a1).str());
+ BOOST_CHECK_EQUAL(mpz_int(si ^ a).str(), packed_type(si ^ a1).str());
BOOST_CHECK_EQUAL(mpz_int(gcd(a, b)).str(), packed_type(gcd(a1, b1)).str());
BOOST_CHECK_EQUAL(mpz_int(lcm(c, d)).str(), packed_type(lcm(c1, d1)).str());
- if(last_error_count != boost::detail::test_errors())
+ if (last_error_count != boost::detail::test_errors())
{
last_error_count = boost::detail::test_errors();
std::cout << std::hex << std::showbase;
std::cout << "c1 = " << c1 << std::endl;
std::cout << "d = " << d << std::endl;
std::cout << "d1 = " << d1 << std::endl;
- std::cout << "a + b = " << a+b << std::endl;
- std::cout << "a1 + b1 = " << a1+b1 << std::endl;
+ std::cout << "a + b = " << a + b << std::endl;
+ std::cout << "a1 + b1 = " << a1 + b1 << std::endl;
std::cout << std::dec;
- std::cout << "a - b = " << a-b << std::endl;
- std::cout << "a1 - b1 = " << a1-b1 << std::endl;
- std::cout << "-a + b = " << mpz_int(-a)+b << std::endl;
- std::cout << "-a1 + b1 = " << packed_type(-a1)+b1 << std::endl;
- std::cout << "-a - b = " << mpz_int(-a)-b << std::endl;
- std::cout << "-a1 - b1 = " << packed_type(-a1)-b1 << std::endl;
- std::cout << "c*d = " << c*d << std::endl;
- std::cout << "c1*d1 = " << c1*d1 << std::endl;
- std::cout << "b*c = " << b*c << std::endl;
- std::cout << "b1*c1 = " << b1*c1 << std::endl;
- std::cout << "a/b = " << a/b << std::endl;
- std::cout << "a1/b1 = " << a1/b1 << std::endl;
- std::cout << "a/d = " << a/d << std::endl;
- std::cout << "a1/d1 = " << a1/d1 << std::endl;
- std::cout << "a%b = " << a%b << std::endl;
- std::cout << "a1%b1 = " << a1%b1 << std::endl;
- std::cout << "a%d = " << a%d << std::endl;
- std::cout << "a1%d1 = " << a1%d1 << std::endl;
+ std::cout << "a - b = " << a - b << std::endl;
+ std::cout << "a1 - b1 = " << a1 - b1 << std::endl;
+ std::cout << "-a + b = " << mpz_int(-a) + b << std::endl;
+ std::cout << "-a1 + b1 = " << packed_type(-a1) + b1 << std::endl;
+ std::cout << "-a - b = " << mpz_int(-a) - b << std::endl;
+ std::cout << "-a1 - b1 = " << packed_type(-a1) - b1 << std::endl;
+ std::cout << "c*d = " << c * d << std::endl;
+ std::cout << "c1*d1 = " << c1 * d1 << std::endl;
+ std::cout << "b*c = " << b * c << std::endl;
+ std::cout << "b1*c1 = " << b1 * c1 << std::endl;
+ std::cout << "a/b = " << a / b << std::endl;
+ std::cout << "a1/b1 = " << a1 / b1 << std::endl;
+ std::cout << "a/d = " << a / d << std::endl;
+ std::cout << "a1/d1 = " << a1 / d1 << std::endl;
+ std::cout << "a%b = " << a % b << std::endl;
+ std::cout << "a1%b1 = " << a1 % b1 << std::endl;
+ std::cout << "a%d = " << a % d << std::endl;
+ std::cout << "a1%d1 = " << a1 % d1 << std::endl;
}
}
return boost::report_errors();
}
-
-
-