githubmirror/omim

coding/bwt_coder.hpp
125 строк3 KB
Maxim Pimenov
[base] Replaced the namespace my with base.
18 сен 2018, 14:46
18 сен 2018, 14:46412387c
#pragma once

#include "coding/huffman.hpp"
#include "coding/varint.hpp"

#include "base/bwt.hpp"
#include "base/move_to_front.hpp"

#include <algorithm>
#include <cstdint>
#include <iterator>
#include <string>
#include <vector>

namespace coding
{
class BWTCoder
{
public:
  struct Params
  {
    size_t m_blockSize = 32000;
  };

  template <typename Sink>
  static void EncodeAndWriteBlock(Sink & sink, size_t n, uint8_t const * s,
                                  std::vector<uint8_t> & bwtBuffer)
  {
    bwtBuffer.resize(n);
    auto const start = base::BWT(n, s, bwtBuffer.data());

    base::MoveToFront mtf;
    for (auto & b : bwtBuffer)
      b = mtf.Transform(b);

    WriteVarUint(sink, start);

    HuffmanCoder huffman;
    huffman.Init(bwtBuffer.begin(), bwtBuffer.end());
    huffman.WriteEncoding(sink);
    huffman.EncodeAndWrite(sink, bwtBuffer.begin(), bwtBuffer.end());
  }

  template <typename Sink>
  static void EncodeAndWriteBlock(Sink & sink, size_t n, uint8_t const * s)
  {
    std::vector<uint8_t> bwtBuffer;
    EncodeAndWriteBlock(sink, n, s, bwtBuffer);
  }

  template <typename Sink>
  static void EncodeAndWriteBlock(Sink & sink, std::string const & s)
  {
    EncodeAndWriteBlock(sink, s.size(), reinterpret_cast<uint8_t const *>(s.data()));
  }

  template <typename Sink>
  static void EncodeAndWrite(Params const & params, Sink & sink, size_t n, uint8_t const * s)
  {
    CHECK(params.m_blockSize != 0, ());
    CHECK_GREATER(n + params.m_blockSize, n, ());

    std::vector<uint8_t> bwtBuffer;

    size_t const numBlocks = (n + params.m_blockSize - 1) / params.m_blockSize;
    WriteVarUint(sink, numBlocks);
    for (size_t i = 0; i < n; i += params.m_blockSize)
    {
      auto const m = std::min(n - i, params.m_blockSize);
      EncodeAndWriteBlock(sink, m, s + i, bwtBuffer);
    }
  }

  template <typename Source, typename OutIt>
  static OutIt ReadAndDecodeBlock(Source & source, std::vector<uint8_t> & bwtBuffer,
                                  std::vector<uint8_t> & revBuffer, OutIt it)
  {
    auto const start = ReadVarUint<uint64_t, Source>(source);

    HuffmanCoder huffman;
    huffman.ReadEncoding(source);

    bwtBuffer.clear();
    huffman.ReadAndDecode(source, std::back_inserter(bwtBuffer));

    size_t const n = bwtBuffer.size();
    base::MoveToFront mtf;
    for (size_t i = 0; i < n; ++i)
    {
      auto const b = mtf[bwtBuffer[i]];
      bwtBuffer[i] = b;
      mtf.Transform(b);
    }

    if (n != 0)
      CHECK_LESS(start, n, ());

    revBuffer.resize(n);
    base::RevBWT(n, static_cast<size_t>(start), bwtBuffer.data(), revBuffer.data());
    return std::copy(revBuffer.begin(), revBuffer.end(), it);
  }

  template <typename Source, typename OutIt>
  static OutIt ReadAndDecodeBlock(Source & source, OutIt it)
  {
    std::vector<uint8_t> bwtBuffer;
    std::vector<uint8_t> revBuffer;
    return ReadAndDecodeBlock(source, bwtBuffer, revBuffer, it);
  }

  template <typename Source, typename OutIt>
  static OutIt ReadAndDecode(Source & source, OutIt it)
  {
    auto const numBlocks = ReadVarUint<uint64_t, Source>(source);
    CHECK_LESS(numBlocks, std::numeric_limits<size_t>::max(), ());

    std::vector<uint8_t> bwtBuffer;
    std::vector<uint8_t> revBuffer;

    for (size_t i = 0; i < static_cast<size_t>(numBlocks); ++i)
      it = ReadAndDecodeBlock(source, bwtBuffer, revBuffer, it);
    return it;
  }
};
}  // namespace coding