boost/multiprecision
2
0
Fork 0
mirror of https://github.com/boostorg/multiprecision.git synced 2026-01-19 04:22:11 +00:00
Code Issues Projects Releases Wiki Activity

First cut at new rational_adaptor.

Browse Source
This commit is contained in:
jzmaddock
2021-09-02 10:53:35 +01:00
parent 3094249d7c
commit e6d136bc50
4 changed files with 1008 additions and 201 deletions
Download Patch File Download Diff File Expand all files Collapse all files

7
include/boost/multiprecision/cpp_int/misc.hpp
View File

@@ -416,9 +416,12 @@ BOOST_MP_FORCEINLINE BOOST_MP_CXX14_CONSTEXPR limb_type eval_gcd(limb_type u, li
inline BOOST_MP_CXX14_CONSTEXPR double_limb_type eval_gcd(double_limb_type u, double_limb_type v)
{
#if (__cpp_lib_gcd_lcm >= 201606L) && (!defined(BOOST_HAS_INT128) || !defined(__STRICT_ANSI__))
#if 0 && (__cpp_lib_gcd_lcm >= 201606L) && (!defined(BOOST_HAS_INT128) || !defined(__STRICT_ANSI__))
return std::gcd(u, v);
#else
if (u == 0)
return v;
unsigned shift = boost::multiprecision::detail::find_lsb(u | v);
u >>= boost::multiprecision::detail::find_lsb(u);
do
@@ -1110,6 +1113,7 @@ BOOST_MP_FORCEINLINE BOOST_MP_CXX14_CONSTEXPR typename std::enable_if<is_trivial
eval_gcd(cpp_int_backend<MinBits1, MaxBits1, SignType1, Checked1, Allocator1>& result, const cpp_int_backend<MinBits1, MaxBits1, SignType1, Checked1, Allocator1>& a, const cpp_int_backend<MinBits1, MaxBits1, SignType1, Checked1, Allocator1>& b) noexcept
{
*result.limbs() = boost::integer::gcd(*a.limbs(), *b.limbs());
result.sign(false);
}
// This one is only enabled for unchecked cpp_int's, for checked int's we need the checking in the default version:
template <unsigned MinBits1, unsigned MaxBits1, cpp_integer_type SignType1, cpp_int_check_type Checked1, class Allocator1>
@@ -1118,6 +1122,7 @@ eval_lcm(cpp_int_backend<MinBits1, MaxBits1, SignType1, Checked1, Allocator1>& r
{
*result.limbs() = boost::integer::lcm(*a.limbs(), *b.limbs());
result.normalize(); // result may overflow the specified number of bits
result.sign(false);
}
inline void conversion_overflow(const std::integral_constant<int, checked>&)

1198
include/boost/multiprecision/rational_adaptor.hpp
View File

File diff suppressed because it is too large Load Diff

1
test/test_arithmetic_cpp_int_br.cpp
View File

@@ -6,6 +6,7 @@
#define BOOST_MP_NOT_TESTING_NUMBER
#include <boost/multiprecision/cpp_int.hpp>
#include <boost/rational.hpp>
#include "test_arithmetic.hpp"

3
test/test_rational_io.cpp
View File

@@ -25,6 +25,7 @@
#if defined(TEST_MPQ)
#include <boost/multiprecision/gmp.hpp>
#include <boost/multiprecision/rational_adaptor.hpp>
#endif
#if defined(TEST_TOMMATH)
#include <boost/multiprecision/tommath.hpp>
@@ -36,7 +37,7 @@
#include <boost/algorithm/string/case_conv.hpp>
#include <boost/random/mersenne_twister.hpp>
#include <boost/random/uniform_int.hpp>
#include <boost/multiprecision/rational_adaptor.hpp>
#include <boost/rational.hpp>
#include "test.hpp"
#include <iostream>
#include <iomanip>