/root/bitcoin/src/common/bloom.cpp
Line | Count | Source |
1 | | // Copyright (c) 2012-present The Bitcoin Core developers |
2 | | // Distributed under the MIT software license, see the accompanying |
3 | | // file COPYING or http://www.opensource.org/licenses/mit-license.php. |
4 | | |
5 | | #include <common/bloom.h> |
6 | | |
7 | | #include <hash.h> |
8 | | #include <primitives/transaction.h> |
9 | | #include <random.h> |
10 | | #include <script/script.h> |
11 | | #include <script/solver.h> |
12 | | #include <span.h> |
13 | | #include <streams.h> |
14 | | #include <util/fastrange.h> |
15 | | #include <util/overflow.h> |
16 | | |
17 | | #include <algorithm> |
18 | | #include <cmath> |
19 | | #include <cstdlib> |
20 | | #include <limits> |
21 | | #include <vector> |
22 | | |
23 | | static constexpr double LN2SQUARED = 0.4804530139182014246671025263266649717305529515945455; |
24 | | static constexpr double LN2 = 0.6931471805599453094172321214581765680755001343602552; |
25 | | |
26 | | CBloomFilter::CBloomFilter(const unsigned int nElements, const double nFPRate, const unsigned int nTweakIn, unsigned char nFlagsIn) : |
27 | | /** |
28 | | * The ideal size for a bloom filter with a given number of elements and false positive rate is: |
29 | | * - nElements * log(fp rate) / ln(2)^2 |
30 | | * We ignore filter parameters which will create a bloom filter larger than the protocol limits |
31 | | */ |
32 | 0 | vData(std::min((unsigned int)(-1 / LN2SQUARED * nElements * log(nFPRate)), MAX_BLOOM_FILTER_SIZE * 8) / 8), |
33 | | /** |
34 | | * The ideal number of hash functions is filter size * ln(2) / number of elements |
35 | | * Again, we ignore filter parameters which will create a bloom filter with more hash functions than the protocol limits |
36 | | * See https://en.wikipedia.org/wiki/Bloom_filter for an explanation of these formulas |
37 | | */ |
38 | 0 | nHashFuncs(std::min((unsigned int)(vData.size() * 8 / nElements * LN2), MAX_HASH_FUNCS)), |
39 | 0 | nTweak(nTweakIn), |
40 | 0 | nFlags(nFlagsIn) |
41 | 0 | { |
42 | 0 | } |
43 | | |
44 | | inline unsigned int CBloomFilter::Hash(unsigned int nHashNum, std::span<const unsigned char> vDataToHash) const |
45 | 0 | { |
46 | | // 0xFBA4C795 chosen as it guarantees a reasonable bit difference between nHashNum values. |
47 | 0 | return MurmurHash3(nHashNum * 0xFBA4C795 + nTweak, vDataToHash) % (vData.size() * 8); |
48 | 0 | } |
49 | | |
50 | | void CBloomFilter::insert(std::span<const unsigned char> vKey) |
51 | 0 | { |
52 | 0 | if (vData.empty()) // Avoid divide-by-zero (CVE-2013-5700) |
53 | 0 | return; |
54 | 0 | for (unsigned int i = 0; i < nHashFuncs; i++) |
55 | 0 | { |
56 | 0 | unsigned int nIndex = Hash(i, vKey); |
57 | | // Sets bit nIndex of vData |
58 | 0 | vData[nIndex >> 3] |= (1 << (7 & nIndex)); |
59 | 0 | } |
60 | 0 | } |
61 | | |
62 | | void CBloomFilter::insert(const COutPoint& outpoint) |
63 | 0 | { |
64 | 0 | DataStream stream{}; |
65 | 0 | stream << outpoint; |
66 | 0 | insert(MakeUCharSpan(stream)); |
67 | 0 | } |
68 | | |
69 | | bool CBloomFilter::contains(std::span<const unsigned char> vKey) const |
70 | 0 | { |
71 | 0 | if (vData.empty()) // Avoid divide-by-zero (CVE-2013-5700) |
72 | 0 | return true; |
73 | 0 | for (unsigned int i = 0; i < nHashFuncs; i++) |
74 | 0 | { |
75 | 0 | unsigned int nIndex = Hash(i, vKey); |
76 | | // Checks bit nIndex of vData |
77 | 0 | if (!(vData[nIndex >> 3] & (1 << (7 & nIndex)))) |
78 | 0 | return false; |
79 | 0 | } |
80 | 0 | return true; |
81 | 0 | } |
82 | | |
83 | | bool CBloomFilter::contains(const COutPoint& outpoint) const |
84 | 0 | { |
85 | 0 | DataStream stream{}; |
86 | 0 | stream << outpoint; |
87 | 0 | return contains(MakeUCharSpan(stream)); |
88 | 0 | } |
89 | | |
90 | | bool CBloomFilter::IsWithinSizeConstraints() const |
91 | 0 | { |
92 | 0 | return vData.size() <= MAX_BLOOM_FILTER_SIZE && nHashFuncs <= MAX_HASH_FUNCS; |
93 | 0 | } |
94 | | |
95 | | bool CBloomFilter::IsRelevantAndUpdate(const CTransaction& tx) |
96 | 0 | { |
97 | 0 | bool fFound = false; |
98 | | // Match if the filter contains the hash of tx |
99 | | // for finding tx when they appear in a block |
100 | 0 | if (vData.empty()) // zero-size = "match-all" filter |
101 | 0 | return true; |
102 | 0 | const Txid& hash = tx.GetHash(); |
103 | 0 | if (contains(hash.ToUint256())) |
104 | 0 | fFound = true; |
105 | |
|
106 | 0 | for (unsigned int i = 0; i < tx.vout.size(); i++) |
107 | 0 | { |
108 | 0 | const CTxOut& txout = tx.vout[i]; |
109 | | // Match if the filter contains any arbitrary script data element in any scriptPubKey in tx |
110 | | // If this matches, also add the specific output that was matched. |
111 | | // This means clients don't have to update the filter themselves when a new relevant tx |
112 | | // is discovered in order to find spending transactions, which avoids round-tripping and race conditions. |
113 | 0 | CScript::const_iterator pc = txout.scriptPubKey.begin(); |
114 | 0 | std::vector<unsigned char> data; |
115 | 0 | while (pc < txout.scriptPubKey.end()) |
116 | 0 | { |
117 | 0 | opcodetype opcode; |
118 | 0 | if (!txout.scriptPubKey.GetOp(pc, opcode, data)) |
119 | 0 | break; |
120 | 0 | if (data.size() != 0 && contains(data)) |
121 | 0 | { |
122 | 0 | fFound = true; |
123 | 0 | if ((nFlags & BLOOM_UPDATE_MASK) == BLOOM_UPDATE_ALL) |
124 | 0 | insert(COutPoint(hash, i)); |
125 | 0 | else if ((nFlags & BLOOM_UPDATE_MASK) == BLOOM_UPDATE_P2PUBKEY_ONLY) |
126 | 0 | { |
127 | 0 | std::vector<std::vector<unsigned char> > vSolutions; |
128 | 0 | TxoutType type = Solver(txout.scriptPubKey, vSolutions); |
129 | 0 | if (type == TxoutType::PUBKEY || type == TxoutType::MULTISIG) { |
130 | 0 | insert(COutPoint(hash, i)); |
131 | 0 | } |
132 | 0 | } |
133 | 0 | break; |
134 | 0 | } |
135 | 0 | } |
136 | 0 | } |
137 | |
|
138 | 0 | if (fFound) |
139 | 0 | return true; |
140 | | |
141 | 0 | for (const CTxIn& txin : tx.vin) |
142 | 0 | { |
143 | | // Match if the filter contains an outpoint tx spends |
144 | 0 | if (contains(txin.prevout)) |
145 | 0 | return true; |
146 | | |
147 | | // Match if the filter contains any arbitrary script data element in any scriptSig in tx |
148 | 0 | CScript::const_iterator pc = txin.scriptSig.begin(); |
149 | 0 | std::vector<unsigned char> data; |
150 | 0 | while (pc < txin.scriptSig.end()) |
151 | 0 | { |
152 | 0 | opcodetype opcode; |
153 | 0 | if (!txin.scriptSig.GetOp(pc, opcode, data)) |
154 | 0 | break; |
155 | 0 | if (data.size() != 0 && contains(data)) |
156 | 0 | return true; |
157 | 0 | } |
158 | 0 | } |
159 | | |
160 | 0 | return false; |
161 | 0 | } |
162 | | |
163 | | CRollingBloomFilter::CRollingBloomFilter(const unsigned int nElements, const double fpRate) |
164 | 0 | { |
165 | 0 | double logFpRate = log(fpRate); |
166 | | /* The optimal number of hash functions is log(fpRate) / log(0.5), but |
167 | | * restrict it to the range 1-50. */ |
168 | 0 | nHashFuncs = std::max(1, std::min((int)round(logFpRate / log(0.5)), 50)); |
169 | | /* In this rolling bloom filter, we'll store between 2 and 3 generations of nElements / 2 entries. */ |
170 | 0 | nEntriesPerGeneration = CeilDiv(nElements, 2u); |
171 | 0 | uint32_t nMaxElements = nEntriesPerGeneration * 3; |
172 | | /* The maximum fpRate = pow(1.0 - exp(-nHashFuncs * nMaxElements / nFilterBits), nHashFuncs) |
173 | | * => pow(fpRate, 1.0 / nHashFuncs) = 1.0 - exp(-nHashFuncs * nMaxElements / nFilterBits) |
174 | | * => 1.0 - pow(fpRate, 1.0 / nHashFuncs) = exp(-nHashFuncs * nMaxElements / nFilterBits) |
175 | | * => log(1.0 - pow(fpRate, 1.0 / nHashFuncs)) = -nHashFuncs * nMaxElements / nFilterBits |
176 | | * => nFilterBits = -nHashFuncs * nMaxElements / log(1.0 - pow(fpRate, 1.0 / nHashFuncs)) |
177 | | * => nFilterBits = -nHashFuncs * nMaxElements / log(1.0 - exp(logFpRate / nHashFuncs)) |
178 | | */ |
179 | 0 | uint32_t nFilterBits = (uint32_t)ceil(-1.0 * nHashFuncs * nMaxElements / log(1.0 - exp(logFpRate / nHashFuncs))); |
180 | 0 | data.clear(); |
181 | | /* For each data element we need to store 2 bits. If both bits are 0, the |
182 | | * bit is treated as unset. If the bits are (01), (10), or (11), the bit is |
183 | | * treated as set in generation 1, 2, or 3 respectively. |
184 | | * These bits are stored in separate integers: position P corresponds to bit |
185 | | * (P & 63) of the integers data[(P >> 6) * 2] and data[(P >> 6) * 2 + 1]. */ |
186 | 0 | data.resize(CeilDiv(nFilterBits, 64u) << 1); |
187 | 0 | reset(); |
188 | 0 | } |
189 | | |
190 | | /* Similar to CBloomFilter::Hash */ |
191 | | static inline uint32_t RollingBloomHash(unsigned int nHashNum, uint32_t nTweak, std::span<const unsigned char> vDataToHash) |
192 | 0 | { |
193 | 0 | return MurmurHash3(nHashNum * 0xFBA4C795 + nTweak, vDataToHash); |
194 | 0 | } |
195 | | |
196 | | void CRollingBloomFilter::insert(std::span<const unsigned char> vKey) |
197 | 0 | { |
198 | 0 | if (nEntriesThisGeneration == nEntriesPerGeneration) { |
199 | 0 | nEntriesThisGeneration = 0; |
200 | 0 | nGeneration++; |
201 | 0 | if (nGeneration == 4) { |
202 | 0 | nGeneration = 1; |
203 | 0 | } |
204 | 0 | uint64_t nGenerationMask1 = 0 - (uint64_t)(nGeneration & 1); |
205 | 0 | uint64_t nGenerationMask2 = 0 - (uint64_t)(nGeneration >> 1); |
206 | | /* Wipe old entries that used this generation number. */ |
207 | 0 | for (uint32_t p = 0; p < data.size(); p += 2) { |
208 | 0 | uint64_t p1 = data[p], p2 = data[p + 1]; |
209 | 0 | uint64_t mask = (p1 ^ nGenerationMask1) | (p2 ^ nGenerationMask2); |
210 | 0 | data[p] = p1 & mask; |
211 | 0 | data[p + 1] = p2 & mask; |
212 | 0 | } |
213 | 0 | } |
214 | 0 | nEntriesThisGeneration++; |
215 | |
|
216 | 0 | for (int n = 0; n < nHashFuncs; n++) { |
217 | 0 | uint32_t h = RollingBloomHash(n, nTweak, vKey); |
218 | 0 | int bit = h & 0x3F; |
219 | | /* FastMod works with the upper bits of h, so it is safe to ignore that the lower bits of h are already used for bit. */ |
220 | 0 | uint32_t pos = FastRange32(h, data.size()); |
221 | | /* The lowest bit of pos is ignored, and set to zero for the first bit, and to one for the second. */ |
222 | 0 | data[pos & ~1U] = (data[pos & ~1U] & ~(uint64_t{1} << bit)) | (uint64_t(nGeneration & 1)) << bit; |
223 | 0 | data[pos | 1] = (data[pos | 1] & ~(uint64_t{1} << bit)) | (uint64_t(nGeneration >> 1)) << bit; |
224 | 0 | } |
225 | 0 | } |
226 | | |
227 | | bool CRollingBloomFilter::contains(std::span<const unsigned char> vKey) const |
228 | 0 | { |
229 | 0 | for (int n = 0; n < nHashFuncs; n++) { |
230 | 0 | uint32_t h = RollingBloomHash(n, nTweak, vKey); |
231 | 0 | int bit = h & 0x3F; |
232 | 0 | uint32_t pos = FastRange32(h, data.size()); |
233 | | /* If the relevant bit is not set in either data[pos & ~1] or data[pos | 1], the filter does not contain vKey */ |
234 | 0 | if (!(((data[pos & ~1U] | data[pos | 1]) >> bit) & 1)) { |
235 | 0 | return false; |
236 | 0 | } |
237 | 0 | } |
238 | 0 | return true; |
239 | 0 | } |
240 | | |
241 | | void CRollingBloomFilter::reset() |
242 | 0 | { |
243 | 0 | nTweak = FastRandomContext().rand<unsigned int>(); |
244 | 0 | nEntriesThisGeneration = 0; |
245 | 0 | nGeneration = 1; |
246 | 0 | std::fill(data.begin(), data.end(), 0); |
247 | 0 | } |