aahash.hpp

#ifndef BTLLIB_AAHASH_HPP
#define BTLLIB_AAHASH_HPP
 
#include <cstdint>
#include <cstring>
#include <limits>
#include <memory>
#include <string>
#include <string_view>
#include <vector>
 
namespace btllib {
 
using SpacedSeed = std::vector<unsigned>;
 
inline std::vector<SpacedSeed>
aa_parse_seeds(const std::vector<std::string>& seeds)
{
  std::vector<SpacedSeed> seed_vec;
  for (const auto& seed : seeds) {
    SpacedSeed seed_vec_tmp;
    for (const auto& c : seed) {
      seed_vec_tmp.push_back((unsigned)(c - '0'));
    }
    seed_vec.push_back(seed_vec_tmp);
  }
  return seed_vec;
}
 
class AAHash;
 
class SeedAAHash;
 
class AAHash
{
  static constexpr const char* HASH_FN_NAME = "aahash1";
 
private:
  friend class SeedAAHash;
 
  bool init();
 
  const char* seq;
  size_t seq_len;
  const uint8_t hash_num;
  const uint16_t k;
  unsigned level;
 
  size_t pos;
  bool initialized = false;
  std::unique_ptr<uint64_t[]> hashes_array;
 
public:
  AAHash(std::string_view seq,
         uint8_t hash_num,
         uint16_t k,
         unsigned level,
         size_t pos = 0)
    : seq(seq.data())
    , seq_len(seq.size())
    , hash_num(hash_num)
    , k(k)
    , level(level)
    , pos(pos)
    , hashes_array(new uint64_t[hash_num])
  {
  }
 
  AAHash(const AAHash& aahash)
    : seq(aahash.seq)
    , seq_len(aahash.seq_len)
    , hash_num(aahash.hash_num)
    , k(aahash.k)
    , level(aahash.level)
    , pos(aahash.pos)
    , initialized(aahash.initialized)
    , hashes_array(new uint64_t[hash_num])
  {
    std::memcpy(hashes_array.get(),
                aahash.hashes_array.get(),
                hash_num * sizeof(uint64_t));
  }
 
  AAHash(AAHash&&) = default;
 
  bool roll();
 
  const uint64_t* hashes() const { return hashes_array.get(); }
  size_t get_pos() const { return pos; }
  unsigned get_hash_num() const { return hash_num; }
  unsigned get_k() const { return k; }
  uint64_t get_forward_hash() const { return hashes_array[0]; }
  unsigned get_level() const { return level; }
  const char* get_seq() const { return seq; }
};
 
class SeedAAHash
{
private:
  AAHash aahash;
  const unsigned hash_num_per_seed;
  std::unique_ptr<uint64_t[]> hashes_array;
  std::vector<SpacedSeed> seeds;
  bool verify_seed();
  void init();
 
public:
  SeedAAHash(const char* seq,
             const std::vector<SpacedSeed>& seeds,
             unsigned hash_num_per_seed,
             unsigned k,
             size_t pos = 0)
    : aahash(seq, 1, k, 1, pos)
    , hash_num_per_seed(hash_num_per_seed)
    , hashes_array(new uint64_t[hash_num_per_seed * seeds.size()])
    , seeds(seeds)
  {
    init();
  }
  SeedAAHash(const std::string& seq,
             const std::vector<SpacedSeed>& seeds,
             unsigned hash_num_per_seed,
             unsigned k,
             size_t pos = 0)
    : aahash(seq, 1, k, 1, pos)
    , hash_num_per_seed(hash_num_per_seed)
    , hashes_array(new uint64_t[hash_num_per_seed * seeds.size()])
    , seeds(seeds)
  {
    init();
  }
  SeedAAHash(const char* seq,
             const std::vector<std::string>& seeds,
             unsigned hash_num_per_seed,
             unsigned k,
             size_t pos = 0)
    : aahash(seq, 1, k, 1, pos)
    , hash_num_per_seed(hash_num_per_seed)
    , hashes_array(new uint64_t[hash_num_per_seed * seeds.size()])
    , seeds(aa_parse_seeds(seeds))
  {
    init();
  }
  SeedAAHash(const std::string& seq,
             const std::vector<std::string>& seeds,
             unsigned hash_num_per_seed,
             unsigned k,
             size_t pos = 0)
    : aahash(seq, 1, k, 1, pos)
    , hash_num_per_seed(hash_num_per_seed)
    , hashes_array(new uint64_t[hash_num_per_seed * seeds.size()])
    , seeds(aa_parse_seeds(seeds))
  {
    init();
  }
 
  SeedAAHash(const SeedAAHash& seed_aahash)
    : aahash(seed_aahash.aahash)
    , hash_num_per_seed(seed_aahash.hash_num_per_seed)
    , hashes_array(new uint64_t[hash_num_per_seed * seed_aahash.seeds.size()])
    , seeds(seed_aahash.seeds)
  {
    std::memcpy(hashes_array.get(),
                seed_aahash.hashes_array.get(),
                hash_num_per_seed * seeds.size() * sizeof(uint64_t));
  }
  SeedAAHash(SeedAAHash&&) = default;
 
  bool roll();
 
  const uint64_t* hashes() const { return hashes_array.get(); }
 
  size_t get_pos() const { return aahash.get_pos(); }
  unsigned get_hash_num() const { return aahash.get_hash_num(); }
  unsigned get_hash_num_per_seed() const { return hash_num_per_seed; }
  unsigned get_k() const { return aahash.get_k(); }
};
 
} // namespace btllib
 
#endif