-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy pathdouble_integer.hpp
177 lines (144 loc) · 8.27 KB
/
double_integer.hpp
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
#ifndef DOUBLE_INTEGER_HPP
#define DOUBLE_INTEGER_HPP
/* Double Integer template
// double_integer.hpp
//
// Author: Jan Ringos, http://Tringi.TrimCore.cz, tringi@trimcore.cz
// Version: 2.3
// Description: Implements common functionality to double the precission
// of an little-endian two-complement integer type, both signed
// and unsigned.
//
// License:
// This software is provided 'as-is', without any express or implied warranty.
// In no event will the author be held liable for any damages arising from the
// use of this software.
//
// Permission is granted to anyone to use this software for any purpose,
// including commercial applications, and to alter it and redistribute it
// freely, subject to the following restrictions:
//
// 1. The origin of this software must not be misrepresented; you must not
// claim that you wrote the original software. If you use this software
// in a product, an acknowledgment in the product documentation would
// be appreciated but is not required.
//
// 2. Altered source versions must be plainly marked as such, and must not
// be misrepresented as being the original software.
//
// 3. This notice may not be removed or altered from any source
// distribution.
//
// Changelog:
// 29.01.2011 - initial version
// 28.05.2013 - now requires C++11 (to enable singed/unsigned conversions)
// - fixed conversion from negative floating point
// - C style fmod, fmodf, fmodl replaced by std::fmod
// 20.03.2015 - completed C++11 release
// 29.09.2023 - removed deprecated numeric_limits denorm_absent/has_denorm_loss
// - silenced several comparison warnings
// 30.05.2024 - added explicit conversion to LO and double, tweaked print
// 01.06.2024 - fixed numeric_limits::min, and modernized the numeric_limits
*/
#include <type_traits>
#include <cstddef>
#include <cstdio>
#include <limits>
// double_integer <LO, HI>
// - unlike the intrinsic types, this one is always zero initialized
// - arguments: LO - always an unsigned integer type
// HI - defaults to LO, type of the upper half
// - notes: - use of types with different width for LO and HI is unsupported
// - to implement a signed type, use a signed integer as HI type
template <typename LO, typename HI = LO>
class double_integer {
// check assumptions
static_assert (sizeof (LO) == sizeof (HI), "both types must have the same size");
static_assert (!std::numeric_limits <LO> ::is_signed, "type of lower half must be unsigned");
public:
// the only two data members
LO lo;
HI hi;
// constructors
constexpr double_integer () : lo (LO ()), hi (HI ()) {};
constexpr double_integer (const double_integer & a) : lo (a.lo), hi (a.hi) {};
constexpr double_integer (const LO & _lo, const HI & _hi) : lo (_lo), hi (_hi) {};
double_integer (signed char);
double_integer (signed short);
double_integer (signed int);
double_integer (signed long);
double_integer (signed long long);
double_integer (unsigned char);
double_integer (unsigned short);
double_integer (unsigned int);
double_integer (unsigned long);
double_integer (unsigned long long);
double_integer (float);
double_integer (double);
double_integer (long double);
// string constructor
// - skips std::isspace characters at string beginning
// - when radix is zero, is detected by following rules:
// "0" -> 8, "0x" -> 16, "0b" -> 2, otherwise 10
explicit double_integer (const char *, unsigned int radix = 0);
// assignment operators
// - only between double_integers
// and unsigned version of this one (if this one is signed)
double_integer & operator = (const double_integer &);
template <class = std::enable_if <!std::is_same <LO, HI>::value>>
double_integer & operator = (const double_integer <LO, LO> &);
// integer operators
constexpr bool operator ! () const;
constexpr explicit operator bool () const;
constexpr explicit operator LO () const;
constexpr explicit operator double () const;
constexpr operator double_integer <LO, LO> () const;
double_integer operator + () const { return *this; };
double_integer operator - () const;
double_integer operator ~ () const;
double_integer operator ++ (int);
double_integer operator -- (int);
double_integer & operator ++ ();
double_integer & operator -- ();
double_integer & operator += (const double_integer &);
double_integer & operator -= (const double_integer &);
double_integer & operator *= (const double_integer &);
double_integer & operator /= (const double_integer &);
double_integer & operator %= (const double_integer &);
double_integer & operator |= (const double_integer &);
double_integer & operator &= (const double_integer &);
double_integer & operator ^= (const double_integer &);
double_integer & operator >>= (unsigned int);
double_integer & operator <<= (unsigned int);
// division
// - returns result, second parameters (if provided) receives remainder
// - throws std::logic_error on division by zero
double_integer div (const double_integer &, double_integer * = 0) const;
// rendering
const char * print (char * buffer, std::size_t length,
unsigned int radix = 10u) const;
const char * c_str (unsigned int radix = 10u) const; // thread unsafe!
// internal implementation details and helpers follow
private:
template <typename T> static T mask (unsigned int);
};
template <typename LO, typename HI> bool operator < (const double_integer <LO, HI> &, const double_integer <LO, HI> &);
template <typename LO, typename HI> bool operator > (const double_integer <LO, HI> &, const double_integer <LO, HI> &);
template <typename LO, typename HI> bool operator <= (const double_integer <LO, HI> &, const double_integer <LO, HI> &);
template <typename LO, typename HI> bool operator <= (const double_integer <LO, HI> &, const double_integer <LO, HI> &);
template <typename LO, typename HI> bool operator == (const double_integer <LO, HI> &, const double_integer <LO, HI> &);
template <typename LO, typename HI> bool operator != (const double_integer <LO, HI> &, const double_integer <LO, HI> &);
template <typename LO, typename HI> bool operator || (const double_integer <LO, HI> &, const double_integer <LO, HI> &);
template <typename LO, typename HI> bool operator && (const double_integer <LO, HI> &, const double_integer <LO, HI> &);
template <typename LO, typename HI> double_integer <LO, HI> operator + (const double_integer <LO, HI> &, const double_integer <LO, HI> &);
template <typename LO, typename HI> double_integer <LO, HI> operator - (const double_integer <LO, HI> &, const double_integer <LO, HI> &);
template <typename LO, typename HI> double_integer <LO, HI> operator * (const double_integer <LO, HI> &, const double_integer <LO, HI> &);
template <typename LO, typename HI> double_integer <LO, HI> operator / (const double_integer <LO, HI> &, const double_integer <LO, HI> &);
template <typename LO, typename HI> double_integer <LO, HI> operator % (const double_integer <LO, HI> &, const double_integer <LO, HI> &);
template <typename LO, typename HI> double_integer <LO, HI> operator & (const double_integer <LO, HI> &, const double_integer <LO, HI> &);
template <typename LO, typename HI> double_integer <LO, HI> operator | (const double_integer <LO, HI> &, const double_integer <LO, HI> &);
template <typename LO, typename HI> double_integer <LO, HI> operator ^ (const double_integer <LO, HI> &, const double_integer <LO, HI> &);
template <typename LO, typename HI> double_integer <LO, HI> operator << (const double_integer <LO, HI> &, unsigned int);
template <typename LO, typename HI> double_integer <LO, HI> operator >> (const double_integer <LO, HI> &, unsigned int);
#include "double_integer.tcc"
#endif