• Home
  • Features
  • Pricing
  • Docs
  • Announcements
  • Sign In

randombit / botan / 26735862306

31 May 2026 11:43PM UTC coverage: 89.37%. Remained the same
26735862306

push

github

web-flow
Merge pull request #5631 from randombit/jack/overflow-checks

Systematically eliminate any possible integer overflows

110295 of 123414 relevant lines covered (89.37%)

11075877.11 hits per line

Source File
Press 'n' to go to next uncovered line, 'b' for previous

98.37
/src/lib/utils/bitvector/bitvector.h
1
/*
2
 * An abstraction for an arbitrarily large bitvector that can
3
 * optionally use the secure_allocator. All bitwise accesses and all
4
 * constructors are implemented in constant time. Otherwise, only methods
5
 * with the "ct_" pre-fix run in constant time.
6
 *
7
 * (C) 2023-2024 Jack Lloyd
8
 * (C) 2023-2024 René Meusel, Rohde & Schwarz Cybersecurity
9
 *
10
 * Botan is released under the Simplified BSD License (see license.txt)
11
 */
12

13
#ifndef BOTAN_BIT_VECTOR_H_
14
#define BOTAN_BIT_VECTOR_H_
15

16
#include <botan/concepts.h>
17
#include <botan/exceptn.h>
18
#include <botan/mem_ops.h>
19
#include <botan/secmem.h>
20
#include <botan/strong_type.h>
21
#include <botan/internal/bit_ops.h>
22
#include <botan/internal/ct_utils.h>
23
#include <botan/internal/int_utils.h>
24
#include <botan/internal/loadstor.h>
25
#include <botan/internal/stl_util.h>
26

27
#include <memory>
28
#include <optional>
29
#include <span>
30
#include <sstream>
31
#include <string>
32
#include <utility>
33
#include <vector>
34

