123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101 |
- ///////////////////////////////////////////////////////////////
- // Copyright 2012 John Maddock. Distributed under the Boost
- // Software License, Version 1.0. (See accompanying file
- // LICENSE_1_0.txt or copy at https://www.boost.org/LICENSE_1_0.txt
- //
- // Compare arithmetic results using fixed_int to GMP results.
- //
- #ifdef _MSC_VER
- #define _SCL_SECURE_NO_WARNINGS
- #endif
- #include <boost/multiprecision/integer.hpp>
- #include "test.hpp"
- #ifdef BOOST_MSVC
- #pragma warning(disable : 4146)
- #endif
- template <class I, class H>
- void test()
- {
- using namespace boost::multiprecision;
- I i(0);
- #ifndef BOOST_NO_EXCEPTIONS
- BOOST_CHECK_THROW(lsb(i), std::range_error);
- #endif
- BOOST_CHECK(bit_test(bit_set(i, 0), 0));
- BOOST_CHECK_EQUAL(bit_set(i, 0), 1);
- BOOST_CHECK_EQUAL(bit_unset(i, 0), 0);
- BOOST_CHECK_EQUAL(bit_flip(bit_set(i, 0), 0), 0);
- unsigned max_index = (std::numeric_limits<I>::digits) - 1;
- BOOST_CHECK(bit_test(bit_set(i, max_index), max_index));
- BOOST_CHECK_EQUAL(bit_unset(i, max_index), 0);
- BOOST_CHECK_EQUAL(bit_flip(bit_set(i, max_index), max_index), 0);
- i = 0;
- bit_set(i, max_index);
- BOOST_CHECK_EQUAL(lsb(i), max_index);
- BOOST_CHECK_EQUAL(msb(i), max_index);
- bit_set(i, max_index / 2);
- BOOST_CHECK_EQUAL(lsb(i), max_index / 2);
- BOOST_CHECK_EQUAL(msb(i), max_index);
- #ifndef BOOST_NO_EXCEPTIONS
- if (std::numeric_limits<I>::is_signed)
- {
- i = static_cast<I>(-1);
- BOOST_CHECK_THROW(lsb(i), std::range_error);
- }
- #endif
- H mx = (std::numeric_limits<H>::max)();
- BOOST_CHECK_EQUAL(multiply(i, mx, mx), static_cast<I>(mx) * static_cast<I>(mx));
- BOOST_CHECK_EQUAL(add(i, mx, mx), static_cast<I>(mx) + static_cast<I>(mx));
- if (std::numeric_limits<I>::is_signed)
- {
- BOOST_CHECK_EQUAL(subtract(i, mx, static_cast<H>(-mx)), static_cast<I>(mx) - static_cast<I>(-mx));
- BOOST_CHECK_EQUAL(add(i, static_cast<H>(-mx), static_cast<H>(-mx)), static_cast<I>(-mx) + static_cast<I>(-mx));
- }
- i = (std::numeric_limits<I>::max)();
- I j = 12345;
- I r, q;
- divide_qr(i, j, q, r);
- BOOST_CHECK_EQUAL(q, i / j);
- BOOST_CHECK_EQUAL(r, i % j);
- BOOST_CHECK_EQUAL(integer_modulus(i, j), i % j);
- I p = 456;
- BOOST_CHECK_EQUAL(powm(i, p, j), pow(cpp_int(i), static_cast<unsigned>(p)) % j);
- for (I i = 0; i < (2 < 8) - 1; ++i)
- {
- I j = i * i;
- I s, r;
- s = sqrt(j, r);
- BOOST_CHECK_EQUAL(s, i);
- BOOST_CHECK(r == 0);
- j += 3;
- s = sqrt(i, r);
- BOOST_CHECK_EQUAL(s, i);
- BOOST_CHECK(r == 3);
- }
- }
- int main()
- {
- using namespace boost::multiprecision;
- test<boost::int16_t, boost::int8_t>();
- test<boost::int32_t, boost::int16_t>();
- test<boost::int64_t, boost::int32_t>();
- test<boost::uint16_t, boost::uint8_t>();
- test<boost::uint32_t, boost::uint16_t>();
- test<boost::uint64_t, boost::uint32_t>();
- return boost::report_errors();
- }
|