File: ScoreAnchors.hpp

package info (click to toggle)
pbseqlib 5.3.4%2Bdfsg-3
  • links: PTS, VCS
  • area: main
  • in suites: bullseye
  • size: 7,020 kB
  • sloc: cpp: 77,246; python: 331; sh: 103; makefile: 42
file content (40 lines) | stat: -rw-r--r-- 1,642 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
#ifndef _BLASR_SCORE_ANCHORS_HPP_
#define _BLASR_SCORE_ANCHORS_HPP_

#include <cmath>

#include <alignment/statistics/cdfs.hpp>
#include <alignment/statistics/pdfs.hpp>
#include <alignment/tuples/DNATuple.hpp>
#include <alignment/tuples/TupleCountTable.hpp>
#include <alignment/tuples/TupleMetrics.hpp>

template <typename TSequence, typename T_Tuple>
int GetTupleCount(TSequence &seq, DNALength startPos, TupleMetrics &tm,
                  TupleCountTable<TSequence, T_Tuple> &ct, int &count);

template <typename TSequence, typename T_Tuple>
int PMatch(TSequence &seq, DNALength startPos, DNALength length, TupleMetrics &tm,
           TupleCountTable<TSequence, T_Tuple> &ct, float &pMatch);

template <typename TSequence, typename T_Tuple>
int POneOrMoreMatches(TSequence &seq, DNALength startPos, DNALength length, TupleMetrics &tm,
                      TupleCountTable<TSequence, T_Tuple> &ct, float &pValue);

template <typename TSequence, typename T_Tuple>
int SumRightShiftMarginalTupleCounts(TupleMetrics &tm, TupleCountTable<TSequence, T_Tuple> &ct,
                                     T_Tuple curTuple, int nextNuc, int &nextSeqCount);

template <typename TSequence, typename T_Tuple>
int ComputeTotalTupleCount(TupleMetrics &tm, TupleCountTable<TSequence, T_Tuple> &ct,
                           TSequence &seq, int start = 0, int end = -1);

template <typename TSequence, typename T_Tuple>
int ComputeAverageTupleCount(TupleMetrics &tm, TupleCountTable<TSequence, T_Tuple> &ct,
                             TSequence &seq);

inline int ComputeExpectedFirstWaitingTime(float lambda);

#include "ScoreAnchorsImpl.hpp"

#endif