6 #include "misc/utils.h"
13 main(int argc, char** argv)
15 std::cout << "LCP sampler" << std::endl;
18 std::cout << "Usage: sample_plcp base_name sample_rate" << std::endl;
22 std::string base_name = argv[1];
23 std::string lcp_name = base_name + LCP_SAMPLES_EXTENSION;
24 std::cout << "LCP samples: " << lcp_name << std::endl;
25 std::ofstream lcp_file(lcp_name.c_str(), std::ios_base::binary);
28 std::cerr << "Error creating LCP sample file!" << std::endl;
32 usint sample_rate = atoi(argv[2]);
33 std::cout << "Sample rate: " << sample_rate << std::endl;
34 std::cout << std::endl;
35 RLCSA* rlcsa = new RLCSA(base_name);
37 double start = readTimer();
38 pair_type* sampled_values = 0;
39 usint samples = rlcsa->sampleLCP(sample_rate, sampled_values, true);
40 usint data_size = rlcsa->getSize();
41 delete rlcsa; // Saves memory.
42 LCPSamples lcp(sampled_values, data_size, samples, true, true);
43 lcp.writeTo(lcp_file);
45 double time = readTimer() - start;
46 double megabytes = data_size / (double)MEGABYTE;
47 double size = lcp.reportSize() / (double)MEGABYTE;
48 std::cout << megabytes << " megabytes in " << time << " seconds (" << (megabytes / time) << " MB/s)" << std::endl;
49 std::cout << "Sampled LCP size: " << size << " MB" << std::endl;
50 std::cout << std::endl;