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
|
#pragma once
#include "data_factory.hpp"
#include "interval_index_iface.hpp"
#include "../coding/var_serial_vector.hpp"
#include "../base/assert.hpp"
#include "../base/base.hpp"
#include "../base/macros.hpp"
#include "../base/stl_add.hpp"
#include "../std/algorithm.hpp"
class ScaleIndexBase
{
public:
enum { NUM_BUCKETS = 18 };
ScaleIndexBase()
{
#ifdef DEBUG
for (size_t i = 0; i < ARRAY_SIZE(kScaleBuckets); ++i)
{
ASSERT_LESS(kScaleBuckets[i], static_cast<uint32_t>(NUM_BUCKETS), (i));
ASSERT(i == 0 || kScaleBuckets[i] >= kScaleBuckets[i-1],
(i, kScaleBuckets[i-1], kScaleBuckets[i]));
}
#endif
}
static uint32_t BucketByScale(uint32_t scale)
{
ASSERT_LESS(scale, ARRAY_SIZE(kScaleBuckets), ());
return scale >= ARRAY_SIZE(kScaleBuckets) ? NUM_BUCKETS - 1 : kScaleBuckets[scale];
}
static pair<uint32_t, uint32_t> ScaleRangeForBucket(uint32_t bucket)
{
// TODO: Cache ScaleRangeForBucket in class member?
ASSERT_LESS(bucket, static_cast<uint32_t>(NUM_BUCKETS), ());
pair<uint32_t, uint32_t> res(ARRAY_SIZE(kScaleBuckets), 0);
for (uint32_t i = 0; i < ARRAY_SIZE(kScaleBuckets); ++i)
{
if (kScaleBuckets[i] == bucket)
{
res.first = min(res.first, i);
res.second = max(res.second, i + 1);
}
}
return res;
}
private:
static uint32_t const kScaleBuckets[18];
};
template <class ReaderT>
class ScaleIndex : public ScaleIndexBase
{
public:
typedef ReaderT ReaderType;
ScaleIndex() {}
explicit ScaleIndex(ReaderT const & reader, IndexFactory & factory)
{
Attach(reader, factory);
}
~ScaleIndex()
{
Clear();
}
void Clear()
{
for_each(m_IndexForScale.begin(), m_IndexForScale.end(), DeleteFunctor());
m_IndexForScale.clear();
}
void Attach(ReaderT const & reader, IndexFactory & factory)
{
Clear();
ReaderSource<ReaderT> source(reader);
VarSerialVectorReader<ReaderT> treesReader(source);
for (size_t i = 0; i < treesReader.Size(); ++i)
m_IndexForScale.push_back(factory.CreateIndex(treesReader.SubReader(i)));
}
template <typename F>
void ForEachInIntervalAndScale(F const & f, uint64_t beg, uint64_t end, uint32_t scale) const
{
size_t const scaleBucket = BucketByScale(scale);
if (scaleBucket < m_IndexForScale.size())
{
IntervalIndexIFace::FunctionT f1(f);
for (size_t i = 0; i <= scaleBucket; ++i)
m_IndexForScale[i]->DoForEach(f1, beg, end);
}
}
private:
vector<IntervalIndexIFace *> m_IndexForScale;
};
|