isotree 0.2.2 → 0.3.0

This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
Files changed (151) hide show
  1. checksums.yaml +4 -4
  2. data/CHANGELOG.md +8 -1
  3. data/LICENSE.txt +2 -2
  4. data/README.md +32 -14
  5. data/ext/isotree/ext.cpp +144 -31
  6. data/ext/isotree/extconf.rb +7 -7
  7. data/lib/isotree/isolation_forest.rb +110 -30
  8. data/lib/isotree/version.rb +1 -1
  9. data/vendor/isotree/LICENSE +1 -1
  10. data/vendor/isotree/README.md +165 -27
  11. data/vendor/isotree/include/isotree.hpp +2111 -0
  12. data/vendor/isotree/include/isotree_oop.hpp +394 -0
  13. data/vendor/isotree/inst/COPYRIGHTS +62 -0
  14. data/vendor/isotree/src/RcppExports.cpp +525 -52
  15. data/vendor/isotree/src/Rwrapper.cpp +1931 -268
  16. data/vendor/isotree/src/c_interface.cpp +953 -0
  17. data/vendor/isotree/src/crit.hpp +4232 -0
  18. data/vendor/isotree/src/dist.hpp +1886 -0
  19. data/vendor/isotree/src/exp_depth_table.hpp +134 -0
  20. data/vendor/isotree/src/extended.hpp +1444 -0
  21. data/vendor/isotree/src/external_facing_generic.hpp +399 -0
  22. data/vendor/isotree/src/fit_model.hpp +2401 -0
  23. data/vendor/isotree/src/{dealloc.cpp → headers_joined.hpp} +38 -22
  24. data/vendor/isotree/src/helpers_iforest.hpp +813 -0
  25. data/vendor/isotree/src/{impute.cpp → impute.hpp} +353 -122
  26. data/vendor/isotree/src/indexer.cpp +515 -0
  27. data/vendor/isotree/src/instantiate_template_headers.cpp +118 -0
  28. data/vendor/isotree/src/instantiate_template_headers.hpp +240 -0
  29. data/vendor/isotree/src/isoforest.hpp +1659 -0
  30. data/vendor/isotree/src/isotree.hpp +1804 -392
  31. data/vendor/isotree/src/isotree_exportable.hpp +99 -0
  32. data/vendor/isotree/src/merge_models.cpp +159 -16
  33. data/vendor/isotree/src/mult.hpp +1321 -0
  34. data/vendor/isotree/src/oop_interface.cpp +842 -0
  35. data/vendor/isotree/src/oop_interface.hpp +278 -0
  36. data/vendor/isotree/src/other_helpers.hpp +219 -0
  37. data/vendor/isotree/src/predict.hpp +1932 -0
  38. data/vendor/isotree/src/python_helpers.hpp +134 -0
  39. data/vendor/isotree/src/ref_indexer.hpp +154 -0
  40. data/vendor/isotree/src/robinmap/LICENSE +21 -0
  41. data/vendor/isotree/src/robinmap/README.md +483 -0
  42. data/vendor/isotree/src/robinmap/include/tsl/robin_growth_policy.h +406 -0
  43. data/vendor/isotree/src/robinmap/include/tsl/robin_hash.h +1620 -0
  44. data/vendor/isotree/src/robinmap/include/tsl/robin_map.h +807 -0
  45. data/vendor/isotree/src/robinmap/include/tsl/robin_set.h +660 -0
  46. data/vendor/isotree/src/serialize.cpp +4300 -139
  47. data/vendor/isotree/src/sql.cpp +141 -59
  48. data/vendor/isotree/src/subset_models.cpp +174 -0
  49. data/vendor/isotree/src/utils.hpp +3808 -0
  50. data/vendor/isotree/src/xoshiro.hpp +467 -0
  51. data/vendor/isotree/src/ziggurat.hpp +405 -0
  52. metadata +38 -104
  53. data/vendor/cereal/LICENSE +0 -24
  54. data/vendor/cereal/README.md +0 -85
  55. data/vendor/cereal/include/cereal/access.hpp +0 -351
  56. data/vendor/cereal/include/cereal/archives/adapters.hpp +0 -163
  57. data/vendor/cereal/include/cereal/archives/binary.hpp +0 -169
  58. data/vendor/cereal/include/cereal/archives/json.hpp +0 -1019
  59. data/vendor/cereal/include/cereal/archives/portable_binary.hpp +0 -334
  60. data/vendor/cereal/include/cereal/archives/xml.hpp +0 -956
  61. data/vendor/cereal/include/cereal/cereal.hpp +0 -1089
  62. data/vendor/cereal/include/cereal/details/helpers.hpp +0 -422
  63. data/vendor/cereal/include/cereal/details/polymorphic_impl.hpp +0 -796
  64. data/vendor/cereal/include/cereal/details/polymorphic_impl_fwd.hpp +0 -65
  65. data/vendor/cereal/include/cereal/details/static_object.hpp +0 -127
  66. data/vendor/cereal/include/cereal/details/traits.hpp +0 -1411
  67. data/vendor/cereal/include/cereal/details/util.hpp +0 -84
  68. data/vendor/cereal/include/cereal/external/base64.hpp +0 -134
  69. data/vendor/cereal/include/cereal/external/rapidjson/allocators.h +0 -284
  70. data/vendor/cereal/include/cereal/external/rapidjson/cursorstreamwrapper.h +0 -78
  71. data/vendor/cereal/include/cereal/external/rapidjson/document.h +0 -2652
  72. data/vendor/cereal/include/cereal/external/rapidjson/encodedstream.h +0 -299
  73. data/vendor/cereal/include/cereal/external/rapidjson/encodings.h +0 -716
  74. data/vendor/cereal/include/cereal/external/rapidjson/error/en.h +0 -74
  75. data/vendor/cereal/include/cereal/external/rapidjson/error/error.h +0 -161
  76. data/vendor/cereal/include/cereal/external/rapidjson/filereadstream.h +0 -99
  77. data/vendor/cereal/include/cereal/external/rapidjson/filewritestream.h +0 -104
  78. data/vendor/cereal/include/cereal/external/rapidjson/fwd.h +0 -151
  79. data/vendor/cereal/include/cereal/external/rapidjson/internal/biginteger.h +0 -290
  80. data/vendor/cereal/include/cereal/external/rapidjson/internal/diyfp.h +0 -271
  81. data/vendor/cereal/include/cereal/external/rapidjson/internal/dtoa.h +0 -245
  82. data/vendor/cereal/include/cereal/external/rapidjson/internal/ieee754.h +0 -78
  83. data/vendor/cereal/include/cereal/external/rapidjson/internal/itoa.h +0 -308
  84. data/vendor/cereal/include/cereal/external/rapidjson/internal/meta.h +0 -186
  85. data/vendor/cereal/include/cereal/external/rapidjson/internal/pow10.h +0 -55
  86. data/vendor/cereal/include/cereal/external/rapidjson/internal/regex.h +0 -740
  87. data/vendor/cereal/include/cereal/external/rapidjson/internal/stack.h +0 -232
  88. data/vendor/cereal/include/cereal/external/rapidjson/internal/strfunc.h +0 -69
  89. data/vendor/cereal/include/cereal/external/rapidjson/internal/strtod.h +0 -290
  90. data/vendor/cereal/include/cereal/external/rapidjson/internal/swap.h +0 -46
  91. data/vendor/cereal/include/cereal/external/rapidjson/istreamwrapper.h +0 -128
  92. data/vendor/cereal/include/cereal/external/rapidjson/memorybuffer.h +0 -70
  93. data/vendor/cereal/include/cereal/external/rapidjson/memorystream.h +0 -71
  94. data/vendor/cereal/include/cereal/external/rapidjson/msinttypes/inttypes.h +0 -316
  95. data/vendor/cereal/include/cereal/external/rapidjson/msinttypes/stdint.h +0 -300
  96. data/vendor/cereal/include/cereal/external/rapidjson/ostreamwrapper.h +0 -81
  97. data/vendor/cereal/include/cereal/external/rapidjson/pointer.h +0 -1414
  98. data/vendor/cereal/include/cereal/external/rapidjson/prettywriter.h +0 -277
  99. data/vendor/cereal/include/cereal/external/rapidjson/rapidjson.h +0 -656
  100. data/vendor/cereal/include/cereal/external/rapidjson/reader.h +0 -2230
  101. data/vendor/cereal/include/cereal/external/rapidjson/schema.h +0 -2497
  102. data/vendor/cereal/include/cereal/external/rapidjson/stream.h +0 -223
  103. data/vendor/cereal/include/cereal/external/rapidjson/stringbuffer.h +0 -121
  104. data/vendor/cereal/include/cereal/external/rapidjson/writer.h +0 -709
  105. data/vendor/cereal/include/cereal/external/rapidxml/license.txt +0 -52
  106. data/vendor/cereal/include/cereal/external/rapidxml/manual.html +0 -406
  107. data/vendor/cereal/include/cereal/external/rapidxml/rapidxml.hpp +0 -2624
  108. data/vendor/cereal/include/cereal/external/rapidxml/rapidxml_iterators.hpp +0 -175
  109. data/vendor/cereal/include/cereal/external/rapidxml/rapidxml_print.hpp +0 -428
  110. data/vendor/cereal/include/cereal/external/rapidxml/rapidxml_utils.hpp +0 -123
  111. data/vendor/cereal/include/cereal/macros.hpp +0 -154
  112. data/vendor/cereal/include/cereal/specialize.hpp +0 -139
  113. data/vendor/cereal/include/cereal/types/array.hpp +0 -79
  114. data/vendor/cereal/include/cereal/types/atomic.hpp +0 -55
  115. data/vendor/cereal/include/cereal/types/base_class.hpp +0 -203
  116. data/vendor/cereal/include/cereal/types/bitset.hpp +0 -176
  117. data/vendor/cereal/include/cereal/types/boost_variant.hpp +0 -164
  118. data/vendor/cereal/include/cereal/types/chrono.hpp +0 -72
  119. data/vendor/cereal/include/cereal/types/common.hpp +0 -129
  120. data/vendor/cereal/include/cereal/types/complex.hpp +0 -56
  121. data/vendor/cereal/include/cereal/types/concepts/pair_associative_container.hpp +0 -73
  122. data/vendor/cereal/include/cereal/types/deque.hpp +0 -62
  123. data/vendor/cereal/include/cereal/types/forward_list.hpp +0 -68
  124. data/vendor/cereal/include/cereal/types/functional.hpp +0 -43
  125. data/vendor/cereal/include/cereal/types/list.hpp +0 -62
  126. data/vendor/cereal/include/cereal/types/map.hpp +0 -36
  127. data/vendor/cereal/include/cereal/types/memory.hpp +0 -425
  128. data/vendor/cereal/include/cereal/types/optional.hpp +0 -66
  129. data/vendor/cereal/include/cereal/types/polymorphic.hpp +0 -483
  130. data/vendor/cereal/include/cereal/types/queue.hpp +0 -132
  131. data/vendor/cereal/include/cereal/types/set.hpp +0 -103
  132. data/vendor/cereal/include/cereal/types/stack.hpp +0 -76
  133. data/vendor/cereal/include/cereal/types/string.hpp +0 -61
  134. data/vendor/cereal/include/cereal/types/tuple.hpp +0 -123
  135. data/vendor/cereal/include/cereal/types/unordered_map.hpp +0 -36
  136. data/vendor/cereal/include/cereal/types/unordered_set.hpp +0 -99
  137. data/vendor/cereal/include/cereal/types/utility.hpp +0 -47
  138. data/vendor/cereal/include/cereal/types/valarray.hpp +0 -89
  139. data/vendor/cereal/include/cereal/types/variant.hpp +0 -109
  140. data/vendor/cereal/include/cereal/types/vector.hpp +0 -112
  141. data/vendor/cereal/include/cereal/version.hpp +0 -52
  142. data/vendor/isotree/src/Makevars +0 -4
  143. data/vendor/isotree/src/crit.cpp +0 -912
  144. data/vendor/isotree/src/dist.cpp +0 -749
  145. data/vendor/isotree/src/extended.cpp +0 -790
  146. data/vendor/isotree/src/fit_model.cpp +0 -1090
  147. data/vendor/isotree/src/helpers_iforest.cpp +0 -324
  148. data/vendor/isotree/src/isoforest.cpp +0 -771
  149. data/vendor/isotree/src/mult.cpp +0 -607
  150. data/vendor/isotree/src/predict.cpp +0 -853
  151. data/vendor/isotree/src/utils.cpp +0 -1566
