-
Notifications
You must be signed in to change notification settings - Fork 3
/
KmerIterator.cpp
executable file
·106 lines (82 loc) · 1.97 KB
/
KmerIterator.cpp
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
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
//
// Created by Nolan Woods on 11/6/2018.
//
#ifndef KMER_ITERATOR_CPP
#define KMER_ITERATOR_CPP
#include "KmerIterator.h"
template<typename T>
KmerIterator<T> get_begin(T *str, int k) {
KmerIterator<T> iter(str, k);
return iter;
}
template<typename T>
KmerIterator<T> get_end(T *str, int length, int k) {
KmerIterator<T> iter(str, k);
return iter + (length - k);
}
template<typename T>
encode_type kmer_encode(KmerIterator<T> start, KmerIterator<T> end) {
encode_type encoding = 0;
for (; start != end; start++) {
encoding = encoding << 2;
encoding = start | encoding;
}
return encoding;
}
template<typename T>
KmerIterator<T>::KmerIterator(const T *str, int k) {
ptr = str;
this->k = k;
}
template<typename T>
KmerIterator<T> KmerIterator<T>::operator+=(int val) {
KmerIterator i = *this;
ptr += val;
return i;
}
template<typename T>
KmerIterator<T> &KmerIterator<T>::operator++() {
ptr += 1;
return *this;
}
template<typename T>
KmerIterator<T> KmerIterator<T>::operator++(int val) {
KmerIterator i = *this;
val = 1;
ptr += val;
return i;
}
template<typename T>
encode_type KmerIterator<T>::operator*() {
KmerIterator i = *this;
return kmer_encode<T>(i, i + k);
}
template<typename T>
encode_type KmerIterator<T>::operator|(encode_type val) {
return ordValue(*ptr) | val;
}
template<typename T>
T KmerIterator<T>::operator->() {
return ptr;
}
template<typename T>
bool KmerIterator<T>::operator==(const KmerIterator &rhs) {
return ptr == rhs.ptr;
}
template<typename T>
bool KmerIterator<T>::operator!=(const KmerIterator &rhs) {
return ptr != rhs.ptr;
}
template<typename T>
KmerIterator<T> KmerIterator<T>::operator+(int integer_value) {
KmerIterator i = *this;
i.ptr += integer_value;
return i;
}
template<typename T>
KmerIterator<T> KmerIterator<T>::operator-(int integer_value) {
KmerIterator i = *this;
i.ptr -= integer_value;
return i;
}
#endif