--- /dev/null
+#ifndef _TURBOPFOR_H
+#define _TURBOPFOR_H 1
+
+// A reimplementation of parts of the TurboPFor codecs, using the same
+// storage format. These are not as fast as the reference implementation
+// (about 1/3 of the performance), and do not support the same breadth of
+// codecs (in particular, only delta-plus-1 is implemented, and only 32-bit
+// docids are tested), but aim to be more portable and easier-to-understand.
+// In particular, they will compile on x86 without SSE4.1 or AVX support.
+//
+// The main reference is https://michael.stapelberg.ch/posts/2019-02-05-turbopfor-analysis/,
+// although some implementation details have been worked out by studying the
+// TurboPFor code.
+
+#include <assert.h>
+#include <endian.h>
+#include <stdint.h>
+#include <string.h>
+#include <limits.h>
+
+#include <algorithm>
+
+template<class Docid>
+Docid read_le(const void *in)
+{
+ Docid val;
+ memcpy(&val, in, sizeof(val));
+ if constexpr (sizeof(Docid) == 8) {
+ return le64toh(val);
+ } else if constexpr (sizeof(Docid) == 4) {
+ return le32toh(val);
+ } else if constexpr (sizeof(Docid) == 2) {
+ return le16toh(val);
+ } else if constexpr (sizeof(Docid) == 1) {
+ return val;
+ } else {
+ assert(false);
+ }
+}
+
+// Reads a single value with an encoding that looks a bit like PrefixVarint.
+// It's unclear why this doesn't use the varbyte encoding.
+template<class Docid>
+const unsigned char *read_baseval(const unsigned char *in, Docid *out)
+{
+ //fprintf(stderr, "baseval: 0x%02x 0x%02x 0x%02x 0x%02x\n", in[0], in[1], in[2], in[3]);
+ if (*in < 128) {
+ *out = *in;
+ return in + 1;
+ } else if (*in < 192) {
+ *out = ((uint32_t(in[0]) << 8) | uint32_t(in[1])) & 0x3fff;
+ return in + 2;
+ } else if (*in < 224) {
+ *out = ((uint32_t(in[0]) << 16) |
+ (uint32_t(in[2]) << 8) |
+ (uint32_t(in[1]))) & 0x1fffff;
+ return in + 3;
+ } else {
+ assert(false); // Not implemented.
+ }
+}
+
+template<class Docid>
+const unsigned char *read_vb(const unsigned char *in, Docid *out)
+{
+ if (*in <= 176) {
+ *out = *in;
+ return in + 1;
+ } else if (*in <= 240) {
+ *out = ((uint32_t(in[0] - 177) << 8) | uint32_t(in[1])) + 177;
+ return in + 2;
+ } else if (*in <= 248) {
+ *out = ((uint32_t(in[0] - 241) << 16) | read_le<uint16_t>(in + 1)) + 16561;
+ return in + 3;
+ } else if (*in == 249) {
+ *out = (uint32_t(in[1])) |
+ (uint32_t(in[2]) << 8) |
+ (uint32_t(in[3]) << 16);
+ return in + 4;
+ } else if (*in == 250) {
+ *out = read_le<uint32_t>(in + 1);
+ return in + 5;
+ } else {
+ assert(false);
+ }
+}
+
+struct BitReader {
+public:
+ BitReader(const unsigned char *in, unsigned bits)
+ : in(in), bits(bits), mask((1U << bits) - 1) {}
+ uint32_t read()
+ {
+ uint32_t val = (read_le<uint32_t>(in) >> bits_used) & mask;
+
+ bits_used += bits;
+ in += bits_used / 8;
+ bits_used %= 8;
+
+ return val;
+ }
+
+private:
+ const unsigned char *in;
+ const unsigned bits;
+ const unsigned mask;
+ unsigned bits_used = 0;
+};
+
+// Does not properly account for overflow.
+inline unsigned div_round_up(unsigned val, unsigned div)
+{
+ return (val + div - 1) / div;
+}
+
+inline unsigned bytes_for_packed_bits(unsigned num, unsigned bit_width)
+{
+ return div_round_up(num * bit_width, CHAR_BIT);
+}
+
+// Constant block. Layout:
+//
+// - Bit width (6 bits) | type << 6
+// - Base values (<bits> bits, rounded up to nearest byte)
+template<class Docid>
+const unsigned char *decode_constant(const unsigned char *in, unsigned num, Docid *out)
+{
+ const unsigned bit_width = *in++ & 0x3f;
+ Docid val = read_le<Docid>(in);
+ if (bit_width < sizeof(Docid) * 8) {
+ val &= ((1U << bit_width) - 1);
+ }
+
+ Docid *prev_out = out - 1;
+ for (unsigned i = 0; i < num; ++i) {
+ out[i] = val + prev_out[i] + 1;
+ }
+ return in + div_round_up(bit_width, 8);
+}
+
+// FOR block (ie., PFor without exceptions). Layout:
+//
+// - Bit width (6 bits) | type << 6
+// - Base values (<num> values of <bits> bits, rounded up to a multiple of 32 values)
+template<class Docid>
+const unsigned char *decode_for(const unsigned char *in, unsigned num, Docid *out)
+{
+ const unsigned bit_width = *in++ & 0x3f;
+
+ Docid *prev_out = out - 1;
+ BitReader bs(in, bit_width);
+ for (unsigned i = 0; i < num; ++i) {
+ out[i] = bs.read() + prev_out[i] + 1;
+ }
+ return in + bytes_for_packed_bits(num, bit_width);
+}
+
+// PFor block with bitmap exceptions. Layout:
+//
+// - Bit width (6 bits) | type << 6
+// - Exception bit width (8 bits)
+// - Bitmap of which values have exceptions (<num> bits, rounded up to a byte)
+// - Exceptions (<num_exc> values of <bits_exc> bits, rounded up to a byte)
+// - Base values (<num> values of <bits> bits, rounded up to a byte)
+template<class Docid>
+const unsigned char *decode_pfor_bitmap(const unsigned char *in, unsigned num, Docid *out)
+{
+ memset(out, 0, num * sizeof(Docid));
+
+ const unsigned bit_width = *in++ & 0x3f;
+ const unsigned exception_bit_width = *in++;
+
+ // Decode exceptions.
+ {
+ const uint64_t *exception_bitmap_ptr = reinterpret_cast<const uint64_t *>(in);
+ in += div_round_up(num, 8);
+
+ int num_exceptions = 0;
+
+ BitReader bs(in, exception_bit_width);
+ for (unsigned i = 0; i < num; i += 64, ++exception_bitmap_ptr) {
+ uint64_t exceptions = read_le<uint64_t>(exception_bitmap_ptr);
+ if (num - i < 64) {
+ // We've read some bytes past the end, so clear out the junk bits.
+ exceptions &= (1ULL << (num - i)) - 1;
+ }
+ for (; exceptions != 0; exceptions &= exceptions - 1, ++num_exceptions) {
+ unsigned idx = (ffsll(exceptions) - 1) + i;
+ out[idx] = bs.read() << bit_width;
+ }
+ }
+ in += bytes_for_packed_bits(num_exceptions, exception_bit_width);
+ }
+
+ // Decode the base values, and delta-decode.
+ Docid *prev_out = out - 1;
+ BitReader bs(in, bit_width);
+ for (unsigned i = 0; i < num; ++i) {
+ out[i] = (out[i] | bs.read()) + prev_out[i] + 1;
+ }
+ return in + bytes_for_packed_bits(num, bit_width);
+}
+
+// PFor block with variable-byte exceptions. Layout:
+//
+// - Bit width (6 bits) | type << 6
+// - Number of exceptions (8 bits)
+// - Base values (<num> values of <bits> bits, rounded up to a byte)
+// - Exceptions:
+// - If first byte is 255, <num_exc> 32-bit values (does not include the 255 byte)
+// - Else, <num_exc> varbyte-encoded values (includes the non-255 byte)
+// - Indexes of exceptions (<num_exc> bytes).
+template<unsigned BlockSize, class Docid>
+const unsigned char *decode_pfor_vb(const unsigned char *in, unsigned num, Docid *out)
+{
+ //fprintf(stderr, "in=%p out=%p num=%u\n", in, out, num);
+
+ const unsigned bit_width = *in++ & 0x3f;
+ unsigned num_exceptions = *in++;
+
+ // Decode the base values.
+ BitReader bs(in, bit_width);
+ for (unsigned i = 0; i < num; ++i) {
+ out[i] = bs.read();
+ }
+ in += bytes_for_packed_bits(num, bit_width);
+
+ // Decode exceptions.
+ Docid exceptions[BlockSize];
+ if (*in == 255) {
+ ++in;
+ for (unsigned i = 0; i < num_exceptions; ++i) {
+ exceptions[i] = read_le<Docid>(in);
+ in += sizeof(Docid);
+ }
+ } else {
+ for (unsigned i = 0; i < num_exceptions; ++i) {
+ in = read_vb(in, &exceptions[i]);
+ }
+ }
+ // Apply exceptions.
+ for (unsigned i = 0; i < num_exceptions; ++i) {
+ unsigned idx = *in++;
+ out[idx] |= exceptions[i] << bit_width;
+ }
+
+ // Delta-decode.
+ Docid *prev_out = out - 1;
+ for (unsigned i = 0; i < num; ++i) {
+ out[i] = out[i] + prev_out[i] + 1;
+ }
+
+ return in;
+}
+
+enum BlockType {
+ FOR = 0,
+ PFOR_VB = 1,
+ PFOR_BITMAP = 2,
+ CONSTANT = 3
+};
+
+template<unsigned BlockSize, class Docid>
+const unsigned char *decode_pfor_delta1(const unsigned char *in, unsigned num, Docid *out)
+{
+ if (num == 0) {
+ return in;
+ }
+ in = read_baseval(in, out++);
+
+ for (unsigned i = 1; i < num; i += BlockSize, out += BlockSize) {
+ const unsigned num_this_block = std::min<unsigned>(num - i, BlockSize);
+ switch (in[0] >> 6) {
+ case BlockType::FOR:
+ dprintf("%d+%d: blocktype=%d (for), bitwidth=%d\n", i, num_this_block, in[0] >> 6, in[0] & 0x3f);
+ in = decode_for(in, num_this_block, out);
+ break;
+ case BlockType::PFOR_VB:
+ dprintf("%d+%d: blocktype=%d (pfor + vb), bitwidth=%d\n", i, num_this_block, in[0] >> 6, in[0] & 0x3f);
+ in = decode_pfor_vb<BlockSize>(in, num_this_block, out);
+ break;
+ case BlockType::PFOR_BITMAP:
+ dprintf("%d+%d: blocktype=%d (pfor + bitmap), bitwidth=%d\n", i, num_this_block, in[0] >> 6, in[0] & 0x3f);
+ in = decode_pfor_bitmap(in, num_this_block, out);
+ break;
+ case BlockType::CONSTANT:
+ dprintf("%d+%d: blocktype=%d (constant), bitwidth=%d\n", i, num_this_block, in[0] >> 6, in[0] & 0x3f);
+ in = decode_constant(in, num_this_block, out);
+ break;
+ }
+ }
+
+ return in;
+}
+
+#endif // !defined(_TURBOPFOR_H)