-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathhash.h
73 lines (57 loc) · 1.08 KB
/
hash.h
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
static size_t hash_simple(const char *kmer, size_t k)
{
size_t return_value = 0;
while (k--) {
char c = *kmer++;
size_t val = 0;
switch (c) {
case 'A': val = 0; break;
case 'C': val = 1; break;
case 'G': val = 2; break;
case 'T': val = 3; break;
}
return_value <<= 2;
return_value |= val;
}
return return_value;
}
static size_t hash_twiddle(const char *kmer, size_t k)
{
size_t return_value = 0;
while (k--) {
char c = *kmer++;
c &= 6;
c ^= c >> 1;
return_value <<= 2;
return_value |= c >> 1;
}
return return_value;
}
static size_t hash_twiddle_length(const char *kmer, size_t K)
{
size_t return_value = 0;
size_t k = 0;
while (k < K) {
char c = kmer[k++];
c &= 6;
c ^= c >> 1;
return_value <<= 2;
return_value |= c >> 1;
}
return return_value;
}
static size_t hash_table(const char *kmer, size_t K)
{
static char table[127];
table['A'] = 0;
table['C'] = 1;
table['G'] = 2;
table['T'] = 3;
size_t return_value = 0;
size_t k = 0;
while (k < K) {
return_value <<= 2;
return_value |= table[kmer[k++]];
}
return return_value;
}