forked from simongog/sdsl-lite
-
Notifications
You must be signed in to change notification settings - Fork 0
/
lcp_construct_test.cpp
119 lines (107 loc) · 4.39 KB
/
lcp_construct_test.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
107
108
109
110
111
112
113
114
115
116
117
118
119
#include <sdsl/suffix_arrays.hpp>
#include <sdsl/construct_lcp.hpp>
#include <sdsl/construct_bwt.hpp>
#include "gtest/gtest.h"
#include <vector>
#include <string>
#include <map>
using namespace sdsl;
using namespace std;
namespace
{
string test_file, temp_dir, test_id,output_file;
typedef map<string, void (*)(cache_config&)> tMSFP;// map <name, lcp method>
// The fixture for testing class int_vector.
class lcp_construct_test : public ::testing::Test
{
protected:
lcp_construct_test():CHECK_KEY("CHECK_"+string(conf::KEY_LCP)) { }
virtual ~lcp_construct_test() { }
// If the constructor and destructor are not enough for setting up
// and cleaning up each test, you can define the following methods:
virtual void SetUp()
{
test_config = cache_config(false, temp_dir, test_id);
lcp_function["bwt_based"] = &construct_lcp_bwt_based;
lcp_function["bwt_based2"] = &construct_lcp_bwt_based2;
lcp_function["PHI"] = &construct_lcp_PHI<8>;
lcp_function["semi_extern_PHI"] = &construct_lcp_semi_extern_PHI;
lcp_function["go"] = &construct_lcp_go;
lcp_function["goPHI"] = &construct_lcp_goPHI;
uint8_t num_bytes = 1;
{
// Prepare Input
int_vector<8> text;
ASSERT_TRUE(load_vector_from_file(text, test_file, num_bytes));
ASSERT_TRUE(contains_no_zero_symbol(text, test_file));
append_zero_symbol(text);
ASSERT_TRUE(store_to_cache(text, conf::KEY_TEXT, test_config));
// Construct SA
int_vector<> sa(text.size(), 0, bits::hi(text.size())+1);
algorithm::calculate_sa((const unsigned char*)text.data(), text.size(), sa);
ASSERT_TRUE(store_to_cache(sa, conf::KEY_SA, test_config));
}
{
// Construct BWT
construct_bwt<8>(test_config);
}
{
// Construct LCP
construct_lcp_kasai<8>(test_config);
std::rename(cache_file_name(conf::KEY_LCP, test_config).c_str(),
cache_file_name(CHECK_KEY, test_config).c_str());
test_config.file_map.erase(conf::KEY_LCP);
}
}
virtual void TearDown()
{
sdsl::remove(cache_file_name(CHECK_KEY, test_config));
}
cache_config test_config;
tMSFP lcp_function;
string CHECK_KEY;
};
TEST_F(lcp_construct_test, construct_lcp)
{
for (tMSFP::const_iterator it = this->lcp_function.begin(), end = this->lcp_function.end(); it != end; ++it) {
string info = "construct_lcp_" + (it->first) + " on test file " + test_file;
// Construct LCP array
(it->second)(this->test_config);
// Check LCP array
int_vector<> lcp_check, lcp;
string lcp_check_file = cache_file_name(CHECK_KEY, this->test_config);
string lcp_file = cache_file_name(conf::KEY_LCP, this->test_config);
ASSERT_TRUE(load_from_file(lcp_check, lcp_check_file))
<< info << " could not load reference lcp array";
ASSERT_TRUE(load_from_file(lcp, lcp_file))
<< info << " could not load created lcp array";
ASSERT_EQ(lcp_check.size(), lcp.size())
<< info << " lcp array size differ";
for (uint64_t j=0; j<lcp.size(); ++j) {
ASSERT_EQ(lcp_check[j], lcp[j])
<< info << " value differ:" << " lcp_check[" << j << "]="
<< lcp_check[j] << "!=" << lcp[j] << "=lcp["<< j << "]";
}
// Clean up LCP array
sdsl::remove(cache_file_name(conf::KEY_LCP, this->test_config));
}
}
} // namespace
int main(int argc, char** argv)
{
::testing::InitGoogleTest(&argc, argv);
if (argc < 4) {
// LCOV_EXCL_START
cout << "Usage: " << argv[0] << " test_file output_file tmp_dir" << endl;
cout << " (1) Generates the SA, BWT and LCP; arrays are stored in tmp_dir." << endl;
cout << " (2) Generates LCP with other algorithm and checks the result." << endl;
cout << " (3) Deletes all generated files." << endl;
return 1;
// LCOV_EXCL_STOP
}
test_file = argv[1];
output_file = argv[2];
temp_dir = argv[3];
test_id = to_string(util::pid());
return RUN_ALL_TESTS();
}