@@ -1,232 +0,0 @@
1
- // Tencent is pleased to support the open source community by making RapidJSON available.
2
- //
3
- // Copyright (C) 2015 THL A29 Limited, a Tencent company, and Milo Yip. All rights reserved.
4
- //
5
- // Licensed under the MIT License (the "License"); you may not use this file except
6
- // in compliance with the License. You may obtain a copy of the License at
7
- //
8
- // http://opensource.org/licenses/MIT
9
- //
10
- // Unless required by applicable law or agreed to in writing, software distributed
11
- // under the License is distributed on an "AS IS" BASIS, WITHOUT WARRANTIES OR
12
- // CONDITIONS OF ANY KIND, either express or implied. See the License for the
13
- // specific language governing permissions and limitations under the License.
14
-
15
- #ifndef CEREAL_RAPIDJSON_INTERNAL_STACK_H_
16
- #define CEREAL_RAPIDJSON_INTERNAL_STACK_H_
17
-
18
- #include "../allocators.h"
19
- #include "swap.h"
20
- #include <cstddef>
21
-
22
- #if defined(__clang__)
23
- CEREAL_RAPIDJSON_DIAG_PUSH
24
- CEREAL_RAPIDJSON_DIAG_OFF(c++98-compat)
25
- #endif
26
-
27
- CEREAL_RAPIDJSON_NAMESPACE_BEGIN
28
- namespace internal {
29
-
30
- ///////////////////////////////////////////////////////////////////////////////
31
- // Stack
32
-
33
- //! A type-unsafe stack for storing different types of data.
34
- /*! \tparam Allocator Allocator for allocating stack memory.
35
- */
36
- template <typename Allocator>
37
- class Stack {
38
- public:
39
- // Optimization note: Do not allocate memory for stack_ in constructor.
40
- // Do it lazily when first Push() -> Expand() -> Resize().
41
- Stack(Allocator* allocator, size_t stackCapacity) : allocator_(allocator), ownAllocator_(0), stack_(0), stackTop_(0), stackEnd_(0), initialCapacity_(stackCapacity) {
42
- }
43
-
44
- #if CEREAL_RAPIDJSON_HAS_CXX11_RVALUE_REFS
45
- Stack(Stack&& rhs)
46
- : allocator_(rhs.allocator_),
47
- ownAllocator_(rhs.ownAllocator_),
48
- stack_(rhs.stack_),
49
- stackTop_(rhs.stackTop_),
50
- stackEnd_(rhs.stackEnd_),
51
- initialCapacity_(rhs.initialCapacity_)
52
- {
53
- rhs.allocator_ = 0;
54
- rhs.ownAllocator_ = 0;
55
- rhs.stack_ = 0;
56
- rhs.stackTop_ = 0;
57
- rhs.stackEnd_ = 0;
58
- rhs.initialCapacity_ = 0;
59
- }
60
- #endif
61
-
62
- ~Stack() {
63
- Destroy();
64
- }
65
-
66
- #if CEREAL_RAPIDJSON_HAS_CXX11_RVALUE_REFS
67
- Stack& operator=(Stack&& rhs) {
68
- if (&rhs != this)
69
- {
70
- Destroy();
71
-
72
- allocator_ = rhs.allocator_;
73
- ownAllocator_ = rhs.ownAllocator_;
74
- stack_ = rhs.stack_;
75
- stackTop_ = rhs.stackTop_;
76
- stackEnd_ = rhs.stackEnd_;
77
- initialCapacity_ = rhs.initialCapacity_;
78
-
79
- rhs.allocator_ = 0;
80
- rhs.ownAllocator_ = 0;
81
- rhs.stack_ = 0;
82
- rhs.stackTop_ = 0;
83
- rhs.stackEnd_ = 0;
84
- rhs.initialCapacity_ = 0;
85
- }
86
- return *this;
87
- }
88
- #endif
89
-
90
- void Swap(Stack& rhs) CEREAL_RAPIDJSON_NOEXCEPT {
91
- internal::Swap(allocator_, rhs.allocator_);
92
- internal::Swap(ownAllocator_, rhs.ownAllocator_);
93
- internal::Swap(stack_, rhs.stack_);
94
- internal::Swap(stackTop_, rhs.stackTop_);
95
- internal::Swap(stackEnd_, rhs.stackEnd_);
96
- internal::Swap(initialCapacity_, rhs.initialCapacity_);
97
- }
98
-
99
- void Clear() { stackTop_ = stack_; }
100
-
101
- void ShrinkToFit() {
102
- if (Empty()) {
103
- // If the stack is empty, completely deallocate the memory.
104
- Allocator::Free(stack_); // NOLINT (+clang-analyzer-unix.Malloc)
105
- stack_ = 0;
106
- stackTop_ = 0;
107
- stackEnd_ = 0;
108
- }
109
- else
110
- Resize(GetSize());
111
- }
112
-
113
- // Optimization note: try to minimize the size of this function for force inline.
114
- // Expansion is run very infrequently, so it is moved to another (probably non-inline) function.
115
- template<typename T>
116
- CEREAL_RAPIDJSON_FORCEINLINE void Reserve(size_t count = 1) {
117
- // Expand the stack if needed
118
- if (CEREAL_RAPIDJSON_UNLIKELY(static_cast<std::ptrdiff_t>(sizeof(T) * count) > (stackEnd_ - stackTop_)))
119
- Expand<T>(count);
120
- }
121
-
122
- template<typename T>
123
- CEREAL_RAPIDJSON_FORCEINLINE T* Push(size_t count = 1) {
124
- Reserve<T>(count);
125
- return PushUnsafe<T>(count);
126
- }
127
-
128
- template<typename T>
129
- CEREAL_RAPIDJSON_FORCEINLINE T* PushUnsafe(size_t count = 1) {
130
- CEREAL_RAPIDJSON_ASSERT(stackTop_);
131
- CEREAL_RAPIDJSON_ASSERT(static_cast<std::ptrdiff_t>(sizeof(T) * count) <= (stackEnd_ - stackTop_));
132
- T* ret = reinterpret_cast<T*>(stackTop_);
133
- stackTop_ += sizeof(T) * count;
134
- return ret;
135
- }
136
-
137
- template<typename T>
138
- T* Pop(size_t count) {
139
- CEREAL_RAPIDJSON_ASSERT(GetSize() >= count * sizeof(T));
140
- stackTop_ -= count * sizeof(T);
141
- return reinterpret_cast<T*>(stackTop_);
142
- }
143
-
144
- template<typename T>
145
- T* Top() {
146
- CEREAL_RAPIDJSON_ASSERT(GetSize() >= sizeof(T));
147
- return reinterpret_cast<T*>(stackTop_ - sizeof(T));
148
- }
149
-
150
- template<typename T>
151
- const T* Top() const {
152
- CEREAL_RAPIDJSON_ASSERT(GetSize() >= sizeof(T));
153
- return reinterpret_cast<T*>(stackTop_ - sizeof(T));
154
- }
155
-
156
- template<typename T>
157
- T* End() { return reinterpret_cast<T*>(stackTop_); }
158
-
159
- template<typename T>
160
- const T* End() const { return reinterpret_cast<T*>(stackTop_); }
161
-
162
- template<typename T>
163
- T* Bottom() { return reinterpret_cast<T*>(stack_); }
164
-
165
- template<typename T>
166
- const T* Bottom() const { return reinterpret_cast<T*>(stack_); }
167
-
168
- bool HasAllocator() const {
169
- return allocator_ != 0;
170
- }
171
-
172
- Allocator& GetAllocator() {
173
- CEREAL_RAPIDJSON_ASSERT(allocator_);
174
- return *allocator_;
175
- }
176
-
177
- bool Empty() const { return stackTop_ == stack_; }
178
- size_t GetSize() const { return static_cast<size_t>(stackTop_ - stack_); }
179
- size_t GetCapacity() const { return static_cast<size_t>(stackEnd_ - stack_); }
180
-
181
- private:
182
- template<typename T>
183
- void Expand(size_t count) {
184
- // Only expand the capacity if the current stack exists. Otherwise just create a stack with initial capacity.
185
- size_t newCapacity;
186
- if (stack_ == 0) {
187
- if (!allocator_)
188
- ownAllocator_ = allocator_ = CEREAL_RAPIDJSON_NEW(Allocator)();
189
- newCapacity = initialCapacity_;
190
- } else {
191
- newCapacity = GetCapacity();
192
- newCapacity += (newCapacity + 1) / 2;
193
- }
194
- size_t newSize = GetSize() + sizeof(T) * count;
195
- if (newCapacity < newSize)
196
- newCapacity = newSize;
197
-
198
- Resize(newCapacity);
199
- }
200
-
201
- void Resize(size_t newCapacity) {
202
- const size_t size = GetSize(); // Backup the current size
203
- stack_ = static_cast<char*>(allocator_->Realloc(stack_, GetCapacity(), newCapacity));
204
- stackTop_ = stack_ + size;
205
- stackEnd_ = stack_ + newCapacity;
206
- }
207
-
208
- void Destroy() {
209
- Allocator::Free(stack_);
210
- CEREAL_RAPIDJSON_DELETE(ownAllocator_); // Only delete if it is owned by the stack
211
- }
212
-
213
- // Prohibit copy constructor & assignment operator.
214
- Stack(const Stack&);
215
- Stack& operator=(const Stack&);
216
-
217
- Allocator* allocator_;
218
- Allocator* ownAllocator_;
219
- char *stack_;
220
- char *stackTop_;
221
- char *stackEnd_;
222
- size_t initialCapacity_;
223
- };
224
-
225
- } // namespace internal
226
- CEREAL_RAPIDJSON_NAMESPACE_END
227
-
228
- #if defined(__clang__)
229
- CEREAL_RAPIDJSON_DIAG_POP
230
- #endif
231
-
232
- #endif // CEREAL_RAPIDJSON_STACK_H_
@@ -1,69 +0,0 @@
1
- // Tencent is pleased to support the open source community by making RapidJSON available.
2
- //
3
- // Copyright (C) 2015 THL A29 Limited, a Tencent company, and Milo Yip. All rights reserved.
4
- //
5
- // Licensed under the MIT License (the "License"); you may not use this file except
6
- // in compliance with the License. You may obtain a copy of the License at
7
- //
8
- // http://opensource.org/licenses/MIT
9
- //
10
- // Unless required by applicable law or agreed to in writing, software distributed
11
- // under the License is distributed on an "AS IS" BASIS, WITHOUT WARRANTIES OR
12
- // CONDITIONS OF ANY KIND, either express or implied. See the License for the
13
- // specific language governing permissions and limitations under the License.
14
-
15
- #ifndef CEREAL_RAPIDJSON_INTERNAL_STRFUNC_H_
16
- #define CEREAL_RAPIDJSON_INTERNAL_STRFUNC_H_
17
-
18
- #include "../stream.h"
19
- #include <cwchar>
20
-
21
- CEREAL_RAPIDJSON_NAMESPACE_BEGIN
22
- namespace internal {
23
-
24
- //! Custom strlen() which works on different character types.
25
- /*! \tparam Ch Character type (e.g. char, wchar_t, short)
26
- \param s Null-terminated input string.
27
- \return Number of characters in the string.
28
- \note This has the same semantics as strlen(), the return value is not number of Unicode codepoints.
29
- */
30
- template <typename Ch>
31
- inline SizeType StrLen(const Ch* s) {
32
- CEREAL_RAPIDJSON_ASSERT(s != 0);
33
- const Ch* p = s;
34
- while (*p) ++p;
35
- return SizeType(p - s);
36
- }
37
-
38
- template <>
39
- inline SizeType StrLen(const char* s) {
40
- return SizeType(std::strlen(s));
41
- }
42
-
43
- template <>
44
- inline SizeType StrLen(const wchar_t* s) {
45
- return SizeType(std::wcslen(s));
46
- }
47
-
48
- //! Returns number of code points in a encoded string.
49
- template<typename Encoding>
50
- bool CountStringCodePoint(const typename Encoding::Ch* s, SizeType length, SizeType* outCount) {
51
- CEREAL_RAPIDJSON_ASSERT(s != 0);
52
- CEREAL_RAPIDJSON_ASSERT(outCount != 0);
53
- GenericStringStream<Encoding> is(s);
54
- const typename Encoding::Ch* end = s + length;
55
- SizeType count = 0;
56
- while (is.src_ < end) {
57
- unsigned codepoint;
58
- if (!Encoding::Decode(is, &codepoint))
59
- return false;
60
- count++;
61
- }
62
- *outCount = count;
63
- return true;
64
- }
65
-
66
- } // namespace internal
67
- CEREAL_RAPIDJSON_NAMESPACE_END
68
-
69
- #endif // CEREAL_RAPIDJSON_INTERNAL_STRFUNC_H_
@@ -1,290 +0,0 @@
1
- // Tencent is pleased to support the open source community by making RapidJSON available.
2
- //
3
- // Copyright (C) 2015 THL A29 Limited, a Tencent company, and Milo Yip. All rights reserved.
4
- //
5
- // Licensed under the MIT License (the "License"); you may not use this file except
6
- // in compliance with the License. You may obtain a copy of the License at
7
- //
8
- // http://opensource.org/licenses/MIT
9
- //
10
- // Unless required by applicable law or agreed to in writing, software distributed
11
- // under the License is distributed on an "AS IS" BASIS, WITHOUT WARRANTIES OR
12
- // CONDITIONS OF ANY KIND, either express or implied. See the License for the
13
- // specific language governing permissions and limitations under the License.
14
-
15
- #ifndef CEREAL_RAPIDJSON_STRTOD_
16
- #define CEREAL_RAPIDJSON_STRTOD_
17
-
18
- #include "ieee754.h"
19
- #include "biginteger.h"
20
- #include "diyfp.h"
21
- #include "pow10.h"
22
- #include <climits>
23
- #include <limits>
24
-
25
- CEREAL_RAPIDJSON_NAMESPACE_BEGIN
26
- namespace internal {
27
-
28
- inline double FastPath(double significand, int exp) {
29
- if (exp < -308)
30
- return 0.0;
31
- else if (exp >= 0)
32
- return significand * internal::Pow10(exp);
33
- else
34
- return significand / internal::Pow10(-exp);
35
- }
36
-
37
- inline double StrtodNormalPrecision(double d, int p) {
38
- if (p < -308) {
39
- // Prevent expSum < -308, making Pow10(p) = 0
40
- d = FastPath(d, -308);
41
- d = FastPath(d, p + 308);
42
- }
43
- else
44
- d = FastPath(d, p);
45
- return d;
46
- }
47
-
48
- template <typename T>
49
- inline T Min3(T a, T b, T c) {
50
- T m = a;
51
- if (m > b) m = b;
52
- if (m > c) m = c;
53
- return m;
54
- }
55
-
56
- inline int CheckWithinHalfULP(double b, const BigInteger& d, int dExp) {
57
- const Double db(b);
58
- const uint64_t bInt = db.IntegerSignificand();
59
- const int bExp = db.IntegerExponent();
60
- const int hExp = bExp - 1;
61
-
62
- int dS_Exp2 = 0, dS_Exp5 = 0, bS_Exp2 = 0, bS_Exp5 = 0, hS_Exp2 = 0, hS_Exp5 = 0;
63
-
64
- // Adjust for decimal exponent
65
- if (dExp >= 0) {
66
- dS_Exp2 += dExp;
67
- dS_Exp5 += dExp;
68
- }
69
- else {
70
- bS_Exp2 -= dExp;
71
- bS_Exp5 -= dExp;
72
- hS_Exp2 -= dExp;
73
- hS_Exp5 -= dExp;
74
- }
75
-
76
- // Adjust for binary exponent
77
- if (bExp >= 0)
78
- bS_Exp2 += bExp;
79
- else {
80
- dS_Exp2 -= bExp;
81
- hS_Exp2 -= bExp;
82
- }
83
-
84
- // Adjust for half ulp exponent
85
- if (hExp >= 0)
86
- hS_Exp2 += hExp;
87
- else {
88
- dS_Exp2 -= hExp;
89
- bS_Exp2 -= hExp;
90
- }
91
-
92
- // Remove common power of two factor from all three scaled values
93
- int common_Exp2 = Min3(dS_Exp2, bS_Exp2, hS_Exp2);
94
- dS_Exp2 -= common_Exp2;
95
- bS_Exp2 -= common_Exp2;
96
- hS_Exp2 -= common_Exp2;
97
-
98
- BigInteger dS = d;
99
- dS.MultiplyPow5(static_cast<unsigned>(dS_Exp5)) <<= static_cast<unsigned>(dS_Exp2);
100
-
101
- BigInteger bS(bInt);
102
- bS.MultiplyPow5(static_cast<unsigned>(bS_Exp5)) <<= static_cast<unsigned>(bS_Exp2);
103
-
104
- BigInteger hS(1);
105
- hS.MultiplyPow5(static_cast<unsigned>(hS_Exp5)) <<= static_cast<unsigned>(hS_Exp2);
106
-
107
- BigInteger delta(0);
108
- dS.Difference(bS, &delta);
109
-
110
- return delta.Compare(hS);
111
- }
112
-
113
- inline bool StrtodFast(double d, int p, double* result) {
114
- // Use fast path for string-to-double conversion if possible
115
- // see http://www.exploringbinary.com/fast-path-decimal-to-floating-point-conversion/
116
- if (p > 22 && p < 22 + 16) {
117
- // Fast Path Cases In Disguise
118
- d *= internal::Pow10(p - 22);
119
- p = 22;
120
- }
121
-
122
- if (p >= -22 && p <= 22 && d <= 9007199254740991.0) { // 2^53 - 1
123
- *result = FastPath(d, p);
124
- return true;
125
- }
126
- else
127
- return false;
128
- }
129
-
130
- // Compute an approximation and see if it is within 1/2 ULP
131
- inline bool StrtodDiyFp(const char* decimals, int dLen, int dExp, double* result) {
132
- uint64_t significand = 0;
133
- int i = 0; // 2^64 - 1 = 18446744073709551615, 1844674407370955161 = 0x1999999999999999
134
- for (; i < dLen; i++) {
135
- if (significand > CEREAL_RAPIDJSON_UINT64_C2(0x19999999, 0x99999999) ||
136
- (significand == CEREAL_RAPIDJSON_UINT64_C2(0x19999999, 0x99999999) && decimals[i] > '5'))
137
- break;
138
- significand = significand * 10u + static_cast<unsigned>(decimals[i] - '0');
139
- }
140
-
141
- if (i < dLen && decimals[i] >= '5') // Rounding
142
- significand++;
143
-
144
- int remaining = dLen - i;
145
- const int kUlpShift = 3;
146
- const int kUlp = 1 << kUlpShift;
147
- int64_t error = (remaining == 0) ? 0 : kUlp / 2;
148
-
149
- DiyFp v(significand, 0);
150
- v = v.Normalize();
151
- error <<= -v.e;
152
-
153
- dExp += remaining;
154
-
155
- int actualExp;
156
- DiyFp cachedPower = GetCachedPower10(dExp, &actualExp);
157
- if (actualExp != dExp) {
158
- static const DiyFp kPow10[] = {
159
- DiyFp(CEREAL_RAPIDJSON_UINT64_C2(0xa0000000, 0x00000000), -60), // 10^1
160
- DiyFp(CEREAL_RAPIDJSON_UINT64_C2(0xc8000000, 0x00000000), -57), // 10^2
161
- DiyFp(CEREAL_RAPIDJSON_UINT64_C2(0xfa000000, 0x00000000), -54), // 10^3
162
- DiyFp(CEREAL_RAPIDJSON_UINT64_C2(0x9c400000, 0x00000000), -50), // 10^4
163
- DiyFp(CEREAL_RAPIDJSON_UINT64_C2(0xc3500000, 0x00000000), -47), // 10^5
164
- DiyFp(CEREAL_RAPIDJSON_UINT64_C2(0xf4240000, 0x00000000), -44), // 10^6
165
- DiyFp(CEREAL_RAPIDJSON_UINT64_C2(0x98968000, 0x00000000), -40) // 10^7
166
- };
167
- int adjustment = dExp - actualExp;
168
- CEREAL_RAPIDJSON_ASSERT(adjustment >= 1 && adjustment < 8);
169
- v = v * kPow10[adjustment - 1];
170
- if (dLen + adjustment > 19) // has more digits than decimal digits in 64-bit
171
- error += kUlp / 2;
172
- }
173
-
174
- v = v * cachedPower;
175
-
176
- error += kUlp + (error == 0 ? 0 : 1);
177
-
178
- const int oldExp = v.e;
179
- v = v.Normalize();
180
- error <<= oldExp - v.e;
181
-
182
- const int effectiveSignificandSize = Double::EffectiveSignificandSize(64 + v.e);
183
- int precisionSize = 64 - effectiveSignificandSize;
184
- if (precisionSize + kUlpShift >= 64) {
185
- int scaleExp = (precisionSize + kUlpShift) - 63;
186
- v.f >>= scaleExp;
187
- v.e += scaleExp;
188
- error = (error >> scaleExp) + 1 + kUlp;
189
- precisionSize -= scaleExp;
190
- }
191
-
192
- DiyFp rounded(v.f >> precisionSize, v.e + precisionSize);
193
- const uint64_t precisionBits = (v.f & ((uint64_t(1) << precisionSize) - 1)) * kUlp;
194
- const uint64_t halfWay = (uint64_t(1) << (precisionSize - 1)) * kUlp;
195
- if (precisionBits >= halfWay + static_cast<unsigned>(error)) {
196
- rounded.f++;
197
- if (rounded.f & (DiyFp::kDpHiddenBit << 1)) { // rounding overflows mantissa (issue #340)
198
- rounded.f >>= 1;
199
- rounded.e++;
200
- }
201
- }
202
-
203
- *result = rounded.ToDouble();
204
-
205
- return halfWay - static_cast<unsigned>(error) >= precisionBits || precisionBits >= halfWay + static_cast<unsigned>(error);
206
- }
207
-
208
- inline double StrtodBigInteger(double approx, const char* decimals, int dLen, int dExp) {
209
- CEREAL_RAPIDJSON_ASSERT(dLen >= 0);
210
- const BigInteger dInt(decimals, static_cast<unsigned>(dLen));
211
- Double a(approx);
212
- int cmp = CheckWithinHalfULP(a.Value(), dInt, dExp);
213
- if (cmp < 0)
214
- return a.Value(); // within half ULP
215
- else if (cmp == 0) {
216
- // Round towards even
217
- if (a.Significand() & 1)
218
- return a.NextPositiveDouble();
219
- else
220
- return a.Value();
221
- }
222
- else // adjustment
223
- return a.NextPositiveDouble();
224
- }
225
-
226
- inline double StrtodFullPrecision(double d, int p, const char* decimals, size_t length, size_t decimalPosition, int exp) {
227
- CEREAL_RAPIDJSON_ASSERT(d >= 0.0);
228
- CEREAL_RAPIDJSON_ASSERT(length >= 1);
229
-
230
- double result = 0.0;
231
- if (StrtodFast(d, p, &result))
232
- return result;
233
-
234
- CEREAL_RAPIDJSON_ASSERT(length <= INT_MAX);
235
- int dLen = static_cast<int>(length);
236
-
237
- CEREAL_RAPIDJSON_ASSERT(length >= decimalPosition);
238
- CEREAL_RAPIDJSON_ASSERT(length - decimalPosition <= INT_MAX);
239
- int dExpAdjust = static_cast<int>(length - decimalPosition);
240
-
241
- CEREAL_RAPIDJSON_ASSERT(exp >= INT_MIN + dExpAdjust);
242
- int dExp = exp - dExpAdjust;
243
-
244
- // Make sure length+dExp does not overflow
245
- CEREAL_RAPIDJSON_ASSERT(dExp <= INT_MAX - dLen);
246
-
247
- // Trim leading zeros
248
- while (dLen > 0 && *decimals == '0') {
249
- dLen--;
250
- decimals++;
251
- }
252
-
253
- // Trim trailing zeros
254
- while (dLen > 0 && decimals[dLen - 1] == '0') {
255
- dLen--;
256
- dExp++;
257
- }
258
-
259
- if (dLen == 0) { // Buffer only contains zeros.
260
- return 0.0;
261
- }
262
-
263
- // Trim right-most digits
264
- const int kMaxDecimalDigit = 767 + 1;
265
- if (dLen > kMaxDecimalDigit) {
266
- dExp += dLen - kMaxDecimalDigit;
267
- dLen = kMaxDecimalDigit;
268
- }
269
-
270
- // If too small, underflow to zero.
271
- // Any x <= 10^-324 is interpreted as zero.
272
- if (dLen + dExp <= -324)
273
- return 0.0;
274
-
275
- // If too large, overflow to infinity.
276
- // Any x >= 10^309 is interpreted as +infinity.
277
- if (dLen + dExp > 309)
278
- return std::numeric_limits<double>::infinity();
279
-
280
- if (StrtodDiyFp(decimals, dLen, dExp, &result))
281
- return result;
282
-
283
- // Use approximation from StrtodDiyFp and make adjustment with BigInteger comparison
284
- return StrtodBigInteger(result, decimals, dLen, dExp);
285
- }
286
-
287
- } // namespace internal
288
- CEREAL_RAPIDJSON_NAMESPACE_END
289
-
290
- #endif // CEREAL_RAPIDJSON_STRTOD_
@@ -1,46 +0,0 @@
1
- // Tencent is pleased to support the open source community by making RapidJSON available.
2
- //
3
- // Copyright (C) 2015 THL A29 Limited, a Tencent company, and Milo Yip. All rights reserved.
4
- //
5
- // Licensed under the MIT License (the "License"); you may not use this file except
6
- // in compliance with the License. You may obtain a copy of the License at
7
- //
8
- // http://opensource.org/licenses/MIT
9
- //
10
- // Unless required by applicable law or agreed to in writing, software distributed
11
- // under the License is distributed on an "AS IS" BASIS, WITHOUT WARRANTIES OR
12
- // CONDITIONS OF ANY KIND, either express or implied. See the License for the
13
- // specific language governing permissions and limitations under the License.
14
-
15
- #ifndef CEREAL_RAPIDJSON_INTERNAL_SWAP_H_
16
- #define CEREAL_RAPIDJSON_INTERNAL_SWAP_H_
17
-
18
- #include "../rapidjson.h"
19
-
20
- #if defined(__clang__)
21
- CEREAL_RAPIDJSON_DIAG_PUSH
22
- CEREAL_RAPIDJSON_DIAG_OFF(c++98-compat)
23
- #endif
24
-
25
- CEREAL_RAPIDJSON_NAMESPACE_BEGIN
26
- namespace internal {
27
-
28
- //! Custom swap() to avoid dependency on C++ <algorithm> header
29
- /*! \tparam T Type of the arguments to swap, should be instantiated with primitive C++ types only.
30
- \note This has the same semantics as std::swap().
31
- */
32
- template <typename T>
33
- inline void Swap(T& a, T& b) CEREAL_RAPIDJSON_NOEXCEPT {
34
- T tmp = a;
35
- a = b;
36
- b = tmp;
37
- }
38
-
39
- } // namespace internal
40
- CEREAL_RAPIDJSON_NAMESPACE_END
41
-
42
- #if defined(__clang__)
43
- CEREAL_RAPIDJSON_DIAG_POP
44
- #endif
45
-
46
- #endif // CEREAL_RAPIDJSON_INTERNAL_SWAP_H_