QGramHashMapReadFinder::QGramHashMapReadFinder(THMMFragStore& fragStore,
					       StatsKeeper* stats_keeper,
					       int k)
    : HashBasedReadFinder(k),
      fragStore(fragStore),
      stats_keeper(stats_keeper),
      readCount(length(fragStore.readSeqStore)),
      retrieved_count(static_cast<int*>(calloc(length(fragStore.readSeqStore), sizeof(int)))),
      max_read_length(0) {

    SEQAN_PROTIMESTART(build_index);
    BuildIndex();
    std::cerr << "Building Index in " << SEQAN_PROTIMEDIFF(build_index) << " seconds." << std::endl;
}
Example #2
0
QGramHashMapReadFinder::QGramHashMapReadFinder(THMMFragStore& fragStore,
					       StatsKeeper& stats_keeper,
					       int k)
    : fragStore(fragStore),
      stats_keeper(stats_keeper),
      rnaseq_k(k),
      readCount(length(fragStore.readSeqStore)),
      retrieved_count(static_cast<int*>(calloc(length(fragStore.readSeqStore), sizeof(int)))),
      extra_bits(-1),
      mask ((1ULL << rnaseq_k * BIT_SHIFT) - 1),
      max_read_length(0),
      grammap(NULL) {

    SEQAN_PROTIMESTART(build_index);
    BuildIndex();
    std::cerr << "Building Index in " << SEQAN_PROTIMEDIFF(build_index) << " seconds." << std::endl;
}