35
namespace Botan {
36

37
template <template <typename> typename AllocatorT>
38
class bitvector_base;
39

40
template <typename T>
41
struct is_bitvector : std::false_type {};
42

43
template <template <typename> typename T>
44
struct is_bitvector<bitvector_base<T>> : std::true_type {};
45

46
template <typename T>
47
constexpr static bool is_bitvector_v = is_bitvector<T>::value;
48

49
template <typename T>
50
concept bitvectorish = is_bitvector_v<strong_type_wrapped_type<T>>;
51

52
namespace detail {
53

54
template <typename T0, typename... Ts>
55
struct first_type {
56
      using type = T0;
57
};
58

59
// get the first type from a parameter pack
60
// TODO: C++26 will bring Parameter Pack indexing:
61
//       using first_t = Ts...[0];
62
template <typename... Ts>
63
   requires(sizeof...(Ts) > 0)
64
using first_t = typename first_type<Ts...>::type;
65

66
// get the first object from a parameter pack
67
// TODO: C++26 will bring Parameter Pack indexing:
68
//       auto first = s...[0];
69
template <typename T0, typename... Ts>
70
constexpr static first_t<T0, Ts...> first(T0&& t, Ts&&... /*rest*/) {
71
   return std::forward<T0>(t);
72
}
73

74
template <typename OutT, typename>
75
using as = OutT;
76

77
template <typename FnT, std::unsigned_integral BlockT, typename... ParamTs>
78
using blockwise_processing_callback_return_type = std::invoke_result_t<FnT, as<BlockT, ParamTs>...>;
79

80
template <typename FnT, typename BlockT, typename... ParamTs>
81
concept is_blockwise_processing_callback_return_type =
82
   std::unsigned_integral<BlockT> &&
83
   (std::same_as<BlockT, blockwise_processing_callback_return_type<FnT, BlockT, ParamTs...>> ||
84
    std::same_as<bool, blockwise_processing_callback_return_type<FnT, BlockT, ParamTs...>> ||
85
    std::same_as<void, blockwise_processing_callback_return_type<FnT, BlockT, ParamTs...>>);
86

87
template <typename FnT, typename... ParamTs>
88
concept blockwise_processing_callback_without_mask =
89
   is_blockwise_processing_callback_return_type<FnT, uint8_t, ParamTs...> &&
90
   is_blockwise_processing_callback_return_type<FnT, uint16_t, ParamTs...> &&
91
   is_blockwise_processing_callback_return_type<FnT, uint32_t, ParamTs...> &&
92
   is_blockwise_processing_callback_return_type<FnT, uint64_t, ParamTs...>;
93

94
template <typename FnT, typename... ParamTs>
95
concept blockwise_processing_callback_with_mask =
96
   is_blockwise_processing_callback_return_type<FnT, uint8_t, ParamTs..., uint8_t /* mask */> &&
97
   is_blockwise_processing_callback_return_type<FnT, uint16_t, ParamTs..., uint16_t /* mask */> &&
98
   is_blockwise_processing_callback_return_type<FnT, uint32_t, ParamTs..., uint32_t /* mask */> &&
99
   is_blockwise_processing_callback_return_type<FnT, uint64_t, ParamTs..., uint64_t /* mask */>;
100

101
/**
102
 * Defines the callback constraints for the BitRangeOperator. For further
103
 * details, see bitvector_base::range_operation().
104
 */
105
template <typename FnT, typename... ParamTs>
106
concept blockwise_processing_callback = blockwise_processing_callback_with_mask<FnT, ParamTs...> ||
107
                                        blockwise_processing_callback_without_mask<FnT, ParamTs...>;
108

109
template <typename FnT, typename... ParamTs>
110
concept manipulating_blockwise_processing_callback =
111
   (blockwise_processing_callback_without_mask<FnT, ParamTs...> &&
112
    std::same_as<uint32_t, blockwise_processing_callback_return_type<FnT, uint32_t, ParamTs...>>) ||
113
   (blockwise_processing_callback_with_mask<FnT, ParamTs...> &&
114
    std::same_as<uint32_t, blockwise_processing_callback_return_type<FnT, uint32_t, ParamTs..., first_t<ParamTs...>>>);
115

116
template <typename FnT, typename... ParamTs>
117
concept predicate_blockwise_processing_callback =
118
   (blockwise_processing_callback_without_mask<FnT, ParamTs...> &&
119
    std::same_as<bool, blockwise_processing_callback_return_type<FnT, uint32_t, ParamTs...>>) ||
120
   (blockwise_processing_callback_with_mask<FnT, ParamTs...> &&
121
    std::same_as<bool, blockwise_processing_callback_return_type<FnT, uint32_t, ParamTs..., first_t<ParamTs...>>>);
122

123
template <typename T>
124
class bitvector_iterator {
125
   private:
126
      using size_type = typename T::size_type;
127

128
   public:
129
      using difference_type = std::make_signed_t<size_type>;
130
      using value_type = std::remove_const_t<decltype(std::declval<T>().at(0))>;
131
      using pointer = value_type*;
132
      using reference = value_type&;
133

134
      // TODO: technically, this could be a random access iterator
135
      using iterator_category = std::bidirectional_iterator_tag;
136

137
   public:
138
      bitvector_iterator() = default;
139
      ~bitvector_iterator() = default;
140

141
      bitvector_iterator(T* bitvector, size_t offset) : m_bitvector(bitvector) { update(offset); }
1,589✔
142

143
      bitvector_iterator(const bitvector_iterator& other) noexcept : m_bitvector(other.m_bitvector) {
2,002✔
144
         update(other.m_offset);
2,004✔
145
      }
146

147
      bitvector_iterator(bitvector_iterator&& other) noexcept : m_bitvector(other.m_bitvector) {
148
         update(other.m_offset);
149
      }
150

151
      bitvector_iterator& operator=(const bitvector_iterator& other) noexcept {
152
         if(this != &other) {
153
            m_bitvector = other.m_bitvector;
154
            update(other.m_offset);
155
         }
156
         return *this;
157
      }
158

159
      bitvector_iterator& operator=(bitvector_iterator&& other) noexcept {
160
         m_bitvector = other.m_bitvector;
161
         update(other.m_offset);
162
         return *this;
163
      }
164

165
      bitvector_iterator& operator++() noexcept {
252,192✔
166
         update(signed_offset() + 1);
504,411✔
167
         return *this;
168
      }
169

170
      bitvector_iterator operator++(int) noexcept {
6✔
171
         auto copy = *this;
6✔
172
         update(signed_offset() + 1);
12✔
173
         return copy;
6✔
174
      }
175

176
      bitvector_iterator& operator--() noexcept {
6✔
177
         update(signed_offset() - 1);
12✔
178
         return *this;
179
      }
180

181
      bitvector_iterator operator--(int) noexcept {
182
         auto copy = *this;
183
         update(signed_offset() - 1);
184
         return copy;
185
      }
186

187
      std::partial_ordering operator<=>(const bitvector_iterator& other) const noexcept {
188
         if(m_bitvector == other.m_bitvector) {
189
            return m_offset <=> other.m_offset;
190
         } else {
191
            return std::partial_ordering::unordered;
192
         }
193
      }
194

195
      bool operator==(const bitvector_iterator& other) const noexcept {
253,233✔
196
         return m_bitvector == other.m_bitvector && m_offset == other.m_offset;
252,705✔
197
      }
198

199
      reference operator*() const { return m_bitref.value(); }
2,347✔
200

201
      pointer operator->() const { return &(m_bitref.value()); }
256✔
202

203
   private:
204
      void update(size_type new_offset) {
524✔
205
         m_offset = new_offset;
1,029✔
206
         if(m_offset < m_bitvector->size()) {
1,506✔
207
            m_bitref.emplace((*m_bitvector)[m_offset]);
254,213✔
208
         } else {
209
            // end() iterator
210
            m_bitref.reset();
251,752✔
211
         }
212
      }
213

214
      difference_type signed_offset() const { return static_cast<difference_type>(m_offset); }
6✔
215

216
   private:
217
      T* m_bitvector;
218
      size_type m_offset;
219
      mutable std::optional<value_type> m_bitref;
220
};
221

222
}  // namespace detail
223

224
/**
225
 * An arbitrarily large bitvector with typical bit manipulation and convenient
226
 * bitwise access methods. Don't use `bitvector_base` directly, but the type
227
 * aliases::
228
 *
229
 *    * bitvector         - with a standard allocator
230
 *    * secure_bitvector  - with a secure allocator that auto-scrubs the memory
231
 */
232
template <template <typename> typename AllocatorT>
233
class bitvector_base final {
1,034,991✔
234
   public:
235
      using block_type = uint8_t;
236
      using size_type = size_t;
237
      using allocator_type = AllocatorT<block_type>;
238
      using value_type = block_type;
239
      using iterator = detail::bitvector_iterator<bitvector_base<AllocatorT>>;
240
      using const_iterator = detail::bitvector_iterator<const bitvector_base<AllocatorT>>;
241

242
      static constexpr size_type block_size_bytes = sizeof(block_type);
243
      static constexpr size_type block_size_bits = block_size_bytes * 8;
244
      static constexpr bool uses_secure_allocator = std::is_same_v<allocator_type, secure_allocator<block_type>>;
245

246
   private:
247
      template <template <typename> typename FriendAllocatorT>
248
      friend class bitvector_base;
249

250
      static constexpr block_type one = block_type(1);
251

252
      static constexpr size_type block_offset_shift = size_type(3) + ceil_log2(block_size_bytes);
253
      static constexpr size_type block_index_mask = (one << block_offset_shift) - 1;
254

255
      static constexpr size_type block_index(size_type pos) { return pos >> block_offset_shift; }
21,179,056✔
256

257
      static constexpr size_type block_offset(size_type pos) { return pos & block_index_mask; }
1,336,813,918✔
258

259
   private:
260
      /**
261
       * Internal helper to wrap a single bit in the bitvector and provide
262
       * certain convenience access methods.
263
       */
264
      template <typename BlockT>
265
         requires std::same_as<block_type, std::remove_cv_t<BlockT>>
266
      class bitref_base {
267
         private:
268
            friend class bitvector_base<AllocatorT>;
269

270
            constexpr bitref_base(std::span<BlockT> blocks, size_type pos) noexcept :
1,781,423,254✔
271
                  m_block(blocks[block_index(pos)]), m_mask(one << block_offset(pos)) {}
1,781,423,254✔
272

273
         public:
274
            bitref_base() = delete;
275
            bitref_base(const bitref_base&) noexcept = default;
276
            bitref_base(bitref_base&&) noexcept = default;
277
            bitref_base& operator=(const bitref_base&) = delete;
278
            bitref_base& operator=(bitref_base&&) = delete;
279

280
            ~bitref_base() = default;
281

282
         public:
283
            // NOLINTNEXTLINE(*-explicit-conversions)
284
            constexpr operator bool() const noexcept { return is_set(); }
13,546,893✔
285

286
            constexpr bool is_set() const noexcept { return (m_block & m_mask) > 0; }
13,546,921✔
287

288
            template <std::unsigned_integral T>
289
            constexpr T as() const noexcept {
312,633✔
290
               return static_cast<T>(is_set());
119,250✔
291
            }
292

293
            constexpr CT::Choice as_choice() const noexcept {
443,004,200✔
294
               return CT::Choice::from_int(static_cast<BlockT>(m_block & m_mask));
443,004,200✔
295
            }
296

297
         protected:
298
            BlockT& m_block;  // NOLINT(*-non-private-member-variable*)
299
            BlockT m_mask;    // NOLINT(*-non-private-member-variable*)
300
      };
301

302
   public:
303
      /**
304
       * Wraps a constant reference into the bitvector. Bit can be accessed
305
       * but not modified.
306
       */
307
      template <typename BlockT>
308
      class bitref final : public bitref_base<BlockT> {
309
         public:
310
            using bitref_base<BlockT>::bitref_base;
12,641,743✔
311
      };
312

313
      /**
314
       * Wraps a modifiable reference into the bitvector. Bit may be accessed
315
       * and modified (e.g. flipped or XOR'ed).
316
       *
317
       * Constant-time operations are used for the bit manipulations. The
318
       * location of the bit in the bit vector may be leaked, though.
319
       */
320
      template <typename BlockT>
321
         requires(!std::is_const_v<BlockT>)
322
      class bitref<BlockT> : public bitref_base<BlockT> {
323
         public:
324
            using bitref_base<BlockT>::bitref_base;
1,324,445,183✔
325

326
            ~bitref() = default;
327
            bitref(const bitref&) noexcept = default;
328
            bitref(bitref&&) noexcept = default;
329

330
            constexpr bitref& set() noexcept {
354✔
331
               this->m_block |= this->m_mask;
8✔
332
               return *this;
250✔
333
            }
334

335
            constexpr bitref& unset() noexcept {
6✔
336
               this->m_block &= ~this->m_mask;
1✔
337
               return *this;
338
            }
339

340
            constexpr bitref& flip() noexcept {
36✔
341
               this->m_block ^= this->m_mask;
7✔
342
               return *this;
343
            }
344

345
            // NOLINTBEGIN
346

347
            constexpr bitref& operator=(bool bit) noexcept {
1,324,228,462✔
348
               this->m_block =
1,324,228,462✔
349
                  CT::Mask<BlockT>::expand(bit).select(this->m_mask | this->m_block, this->m_block & ~this->m_mask);
1,324,228,462✔
350
               return *this;
1,324,228,462✔
351
            }
352

353
            constexpr bitref& operator=(const bitref& bit) noexcept { return *this = bit.is_set(); }
354

355
            constexpr bitref& operator=(bitref&& bit) noexcept { return *this = bit.is_set(); }
356

357
            // NOLINTEND
358

359
            constexpr bitref& operator&=(bool other) noexcept {
4✔
360
               this->m_block &= ~CT::Mask<BlockT>::expand(other).if_not_set_return(this->m_mask);
4✔
361
               return *this;
4✔
362
            }
363

364
            constexpr bitref& operator|=(bool other) noexcept {
1✔
365
               this->m_block |= CT::Mask<BlockT>::expand(other).if_set_return(this->m_mask);
1✔
366
               return *this;
367
            }
368

369
            constexpr bitref& operator^=(bool other) noexcept {
79,045✔
370
               this->m_block ^= CT::Mask<BlockT>::expand(other).if_set_return(this->m_mask);
79,045✔
371
               return *this;
372
            }
373
      };
374

375
   public:
376
      bitvector_base() : m_bits(0) {}
78✔
377

378
      explicit bitvector_base(size_type bits) : m_bits(bits), m_blocks(ceil_toblocks(bits)) {}
386,171✔
379

380
      /**
381
       * Initialize the bitvector from a byte-array. Bits are taken byte-wise
382
       * from least significant to most significant. Example::
383
       *
384
       *    bitvector[0] -> LSB(Byte[0])
385
       *    bitvector[1] -> LSB+1(Byte[0])
386
       *    ...
387
       *    bitvector[8] -> LSB(Byte[1])
388
       *
389
       * @param bytes The byte vector to be loaded
390
       * @param bits  The number of bits to be loaded. This must not be more
391
       *              than the number of bytes in @p bytes.
392
       */
393
      bitvector_base(std::span<const uint8_t> bytes, /* NOLINT(*-explicit-conversions) */
80,942✔
394
                     std::optional<size_type> bits = std::nullopt) :
395
            m_bits() {
80,942✔
396
         from_bytes(bytes, bits);
80,942✔
397
      }
80,942✔
398

399
      bitvector_base(std::initializer_list<block_type> blocks, std::optional<size_type> bits = std::nullopt) :
196✔
400
            m_bits(bits.value_or(blocks.size() * block_size_bits)), m_blocks(blocks.begin(), blocks.end()) {}
196✔
401

402
      bool empty() const { return m_bits == 0; }
6✔
403

404
      size_type size() const { return m_bits; }
8,455,994✔
405

406
      /**
407
       * @returns true iff the number of 1-bits in this is odd, false otherwise (constant time)
408
       */
409
      CT::Choice has_odd_hamming_weight() const {
79,047✔
410
         uint64_t acc = 0;
79,047✔
411
         full_range_operation([&](std::unsigned_integral auto block) { acc ^= block; }, *this);
79,047✔
412

413
         for(size_t i = (sizeof(acc) * 8) >> 1; i > 0; i >>= 1) {
553,329✔
414
            acc ^= acc >> i;
474,282✔
415
         }
416

417
         return CT::Choice::from_int(acc & one);
79,047✔
418
      }
419

420
      /**
421
       * Counts the number of 1-bits in the bitvector in constant time.
422
       * @returns the "population count" (or hamming weight) of the bitvector
423
       */
424
      size_type hamming_weight() const {
155✔
425
         size_type acc = 0;
155✔
426
         full_range_operation([&](std::unsigned_integral auto block) { acc += ct_popcount(block); }, *this);
93✔
427
         return acc;
138✔
428
      }
429

430
      /**
431
       * @returns copies this bitvector into a new bitvector of type @p OutT
432
       */
433
      template <bitvectorish OutT>
434
      OutT as() const {
135✔
435
         return subvector<OutT>(0, size());
234✔
436
      }
437

438
      /**
439
       * @returns true if @p other contains the same bit pattern as this
440
       */
441
      template <bitvectorish OtherT>
442
      bool equals_vartime(const OtherT& other) const noexcept {
12✔
443
         return size() == other.size() &&
14✔
444
                full_range_operation([]<std::unsigned_integral BlockT>(BlockT lhs, BlockT rhs) { return lhs == rhs; },
14✔
445
                                     *this,
446
                                     unwrap_strong_type(other));
447
      }
448

449
      /**
450
       * @returns true if @p other contains the same bit pattern as this
451
       */
452
      template <bitvectorish OtherT>
453
      bool equals(const OtherT& other) const noexcept {
47✔
454
         return (*this ^ other).none();
47✔
455
      }
456

457
      /// @name Serialization
458
      /// @{
459

460
      /**
461
       * Re-initialize the bitvector with the given bytes. See the respective
462
       * constructor for details. This should be used only when trying to save
463
       * allocations. Otherwise, use the constructor.
464
       *
465
       * @param bytes  the byte range to load bits from
466
       * @param bits   (optional) if not all @p bytes should be loaded in full
467
       */
468
      void from_bytes(std::span<const uint8_t> bytes, std::optional<size_type> bits = std::nullopt) {
80,942✔
469
         const size_type new_bits = bits.has_value()
161,884✔
470
                                       ? bits.value()
80,942✔
471
                                       : mul_or_throw<size_t>(8, bytes.size_bytes(), "bitvector input is too large");
1,766✔
472
         const size_type bytes_needed = (new_bits / 8) + (new_bits % 8 != 0 ? 1 : 0);
143,306✔
473
         BOTAN_ARG_CHECK(bytes_needed <= bytes.size_bytes(), "not enough data to load so many bits");
80,942✔
474
         resize(new_bits);
80,942✔
475

476
         // load as much aligned data as possible
477
         const auto verbatim_blocks = new_bits / block_size_bits;
80,942✔
478
         const auto verbatim_bytes = verbatim_blocks * block_size_bytes;
80,942✔
479
         if(verbatim_blocks > 0) {
80,942✔
480
            typecast_copy(std::span{m_blocks}.first(verbatim_blocks), bytes.first(verbatim_bytes));
80,941✔
481
         }
482

483
         // load remaining unaligned data
484
         for(size_type i = verbatim_bytes * 8; i < new_bits; ++i) {
173,806✔
485
            ref(i) = ((bytes[i >> 3] & (uint8_t(1) << (i & 7))) != 0);
92,864✔
486
         }
487
      }
80,942✔
488

489
      /**
490
       * Renders the bitvector into a byte array. By default, this will use
491
       * `std::vector<uint8_t>` or `Botan::secure_vector<uint8_t>`, depending on
492
       * the allocator used by the bitvector. The rendering is compatible with
493
       * the bit layout explained in the respective constructor.
494
       */
495
      template <concepts::resizable_byte_buffer OutT =
496
                   std::conditional_t<uses_secure_allocator, secure_vector<uint8_t>, std::vector<uint8_t>>>
497
      OutT to_bytes() const {
328✔
498
         OutT out(ceil_tobytes(m_bits));
328✔
499
         to_bytes(out);
328✔
500
         return out;
328✔
501
      }
×
502

503
      /**
504
       * Renders the bitvector into a properly sized byte range.
505
       *
506
       * @param out  a byte range that has a length of at least `ceil_tobytes(size())`.
507
       */
508
      void to_bytes(std::span<uint8_t> out) const {
128,535✔
509
         const auto bytes_needed = ceil_tobytes(m_bits);
128,535✔
510
         BOTAN_ARG_CHECK(bytes_needed <= out.size_bytes(), "Not enough space to render bitvector");
128,535✔
511

512
         // copy as much aligned data as possible
513
         const auto verbatim_blocks = m_bits / block_size_bits;
128,535✔
514
         const auto verbatim_bytes = verbatim_blocks * block_size_bytes;
128,535✔
515
         if(verbatim_blocks > 0) {
128,535✔
516
            typecast_copy(out.first(verbatim_bytes), std::span{m_blocks}.first(verbatim_blocks));
128,532✔
517
         }
518

519
         // copy remaining unaligned data
520
         clear_mem(out.subspan(verbatim_bytes));
128,535✔
521
         for(size_type i = verbatim_bytes * 8; i < m_bits; ++i) {
321,954✔
522
            out[i >> 3] |= ref(i).template as<uint8_t>() << (i & 7);
193,419✔
523
         }
524
      }
128,535✔
525

526
      /**
527
       * Renders this bitvector into a sequence of "0"s and "1"s.
528
       * This is meant for debugging purposes and is not efficient.
529
       */
530
      std::string to_string() const {
1✔
531
         std::stringstream ss;
1✔
532
         for(size_type i = 0; i < size(); ++i) {
6✔
533
            ss << ref(i);
5✔
534
         }
535
         return ss.str();
1✔
536
      }
1✔
537

538
      /// @}
539

540
      /// @name Capacity Accessors and Modifiers
541
      /// @{
542

543
      size_type capacity() const { return m_blocks.capacity() * block_size_bits; }
6✔
544

545
      void reserve(size_type bits) { m_blocks.reserve(ceil_toblocks(bits)); }
219✔
546

547
      void resize(size_type bits) {
478,325✔
548
         const auto new_number_of_blocks = ceil_toblocks(bits);
478,325✔
549
         if(new_number_of_blocks != m_blocks.size()) {
478,325✔
550
            m_blocks.resize(new_number_of_blocks);
130,684✔
551
         }
552

553
         m_bits = bits;
478,325✔
554
         zero_unused_bits();
478,325✔
555
      }
478,325✔
556

557
      void push_back(bool bit) {
397,295✔
558
         const auto i = size();
397,295✔
559
         resize(i + 1);
397,295✔
560
         ref(i) = bit;
794,590✔
561
      }
397,295✔
562

563
      void pop_back() {
9✔
564
         if(!empty()) {
9✔
565
            resize(size() - 1);
9✔
566
         }
567
      }
568

569
      /// @}
570

571
      /// @name Bitwise and Global Accessors and Modifiers
572
      /// @{
573

574
      auto at(size_type pos) {
451,893,388✔
575
         check_offset(pos);
451,893,384✔
576
         return ref(pos);
451,893,144✔
577
      }
578

579
      // TODO C++23: deducing this
580
      auto at(size_type pos) const {
890,817✔
581
         check_offset(pos);
890,816✔
582
         return ref(pos);
890,816✔
583
      }
584

585
      auto front() { return ref(0); }
2✔
586

587
      // TODO C++23: deducing this
588
      auto front() const { return ref(0); }
589

590
      auto back() { return ref(size() - 1); }
5✔
591

592
      // TODO C++23: deducing this
593
      auto back() const { return ref(size() - 1); }
594

595
      /**
596
       * Sets the bit at position @p pos.
597
       * @throws Botan::Invalid_Argument if @p pos is out of range
598
       */
599
      bitvector_base& set(size_type pos) {
97✔
600
         check_offset(pos);
96✔
601
         ref(pos).set();
62✔
602
         return *this;
43✔
603
      }
604

605
      /**
606
       * Sets all currently allocated bits.
607
       */
608
      bitvector_base& set() {
2✔
609
         full_range_operation(
2✔
610
            [](std::unsigned_integral auto block) -> decltype(block) {
2✔
611
               return static_cast<decltype(block)>(~static_cast<decltype(block)>(0));
612
            },
613
            *this);
614
         zero_unused_bits();
2✔
615
         return *this;
2✔
616
      }
617

618
      /**
619
       * Unsets the bit at position @p pos.
620
       * @throws Botan::Invalid_Argument if @p pos is out of range
621
       */
622
      bitvector_base& unset(size_type pos) {
6✔
623
         check_offset(pos);
5✔
624
         ref(pos).unset();
5✔
625
         return *this;
5✔
626
      }
627

628
      /**
629
       * Unsets all currently allocated bits.
630
       */
631
      bitvector_base& unset() {
4✔
632
         full_range_operation(
4✔
633
            [](std::unsigned_integral auto block) -> decltype(block) { return static_cast<decltype(block)>(0); },
4✔
634
            *this);
635
         return *this;
636
      }
637

638
      /**
639
       * Flips the bit at position @p pos.
640
       * @throws Botan::Invalid_Argument if @p pos is out of range
641
       */
642
      bitvector_base& flip(size_type pos) {
30✔
643
         check_offset(pos);
29✔
644
         ref(pos).flip();
29✔
645
         return *this;
29✔
646
      }
647

648
      /**
649
       * Flips all currently allocated bits.
650
       */
651
      bitvector_base& flip() {
5✔
652
         full_range_operation([](std::unsigned_integral auto block) -> decltype(block) { return ~block; }, *this);
5✔
653
         zero_unused_bits();
5✔
654
         return *this;
5✔
655
      }
656

657
      /**
658
       * @returns true iff no bit is set
659
       */
660
      bool none_vartime() const {
13✔
661
         return full_range_operation([](std::unsigned_integral auto block) { return block == 0; }, *this);
6✔
662
      }
663

664
      /**
665
       * @returns true iff no bit is set in constant time
666
       */
667
      bool none() const { return hamming_weight() == 0; }
52✔
668

669
      /**
670
       * @returns true iff at least one bit is set
671
       */
672
      bool any_vartime() const { return !none_vartime(); }
7✔
673

674
      /**
675
       * @returns true iff at least one bit is set in constant time
676
       */
677
      bool any() const { return !none(); }
5✔
678

679
      /**
680
       * @returns true iff all bits are set
681
       */
682
      bool all_vartime() const {
7✔
683
         return full_range_operation(
7✔
684
            []<std::unsigned_integral BlockT>(BlockT block, BlockT mask) { return block == mask; }, *this);
6✔
685
      }
686

687
      /**
688
       * @returns true iff all bits are set in constant time
689
       */
690
      bool all() const { return hamming_weight() == m_bits; }
5✔
691

692
      auto operator[](size_type pos) { return ref(pos); }
1,315,507,015✔
693

694
      // TODO C++23: deducing this
695
      auto operator[](size_type pos) const { return ref(pos); }
12,448,319✔
696

697
      /// @}
698

699
      /// @name Subvectors
700
      /// @{
701

702
      /**
703
       * Creates a new bitvector with a subsection of this bitvector starting at
704
       * @p pos copying exactly @p length bits.
705
       */
706
      template <bitvectorish OutT = bitvector_base<AllocatorT>>
707
      auto subvector(size_type pos, std::optional<size_type> length = std::nullopt) const {
128,423✔
708
         const size_type bitlen = length.value_or(size() - pos);
128,423✔
709
         BOTAN_ARG_CHECK(pos + bitlen <= size(), "Not enough bits to copy");
128,423✔
710

711
         OutT newvector(bitlen);
128,419✔
712

713
         // Handle bitvectors that are wrapped in strong types
714
         auto& newvector_unwrapped = unwrap_strong_type(newvector);
128,419✔
715

716
         if(bitlen > 0) {
128,419✔
717
            if(pos % 8 == 0) {
128,416✔
718
               copy_mem(
89,722✔
719
                  newvector_unwrapped.m_blocks,
89,722✔
720
                  std::span{m_blocks}.subspan(block_index(pos), block_index(pos + bitlen - 1) - block_index(pos) + 1));
89,722✔
721
            } else {
722
               const BitRangeOperator<const bitvector_base<AllocatorT>, BitRangeAlignment::no_alignment> from_op(
38,694✔
723
                  *this, pos, bitlen);
724
               const BitRangeOperator<strong_type_wrapped_type<OutT>> to_op(
38,694✔
725
                  unwrap_strong_type(newvector_unwrapped), 0, bitlen);
726
               range_operation([](auto /* to */, auto from) { return from; }, to_op, from_op);
38,694✔
727
            }
728

729
            newvector_unwrapped.zero_unused_bits();
128,419✔
730
         }
731

732
         return newvector;
128,419✔
733
      }
×
734

735
      /**
736
       * Extracts a subvector of bits as an unsigned integral type @p OutT
737
       * starting from bit @p pos and copying exactly sizeof(OutT)*8 bits.
738
       *
739
       * Hint: The bits are in big-endian order, i.e. the least significant bit
740
       *       is the 0th bit and the most significant bit it the n-th. Hence,
741
       *       addressing the bits with bitwise operations is done like so:
742
       *       bool bit = (out_int >> pos) & 1;
743
       */
744
      template <typename OutT>
745
         requires(std::unsigned_integral<strong_type_wrapped_type<OutT>> &&
746
                  !std::same_as<bool, strong_type_wrapped_type<OutT>>)
747
      OutT subvector(size_type pos) const {
74,960✔
748
         using result_t = strong_type_wrapped_type<OutT>;
749
         constexpr size_t bits = sizeof(result_t) * 8;
74,960✔
750
         BOTAN_ARG_CHECK(pos + bits <= size(), "Not enough bits to copy");
74,960✔
751
         result_t out = 0;
74,956✔
752

753
         if(pos % 8 == 0) {
74,956✔
754
            out = load_le<result_t>(std::span{m_blocks}.subspan(block_index(pos)).template first<sizeof(result_t)>());
48,648✔
755
         } else {
756
            const BitRangeOperator<const bitvector_base<AllocatorT>, BitRangeAlignment::no_alignment> op(
26,308✔
757
               *this, pos, bits);
758
            range_operation(
26,308✔
759
               [&](std::unsigned_integral auto integer) {
26,308✔
760
                  if constexpr(std::same_as<result_t, decltype(integer)>) {
761
                     out = integer;
26,308✔
762
                  }
763
               },
764
               op);
765
         }
766

767
         return wrap_strong_type<OutT>(out);
74,956✔
768
      }
769

770
      /**
771
       * Replaces a subvector of bits with the bits of another bitvector @p value
772
       * starting at bit @p pos. The number of bits to replace is determined by
773
       * the size of @p value.
774
       *
775
       * @note This is currently supported for byte-aligned @p pos only.
776
       *
777
       * @throws Not_Implemented when called with @p pos not divisible by 8.
778
       *
779
       * @param pos    the position to start replacing bits
780
       * @param value  the bitvector to copy bits from
781
       */
782
      template <typename InT>
783
         requires(std::unsigned_integral<strong_type_wrapped_type<InT>> && !std::same_as<bool, InT>)
784
      void subvector_replace(size_type pos, InT value) {
74,967✔
785
         using in_t = strong_type_wrapped_type<InT>;
786
         constexpr size_t bits = sizeof(in_t) * 8;
74,967✔
787
         BOTAN_ARG_CHECK(pos + bits <= size(), "Not enough bits to replace");
74,967✔
788

789
         if(pos % 8 == 0) {
74,963✔
790
            store_le(std::span{m_blocks}.subspan(block_index(pos)).template first<sizeof(in_t)>(),
48,650✔
791
                     unwrap_strong_type(value));
792
         } else {
793
            const BitRangeOperator<bitvector_base<AllocatorT>, BitRangeAlignment::no_alignment> op(*this, pos, bits);
26,313✔
794
            range_operation(
26,313✔
795
               [&]<std::unsigned_integral BlockT>(BlockT block) -> BlockT {
26,313✔
796
                  if constexpr(std::same_as<in_t, BlockT>) {
797
                     return unwrap_strong_type(value);
26,313✔
798
                  } else {
799
                     // This should never be reached. BOTAN_ASSERT_UNREACHABLE()
800
                     // caused warning "unreachable code" on MSVC, though. You
801
                     // don't say!
802
                     //
803
                     // Returning the given block back, is the most reasonable
804
                     // thing to do in this case, though.
805
                     return block;
806
                  }
807
               },
808
               op);
809
         }
810
      }
74,963✔
811

812
      /// @}
813

814
      /// @name Operators
815
      ///
816
      /// @{
817

818
      auto operator~() {
1✔
819
         auto newbv = *this;
1✔
820
         newbv.flip();
1✔
821
         return newbv;
1✔
822
      }
×
823

824
      template <bitvectorish OtherT>
825
      auto& operator|=(const OtherT& other) {
6✔
826
         full_range_operation([]<std::unsigned_integral BlockT>(BlockT lhs, BlockT rhs) -> BlockT { return lhs | rhs; },
2✔
827
                              *this,
828
                              unwrap_strong_type(other));
829
         return *this;
830
      }
831

832
      template <bitvectorish OtherT>
833
      auto& operator&=(const OtherT& other) {
79,049✔
834
         full_range_operation([]<std::unsigned_integral BlockT>(BlockT lhs, BlockT rhs) -> BlockT { return lhs & rhs; },
79,045✔
835
                              *this,
836
                              unwrap_strong_type(other));
837
         return *this;
838
      }
839

840
      template <bitvectorish OtherT>
841
      auto& operator^=(const OtherT& other) {
55✔
842
         full_range_operation([]<std::unsigned_integral BlockT>(BlockT lhs, BlockT rhs) -> BlockT { return lhs ^ rhs; },
136✔
843
                              *this,
844
                              unwrap_strong_type(other));
845
         return *this;
846
      }
847

848
      /// @}
849

850
      /// @name Constant Time Operations
851
      ///
852
      /// @{
853

854
      /**
855
       * Implements::
856
       *
857
       *    if(condition) {
858
       *       *this ^= other;
859
       *    }
860
       *
861
       * omitting runtime dependence on any of the parameters.
862
       */
863
      template <bitvectorish OtherT>
864
      void ct_conditional_xor(CT::Choice condition, const OtherT& other) {
447,156,718✔
865
         BOTAN_ASSERT_NOMSG(m_bits == other.m_bits);
447,156,718✔
866
         BOTAN_ASSERT_NOMSG(m_blocks.size() == other.m_blocks.size());
447,156,718✔
867

868
         auto maybe_xor = overloaded{
869
            [m = CT::Mask<uint64_t>::from_choice(condition)](uint64_t lhs, uint64_t rhs) -> uint64_t {
2,147,483,647✔
870
               return lhs ^ m.if_set_return(rhs);
2,147,483,647✔
871
            },
872
            [m = CT::Mask<uint32_t>::from_choice(condition)](uint32_t lhs, uint32_t rhs) -> uint32_t {
718,244,525✔
873
               return lhs ^ m.if_set_return(rhs);
271,087,807✔
874
            },
875
            [m = CT::Mask<uint16_t>::from_choice(condition)](uint16_t lhs, uint16_t rhs) -> uint16_t {
569,109,836✔
876
               return lhs ^ m.if_set_return(rhs);
121,953,118✔
877
            },
878
            [m = CT::Mask<uint8_t>::from_choice(condition)](uint8_t lhs, uint8_t rhs) -> uint8_t {
447,156,724✔
879
               return lhs ^ m.if_set_return(rhs);
6✔
880
            },
881
         };
882

883
         full_range_operation(maybe_xor, *this, unwrap_strong_type(other));
447,156,718✔
884
      }
447,156,718✔
885

886
      constexpr void _const_time_poison() const { CT::poison(m_blocks); }
71✔
887

888
      constexpr void _const_time_unpoison() const { CT::unpoison(m_blocks); }
98✔
889

890
      /// @}
891

892
      /// @name Iterators
893
      ///
894
      /// @{
895

896
      iterator begin() noexcept { return iterator(this, 0); }
1,067✔
897

898
      const_iterator begin() const noexcept { return const_iterator(this, 0); }
10✔
899

900
      const_iterator cbegin() const noexcept { return const_iterator(this, 0); }
10✔
901

902
      iterator end() noexcept { return iterator(this, size()); }
537✔
903

904
      const_iterator end() const noexcept { return const_iterator(this, size()); }
5✔
905

906
      const_iterator cend() noexcept { return const_iterator(this, size()); }
8✔
907

908
      /// @}
909

910
   private:
911
      void check_offset(size_type pos) const {
452,783,959✔
912
         // BOTAN_ASSERT_NOMSG(!CT::is_poisoned(&m_bits, sizeof(m_bits)));
913
         // BOTAN_ASSERT_NOMSG(!CT::is_poisoned(&pos, sizeof(pos)));
914
         BOTAN_ARG_CHECK(pos < m_bits, "Out of range");
452,784,304✔
915
      }
916

917
      void zero_unused_bits() {
606,748✔
918
         const auto first_unused_bit = size();
606,748✔
919

920
         // Zero out any unused bits in the last block
921
         if(first_unused_bit % block_size_bits != 0) {
606,748✔
922
            const block_type mask = (one << block_offset(first_unused_bit)) - one;
404,956✔
923
            m_blocks[block_index(first_unused_bit)] &= mask;
404,956✔
924
         }
925
      }
926

927
      static constexpr size_type ceil_toblocks(size_type bits) {
607,169✔
928
         return add_or_throw(bits, block_size_bits - 1, "bitvector size is too large") / block_size_bits;
607,169✔
929
      }
930

931
      auto ref(size_type pos) const { return bitref<const block_type>(m_blocks, pos); }
12,641,743✔
932

933
      auto ref(size_type pos) { return bitref<block_type>(m_blocks, pos); }
1,324,445,183✔
934

935
   private:
936
      enum class BitRangeAlignment : uint8_t { byte_aligned, no_alignment };
937

938
      /**
939
       * Helper construction to implement bit range operations on the bitvector.
940
       * It basically implements an iterator to read and write blocks of bits
941
       * from the underlying bitvector. Where "blocks of bits" are unsigned
942
       * integers of varying bit lengths.
943
       *
944
       * If the iteration starts at a byte boundary in the underlying bitvector,
945
       * this applies certain optimizations (i.e. loading blocks of bits straight
946
       * from the underlying byte buffer). The optimizations are enabled at
947
       * compile time (with the template parameter `alignment`).
948
       */
949
      template <typename BitvectorT, auto alignment = BitRangeAlignment::byte_aligned>
950
         requires is_bitvector_v<std::remove_cvref_t<BitvectorT>>
951
      class BitRangeOperator {
952
         private:
953
            constexpr static bool is_const() { return std::is_const_v<BitvectorT>; }
954

955
            struct UnalignedDataHelper {
956
                  const uint8_t padding_bits;
957
                  const uint8_t bits_to_byte_alignment;
958
            };
959

960
         public:
961
            BitRangeOperator(BitvectorT& source, size_type start_bitoffset, size_type bitlength) :
894,680,922✔
962
                  m_source(source),
894,680,922✔
963
                  m_start_bitoffset(start_bitoffset),
894,680,922✔
964
                  m_bitlength(bitlength),
894,680,922✔
965
                  m_unaligned_helper({.padding_bits = static_cast<uint8_t>(start_bitoffset % 8),
894,680,922✔
966
                                      .bits_to_byte_alignment = static_cast<uint8_t>(8 - (start_bitoffset % 8))}),
894,680,922✔
967
                  m_read_bitpos(start_bitoffset),
894,680,922✔
968
                  m_write_bitpos(start_bitoffset) {
894,680,922✔
969
               BOTAN_ASSERT(is_byte_aligned() == (m_start_bitoffset % 8 == 0), "byte alignment guarantee");
894,680,922✔
970
               BOTAN_ASSERT(m_source.size() >= m_start_bitoffset + m_bitlength, "enough bytes in underlying source");
894,680,922✔
971
            }
894,680,922✔
972

973
            explicit BitRangeOperator(BitvectorT& source) : BitRangeOperator(source, 0, source.size()) {}
894,550,913✔
974

975
            static constexpr bool is_byte_aligned() { return alignment == BitRangeAlignment::byte_aligned; }
976

977
            /**
978
             * @returns the overall number of bits to be iterated with this operator
979
             */
980
            size_type size() const { return m_bitlength; }
447,274,534✔
981

982
            /**
983
             * @returns the number of bits not yet read from this operator
984
             */
985
            size_type bits_to_read() const { return m_bitlength - m_read_bitpos + m_start_bitoffset; }
1,796,388,135✔
986

987
            /**
988
             * @returns the number of bits still to be written into this operator
989
             */
990
            size_type bits_to_write() const { return m_bitlength - m_write_bitpos + m_start_bitoffset; }
3,398,320✔
991

992
            /**
993
             * Loads the next block of bits from the underlying bitvector. No
994
             * bounds checks are performed. The caller can define the size of
995
             * the resulting unsigned integer block.
996
             */
997
            template <std::unsigned_integral BlockT>
998
            BlockT load_next() const {
6,762,620✔
999
               constexpr size_type block_size = sizeof(BlockT);
6,762,620✔
1000
               constexpr size_type block_bits = block_size * 8;
6,762,620✔
1001
               const auto bits_remaining = bits_to_read();
6,762,620✔
1002

1003
               BlockT result_block = 0;
6,762,620✔
1004
               if constexpr(is_byte_aligned()) {
1005
                  result_block = load_le(m_source.as_byte_span().subspan(read_bytepos()).template first<block_size>());
3,382,906✔
1006
               } else {
1007
                  const size_type byte_pos = read_bytepos();
3,379,714✔
1008
                  const size_type bits_to_collect = std::min(block_bits, bits_to_read());
6,759,428✔
1009

1010
                  const uint8_t first_byte = m_source.as_byte_span()[byte_pos];
3,379,714✔
1011

1012
                  // Initialize the left-most bits from the first byte.
1013
                  result_block = BlockT(first_byte) >> m_unaligned_helper.padding_bits;
3,379,714✔
1014

1015
                  // If more bits are needed, we pull them from the remaining bytes.
1016
                  if(m_unaligned_helper.bits_to_byte_alignment < bits_to_collect) {
3,379,714✔
1017
                     const BlockT block =
1018
                        load_le(m_source.as_byte_span().subspan(byte_pos + 1).template first<block_size>());
6,682,042✔
1019
                     result_block |= block << m_unaligned_helper.bits_to_byte_alignment;
3,341,024✔
1020
                  }
1021
               }
1022

1023
               m_read_bitpos += std::min(block_bits, bits_remaining);
6,762,620✔
1024
               return result_block;
6,668,117✔
1025
            }
1026

1027
            /**
1028
             * Stores the next block of bits into the underlying bitvector.
1029
             * No bounds checks are performed. Storing bit blocks that are not
1030
             * aligned at a byte-boundary in the underlying bitvector is
1031
             * currently not implemented.
1032
             */
1033
            template <std::unsigned_integral BlockT>
1034
               requires(!is_const())
1035
            void store_next(BlockT block) {
3,372,007✔
1036
               constexpr size_type block_size = sizeof(BlockT);
3,372,007✔
1037
               constexpr size_type block_bits = block_size * 8;
3,372,007✔
1038

1039
               if constexpr(is_byte_aligned()) {
1040
                  auto sink = m_source.as_byte_span().subspan(write_bytepos()).template first<block_size>();
3,345,694✔
1041
                  store_le(sink, block);
3,345,694✔
1042
               } else {
1043
                  const size_type byte_pos = write_bytepos();
26,313✔
1044
                  const size_type bits_to_store = std::min(block_bits, bits_to_write());
52,626✔
1045

1046
                  uint8_t& first_byte = m_source.as_byte_span()[byte_pos];
26,313✔
1047

1048
                  // Set the left-most bits in the first byte, leaving all others unchanged
1049
                  first_byte = (first_byte & uint8_t(0xFF >> m_unaligned_helper.bits_to_byte_alignment)) |
26,313✔
1050
                               uint8_t(block << m_unaligned_helper.padding_bits);
26,313✔
1051

1052
                  // If more bits are provided, we store them in the remaining bytes.
1053
                  if(m_unaligned_helper.bits_to_byte_alignment < bits_to_store) {
26,313✔
1054
                     const auto remaining_bytes =
1055
                        m_source.as_byte_span().subspan(byte_pos + 1).template first<block_size>();
26,313✔
1056
                     const BlockT padding_mask = ~(BlockT(-1) >> m_unaligned_helper.bits_to_byte_alignment);
26,313✔
1057
                     const BlockT new_bytes =
26,313✔
1058
                        (load_le(remaining_bytes) & padding_mask) | block >> m_unaligned_helper.bits_to_byte_alignment;
26,313✔
1059
                     store_le(remaining_bytes, new_bytes);
26,313✔
1060
                  }
1061
               }
1062

1063
               m_write_bitpos += std::min(block_bits, bits_to_write());
3,372,007✔
1064
            }
3,372,007✔
1065

1066
            template <std::unsigned_integral BlockT>
1067
               requires(is_byte_aligned() && !is_const())
1068
            std::span<BlockT> span(size_type blocks) const {
1,341,707,517✔
1069
               BOTAN_DEBUG_ASSERT(blocks == 0 || is_memory_aligned_to<BlockT>());
1,341,707,517✔
1070
               BOTAN_DEBUG_ASSERT(read_bytepos() % sizeof(BlockT) == 0);
1,341,707,517✔
1071
               // Intermittently casting to void* to avoid a compiler warning
1072
               void* ptr = reinterpret_cast<void*>(m_source.as_byte_span().data() + read_bytepos());
1,341,707,517✔
1073
               return {reinterpret_cast<BlockT*>(ptr), blocks};
1,341,707,517✔
1074
            }
1075

1076
            template <std::unsigned_integral BlockT>
1077
               requires(is_byte_aligned() && is_const())
1078
            std::span<const BlockT> span(size_type blocks) const {
1,341,945,163✔
1079
               BOTAN_DEBUG_ASSERT(blocks == 0 || is_memory_aligned_to<BlockT>());
1,341,945,163✔
1080
               BOTAN_DEBUG_ASSERT(read_bytepos() % sizeof(BlockT) == 0);
1,341,945,163✔
1081
               // Intermittently casting to void* to avoid a compiler warning
1082
               const void* ptr = reinterpret_cast<const void*>(m_source.as_byte_span().data() + read_bytepos());
1,341,945,163✔
1083
               return {reinterpret_cast<const BlockT*>(ptr), blocks};
1,341,945,163✔
1084
            }
1085

1086
            void advance(size_type bytes)
1,341,945,196✔
1087
               requires(is_byte_aligned())
1088
            {
1089
               m_read_bitpos += bytes * 8;
1,341,945,196✔
1090
               m_write_bitpos += bytes * 8;
1,341,945,196✔
1091
            }
1092

1093
            template <std::unsigned_integral BlockT>
1094
               requires(is_byte_aligned())
1095
            size_t is_memory_aligned_to() const {
894,550,913✔
1096
               const void* cptr = m_source.as_byte_span().data() + read_bytepos();
894,550,913✔
1097
               const void* ptr_before = cptr;
894,550,913✔
1098

1099
               // std::align takes `ptr` as a reference (!), i.e. `void*&` and
1100
               // uses it as an out-param. Though, `cptr` is const because this
1101
               // method is const-qualified, hence the const_cast<>.
1102
               void* ptr = const_cast<void*>(cptr);  // NOLINT(*-const-correctness)
894,550,913✔
1103
               size_t size = sizeof(BlockT);
894,550,913✔
1104
               return ptr_before != nullptr && std::align(alignof(BlockT), size, ptr, size) == ptr_before;
1,789,101,826✔
1105
            }
1106

1107
         private:
1108
            size_type read_bytepos() const { return m_read_bitpos / 8; }
901,313,516✔
1109

1110
            size_type write_bytepos() const { return m_write_bitpos / 8; }
3,372,007✔
1111

1112
         private:
1113
            BitvectorT& m_source;
1114
            size_type m_start_bitoffset;
1115
            size_type m_bitlength;
1116

1117
            UnalignedDataHelper m_unaligned_helper;
1118

1119
            mutable size_type m_read_bitpos;
1120
            mutable size_type m_write_bitpos;
1121
      };
1122

1123
      /**
1124
       * Helper struct for the low-level handling of blockwise operations
1125
       *
1126
       * This has two main code paths: Optimized for byte-aligned ranges that
1127
       * can simply be taken from memory as-is. And a generic implementation
1128
       * that must assemble blocks from unaligned bits before processing.
1129
       */
1130
      template <typename FnT, typename... ParamTs>
1131
         requires detail::blockwise_processing_callback<FnT, ParamTs...>
1132
      class blockwise_processing_callback_trait {
1133
         public:
1134
            constexpr static bool needs_mask = detail::blockwise_processing_callback_with_mask<FnT, ParamTs...>;
1135
            constexpr static bool is_manipulator = detail::manipulating_blockwise_processing_callback<FnT, ParamTs...>;
1136
            constexpr static bool is_predicate = detail::predicate_blockwise_processing_callback<FnT, ParamTs...>;
1137
            static_assert(!is_manipulator || !is_predicate, "cannot be manipulator and predicate at the same time");
1138

1139
            /**
1140
             * Applies @p fn to the blocks provided in @p blocks by simply reading from
1141
             * memory without re-arranging any bits across byte-boundaries.
1142
             */
1143
            template <std::unsigned_integral... BlockTs>
1144
               requires(all_same_v<std::remove_cv_t<BlockTs>...> && sizeof...(BlockTs) == sizeof...(ParamTs))
1145
            constexpr static bool apply_on_full_blocks(FnT fn, std::span<BlockTs>... blocks) {
1,341,945,196✔
1146
               constexpr size_type bits = sizeof(detail::first_t<BlockTs...>) * 8;
1,341,945,196✔
1147
               const size_type iterations = detail::first(blocks...).size();
1,341,707,949✔
1148
               for(size_type i = 0; i < iterations; ++i) {
2,147,483,647✔
1149
                  if constexpr(is_predicate) {
1150
                     if(!apply(fn, bits, blocks[i]...)) {
33✔
1151
                        return false;
1152
                     }
1153
                  } else if constexpr(is_manipulator) {
1154
                     detail::first(blocks...)[i] = apply(fn, bits, blocks[i]...);
2,147,483,647✔
1155
                  } else {
1156
                     apply(fn, bits, blocks[i]...);
6,156,777✔
1157
                  }
1158
               }
1159
               return true;
465✔
1160
            }
1161

1162
            /**
1163
             * Applies @p fn to as many blocks as @p ops provide for the given type.
1164
             */
1165
            template <std::unsigned_integral BlockT, typename... BitRangeOperatorTs>
1166
               requires(sizeof...(BitRangeOperatorTs) == sizeof...(ParamTs))
1167
            constexpr static bool apply_on_unaligned_blocks(FnT fn, BitRangeOperatorTs&... ops) {
447,680,319✔
1168
               constexpr size_type block_bits = sizeof(BlockT) * 8;
447,680,319✔
1169
               auto bits = detail::first(ops...).bits_to_read();
447,680,319✔
1170
               if(bits == 0) {
447,680,319✔
1171
                  return true;
1172
               }
1173

1174
               bits += block_bits;  // avoid unsigned integer underflow in the following loop
244,635✔
1175
               while(bits > block_bits * 2 - 8) {
3,661,557✔
1176
                  bits -= block_bits;
3,416,925✔
1177
                  if constexpr(is_predicate) {
1178
                     if(!apply(fn, bits, ops.template load_next<BlockT>()...)) {
29✔
1179
                        return false;
1180
                     }
1181
                  } else if constexpr(is_manipulator) {
1182
                     detail::first(ops...).store_next(apply(fn, bits, ops.template load_next<BlockT>()...));
3,410,697✔
1183
                  } else {
1184
                     apply(fn, bits, ops.template load_next<BlockT>()...);
44,900✔
1185
                  }
1186
               }
1187
               return true;
1188
            }
1189

1190
         private:
1191
            template <std::unsigned_integral... BlockTs>
1192
               requires(all_same_v<std::remove_cv_t<BlockTs>...>)
1193
            constexpr static auto apply(FnT fn, size_type bits, BlockTs... blocks) {
2,147,483,647✔
1194
               if constexpr(needs_mask) {
1195
                  return fn(blocks..., make_mask<detail::first_t<BlockTs...>>(bits));
3✔
1196
               } else {
1197
                  return fn(blocks...);
2,147,483,647✔
1198
               }
1199
            }
1200
      };
1201

1202
      /**
1203
       * Helper function of `full_range_operation` and `range_operation` that
1204
       * calls @p fn on a given aligned unsigned integer block as long as the
1205
       * underlying bit range contains enough bits to fill the block fully.
1206
       *
1207
       * This uses bare memory access to gain a speed up for aligned data.
1208
       */
1209
      template <std::unsigned_integral BlockT, typename FnT, typename... BitRangeOperatorTs>
1210
         requires(detail::blockwise_processing_callback<FnT, BitRangeOperatorTs...> &&
1211
                  sizeof...(BitRangeOperatorTs) > 0)
1212
      static bool _process_in_fully_aligned_blocks_of(FnT fn, BitRangeOperatorTs&... ops) {
1,341,945,196✔
1213
         constexpr size_type block_bytes = sizeof(BlockT);
1,341,945,196✔
1214
         constexpr size_type block_bits = block_bytes * 8;
1,341,945,196✔
1215
         const size_type blocks = detail::first(ops...).bits_to_read() / block_bits;
1,341,945,196✔
1216

1217
         using callback_trait = blockwise_processing_callback_trait<FnT, BitRangeOperatorTs...>;
1218
         const auto result = callback_trait::apply_on_full_blocks(fn, ops.template span<BlockT>(blocks)...);
2,147,483,647✔
1219
         (ops.advance(block_bytes * blocks), ...);
1,341,945,196✔
1220
         return result;
1,341,945,196✔
1221
      }
1222

1223
      /**
1224
       * Helper function of `full_range_operation` and `range_operation` that
1225
       * calls @p fn on a given unsigned integer block size as long as the
1226
       * underlying bit range contains enough bits to fill the block.
1227
       */
1228
      template <std::unsigned_integral BlockT, typename FnT, typename... BitRangeOperatorTs>
1229
         requires(detail::blockwise_processing_callback<FnT, BitRangeOperatorTs...>)
1230
      static bool _process_in_unaligned_blocks_of(FnT fn, BitRangeOperatorTs&... ops) {
447,680,319✔
1231
         using callback_trait = blockwise_processing_callback_trait<FnT, BitRangeOperatorTs...>;
1232
         return callback_trait::template apply_on_unaligned_blocks<BlockT>(fn, ops...);
447,589,004✔
1233
      }
1234

1235
      /**
1236
       * Apply @p fn to all bits in the ranges defined by @p ops. If more than
1237
       * one range operator is passed to @p ops, @p fn receives corresponding
1238
       * blocks of bits from each operator. Therefore, all @p ops have to define
1239
       * the exact same length of their underlying ranges.
1240
       *
1241
       * @p fn may return a bit block that will be stored into the _first_ bit
1242
       * range passed into @p ops. If @p fn returns a boolean, and its value is
1243
       * `false`, the range operation is cancelled and `false` is returned.
1244
       *
1245
       * The implementation ensures to pull bits in the largest bit blocks
1246
       * possible and reverts to smaller bit blocks only when needed.
1247
       */
1248
      template <typename FnT, typename... BitRangeOperatorTs>
1249
         requires(detail::blockwise_processing_callback<FnT, BitRangeOperatorTs...> &&
1250
                  sizeof...(BitRangeOperatorTs) > 0)
1251
      static bool range_operation(FnT fn, BitRangeOperatorTs... ops) {
447,406,388✔
1252
         BOTAN_ASSERT(has_equal_lengths(ops...), "all BitRangeOperators have the same length");
170,548✔
1253

1254
         if constexpr((BitRangeOperatorTs::is_byte_aligned() && ...)) {
1255
            // Note: At the moment we can assume that this will always be used
1256
            //       on the _entire_ bitvector. Therefore, we can safely assume
1257
            //       that the bitvectors' underlying buffers are properly aligned.
1258
            //       If this assumption changes, we need to add further handling
1259
            //       to process a byte padding at the beginning of the bitvector
1260
            //       until a memory alignment boundary is reached.
1261
            const bool alignment = (ops.template is_memory_aligned_to<uint64_t>() && ...);
894,550,913✔
1262
            BOTAN_ASSERT_NOMSG(alignment);
×
1263

1264
            return _process_in_fully_aligned_blocks_of<uint64_t>(fn, ops...) &&
894,630,136✔
1265
                   _process_in_fully_aligned_blocks_of<uint32_t>(fn, ops...) &&
894,630,123✔
1266
                   _process_in_fully_aligned_blocks_of<uint16_t>(fn, ops...) &&
894,630,133✔
1267
                   _process_in_unaligned_blocks_of<uint8_t>(fn, ops...);
447,315,059✔
1268
         } else {
1269
            return _process_in_unaligned_blocks_of<uint64_t>(fn, ops...) &&
182,630✔
1270
                   _process_in_unaligned_blocks_of<uint32_t>(fn, ops...) &&
91,315✔
1271
                   _process_in_unaligned_blocks_of<uint16_t>(fn, ops...) &&
182,630✔
1272
                   _process_in_unaligned_blocks_of<uint8_t>(fn, ops...);
91,315✔
1273
         }
1274
      }
1275

1276
      /**
1277
       * Apply @p fn to all bit blocks in the bitvector(s).
1278
       */
1279
      template <typename FnT, typename... BitvectorTs>
1280
         requires(detail::blockwise_processing_callback<FnT, BitvectorTs...> &&
1281
                  (is_bitvector_v<std::remove_cvref_t<BitvectorTs>> && ... && true))
1282
      static bool full_range_operation(FnT&& fn, BitvectorTs&... bitvecs) {
447,315,073✔
1283
         BOTAN_ASSERT(has_equal_lengths(bitvecs...), "all bitvectors have the same length");
447,315,073✔
1284
         return range_operation(std::forward<FnT>(fn), BitRangeOperator<BitvectorTs>(bitvecs)...);
447,315,073✔
1285
      }
1286

1287
      template <typename SomeT, typename... SomeTs>
1288
      static bool has_equal_lengths(const SomeT& v, const SomeTs&... vs) {
894,510,374✔
1289
         return ((v.size() == vs.size()) && ... && true);
894,510,374✔
1290
      }
1291

1292
      template <std::unsigned_integral T>
1293
      static constexpr T make_mask(size_type bits) {
3✔
1294
         const bool max = bits >= sizeof(T) * 8;
3✔
1295
         bits &= T(max - 1);
3✔
1296
         return (T(!max) << bits) - 1;
3✔
1297
      }
1298

1299
      auto as_byte_span() { return std::span{m_blocks.data(), m_blocks.size() * sizeof(block_type)}; }
1,795,687,361✔
1300

1301
      auto as_byte_span() const { return std::span{m_blocks.data(), m_blocks.size() * sizeof(block_type)}; }
1,792,650,842✔
1302

1303
   private:
1304
      size_type m_bits;
1305
      std::vector<block_type, allocator_type> m_blocks;
1306
};
1307

1308
using secure_bitvector = bitvector_base<secure_allocator>;
1309
using bitvector = bitvector_base<std::allocator>;
1310

1311
namespace detail {
1312

1313
/**
1314
 * If one of the allocators is a Botan::secure_allocator, this will always
1315
 * prefer it. Otherwise, the allocator of @p lhs will be used as a default.
1316
 */
1317
template <bitvectorish T1, bitvectorish T2>
1318
constexpr auto copy_lhs_allocator_aware(const T1& lhs, const T2& /*rhs*/) {
63✔
1319
   constexpr bool needs_secure_allocator =
63✔
1320
      strong_type_wrapped_type<T1>::uses_secure_allocator || strong_type_wrapped_type<T2>::uses_secure_allocator;
1321

1322
   if constexpr(needs_secure_allocator) {
1323
      return lhs.template as<secure_bitvector>();
56✔
1324
   } else {
1325
      return lhs.template as<bitvector>();
8✔
1326
   }
1327
}
1328

1329
}  // namespace detail
1330

1331
template <bitvectorish T1, bitvectorish T2>
1332
auto operator|(const T1& lhs, const T2& rhs) {
5✔
1333
   auto res = detail::copy_lhs_allocator_aware(lhs, rhs);
5✔
1334
   res |= rhs;
5✔
1335
   return res;
5✔
1336
}
×
1337

1338
template <bitvectorish T1, bitvectorish T2>
1339
auto operator&(const T1& lhs, const T2& rhs) {
4✔
1340
   auto res = detail::copy_lhs_allocator_aware(lhs, rhs);
4✔
1341
   res &= rhs;
4✔
1342
   return res;
4✔
1343
}
×
1344

1345
template <bitvectorish T1, bitvectorish T2>
1346
auto operator^(const T1& lhs, const T2& rhs) {
54✔
1347
   auto res = detail::copy_lhs_allocator_aware(lhs, rhs);
54✔
1348
   res ^= rhs;
54✔
1349
   return res;
54✔
1350
}
×
1351

1352
template <bitvectorish T1, bitvectorish T2>
1353
bool operator==(const T1& lhs, const T2& rhs) {
10✔
1354
   return lhs.equals_vartime(rhs);
10✔
1355
}
1356

1357
namespace detail {
1358

1359
/**
1360
 * A Strong<> adapter for arbitrarily large bitvectors
1361
 */
1362
template <concepts::container T>
1363
   requires is_bitvector_v<T>
1364
class Strong_Adapter<T> : public Container_Strong_Adapter_Base<T> {
612,985✔
1365
   public:
1366
      using size_type = typename T::size_type;
1367

1368
   public:
1369
      using Container_Strong_Adapter_Base<T>::Container_Strong_Adapter_Base;
597✔
1370

1371
      auto at(size_type i) const { return this->get().at(i); }
96,256✔
1372

1373
      auto at(size_type i) { return this->get().at(i); }
443,326,347✔
1374

1375
      auto set(size_type i) { return this->get().set(i); }
1✔
1376

1377
      auto unset(size_type i) { return this->get().unset(i); }
1✔
1378

1379
      auto flip(size_type i) { return this->get().flip(i); }
1✔
1380

1381
      auto flip() { return this->get().flip(); }
1✔
1382

1383
      template <typename OutT>
1384
      auto as() const {
73✔
1385
         return this->get().template as<OutT>();
73✔
1386
      }
1387

1388
      template <bitvectorish OutT = T>
1389
      auto subvector(size_type pos, std::optional<size_type> length = std::nullopt) const {
128,266✔
1390
         return this->get().template subvector<OutT>(pos, length);
128,207✔
1391
      }
1392

1393
      template <typename OutT>
1394
         requires(std::unsigned_integral<strong_type_wrapped_type<OutT>> &&
1395
                  !std::same_as<bool, strong_type_wrapped_type<OutT>>)
1396
      auto subvector(size_type pos) const {
74,940✔
1397
         return this->get().template subvector<OutT>(pos);
74,940✔
1398
      }
1399

1400
      template <typename InT>
1401
         requires(std::unsigned_integral<strong_type_wrapped_type<InT>> && !std::same_as<bool, InT>)
1402
      void subvector_replace(size_type pos, InT value) {
74,939✔
1403
         return this->get().subvector_replace(pos, value);
74,939✔
1404
      }
1405

1406
      template <bitvectorish OtherT>
1407
      auto equals(const OtherT& other) const {
45✔
1408
         return this->get().equals(other);
45✔
1409
      }
1410

1411
      auto push_back(bool b) { return this->get().push_back(b); }
397,281✔
1412

1413
      auto pop_back() { return this->get().pop_back(); }
2✔
1414

1415
      auto front() const { return this->get().front(); }
1416

1417
      auto front() { return this->get().front(); }
1✔
1418

1419
      auto back() const { return this->get().back(); }
1420

1421
      auto back() { return this->get().back(); }
1✔
1422

1423
      auto any_vartime() const { return this->get().any_vartime(); }
1✔
1424

1425
      auto all_vartime() const { return this->get().all_vartime(); }
1✔
1426

1427
      auto none_vartime() const { return this->get().none_vartime(); }
1✔
1428

1429
      auto has_odd_hamming_weight() const { return this->get().has_odd_hamming_weight(); }
1✔
1430

1431
      auto hamming_weight() const { return this->get().hamming_weight(); }
83✔
1432

1433
      auto from_bytes(std::span<const uint8_t> bytes, std::optional<size_type> bits = std::nullopt) {
1434
         return this->get().from_bytes(bytes, bits);
1435
      }
1436

1437
      template <typename OutT = T>
1438
      auto to_bytes() const {
1439
         return this->get().template to_bytes<OutT>();
1440
      }
1441

1442
      auto to_bytes(std::span<uint8_t> out) const { return this->get().to_bytes(out); }
59✔
1443

1444
      auto to_string() const { return this->get().to_string(); }
1✔
1445

1446
      auto capacity() const { return this->get().capacity(); }
1447

1448
      auto reserve(size_type n) { return this->get().reserve(n); }
1449

1450
      constexpr void _const_time_poison() const { this->get()._const_time_poison(); }
71✔
1451

1452
      constexpr void _const_time_unpoison() const { this->get()._const_time_unpoison(); }
98✔
1453
};
1454

1455
}  // namespace detail
1456

1457
}  // namespace Botan
1458

1459
#endif
STATUS · Troubleshooting · Open an Issue · Sales · Support · CAREERS · ENTERPRISE · START FREE · SCHEDULE DEMO
ANNOUNCEMENTS · TWITTER · TOS & SLA · Supported CI Services · What's a CI service? · Automated Testing

© 2026 Coveralls, Inc