mirror of
https://github.com/boostorg/multiprecision.git
synced 2026-02-14 12:52:17 +00:00
126 lines
4.6 KiB
C++
126 lines
4.6 KiB
C++
///////////////////////////////////////////////////////////////
|
|
// Copyright 2012 John Maddock. Distributed under the Boost
|
|
// Software License, Version 1.0. (See accompanying file
|
|
// LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_
|
|
|
|
//
|
|
// Compare arithmetic results using fixed_int to GMP results.
|
|
//
|
|
|
|
#ifdef _MSC_VER
|
|
# define _SCL_SECURE_NO_WARNINGS
|
|
#endif
|
|
|
|
#include <boost/multiprecision/gmp.hpp>
|
|
#include <boost/multiprecision/fixed_int.hpp>
|
|
#include <boost/random/mersenne_twister.hpp>
|
|
#include <boost/random/uniform_int.hpp>
|
|
#include "test.hpp"
|
|
|
|
template <class T>
|
|
T generate_random(unsigned bits_wanted)
|
|
{
|
|
static boost::random::mt19937 gen;
|
|
typedef boost::random::mt19937::result_type random_type;
|
|
|
|
T max_val;
|
|
unsigned 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;
|
|
}
|
|
else
|
|
{
|
|
max_val = T(1) << 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;
|
|
|
|
unsigned terms_needed = digits / bits_per_r_val + 1;
|
|
|
|
T val = 0;
|
|
for(unsigned i = 0; i < terms_needed; ++i)
|
|
{
|
|
val *= (gen.max)();
|
|
val += gen();
|
|
}
|
|
val %= max_val;
|
|
return val;
|
|
}
|
|
|
|
int main()
|
|
{
|
|
using namespace boost::multiprecision;
|
|
typedef mp_number<fixed_int<1024, true> > packed_type;
|
|
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);
|
|
mpz_int c = generate_random<mpz_int>(256);
|
|
mpz_int d = generate_random<mpz_int>(32);
|
|
|
|
packed_type a1 = a.str();
|
|
packed_type b1 = b.str();
|
|
packed_type c1 = c.str();
|
|
packed_type d1 = d.str();
|
|
|
|
BOOST_CHECK_EQUAL(a.str(), a1.str());
|
|
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(c * d).str(), packed_type(c1 * d1).str());
|
|
BOOST_CHECK_EQUAL(mpz_int(b * c).str(), packed_type(b1 * c1).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());
|
|
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());
|
|
|
|
if(last_error_count != boost::detail::test_errors())
|
|
{
|
|
last_error_count = boost::detail::test_errors();
|
|
std::cout << std::hex << std::showbase;
|
|
|
|
std::cout << "a = " << a << std::endl;
|
|
std::cout << "a1 = " << a1 << std::endl;
|
|
std::cout << "b = " << b << std::endl;
|
|
std::cout << "b1 = " << b1 << std::endl;
|
|
std::cout << "c = " << c << std::endl;
|
|
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 << "-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();
|
|
}
|
|
|
|
|
|
|