File: Pos.hpp

package info (click to toggle)
pbseqlib 5.3.5%2Bdfsg-4
  • links: PTS, VCS
  • area: main
  • in suites: bookworm
  • size: 7,020 kB
  • sloc: cpp: 77,250; python: 331; sh: 103; makefile: 41
file content (43 lines) | stat: -rw-r--r-- 1,121 bytes parent folder | download | duplicates (4)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
#ifndef _BLASR_POS_HPP_
#define _BLASR_POS_HPP_

#include <fstream>
#include <vector>

#include <pbdata/Types.h>
#include <alignment/bwt/PackedHash.hpp>
#include <pbdata/DNASequence.hpp>
#include <pbdata/utils/BitUtils.hpp>

template <typename T_BWT_Sequence>
class Pos
{
public:
    static const unsigned int stride = 8;
    PackedHash packedHash;
    std::vector<int> hashCount;
    std::vector<int> fullPos;
    int hasDebugInformation;
    void Write(std::ostream &out) { packedHash.Write(out); }

    void Read(std::istream &in) { packedHash.Read(in); }

    void InitializeFromSuffixArray(DNALength suffixArray[], DNALength suffixArrayLength)
    {
        DNALength p;
        packedHash.Allocate(suffixArrayLength);
        std::fill(hashCount.begin(), hashCount.end(), 0);
        for (p = 0; p < suffixArrayLength; p++) {
            if (suffixArray[p] % stride == 0) {
                packedHash.AddValue(p, suffixArray[p]);
            }
        }
    }

    int Lookup(DNALength bwtPos, DNALength &seqPos)
    {
        return packedHash.LookupValue(bwtPos - 1, seqPos);
    }
};

#endif  // _BLASR_POS_HPP_