114 JParser<> zap(
"Example program to test performance of merge sort.");
121 catch(
const exception &error) {
122 FATAL(error.what() << endl);
128 if (rate_Hz <= 0.0) {
129 FATAL(
"Rate " << rate_Hz <<
" < 0 Hz." << endl);
132 const double period = 1.0e9 / rate_Hz;
134 const double Tmin = 0.0;
135 const double Tmax = 1.0e8;
137 const int NUMBER_OF_PMTS = 31;
138 const int NUMBER_OF_MODULES = 115 * 18;
142 typedef double hit_type;
144 const JElement_t getElement;
151 JBuffer2D input(NUMBER_OF_PMTS);
153 size_t number_of_hits = 0;
155 for (
int i = 0; i != NUMBER_OF_PMTS; ++i) {
157 for (
double t1 = Tmin + gRandom->Exp(period); t1 < Tmax; t1 += gRandom->Exp(period), ++number_of_hits) {
158 input[i].push_back(getElement(t1));
161 NOTICE(
"PMT[" << setw(2) << i <<
"] " << input[i].size() << endl);
163 putEndMarker(input[i], JElement_t::getEndMarker());
169 JTimer(
"std::inplace_merge"),
177 for (
int i = 0; i != NUMBER_OF_MODULES; ++i) {
183 for (
int i = 0; i != NUMBER_OF_PMTS; ++i) {
184 copy(input[i].begin(), input[i].end(), back_inserter(buffer));
187 std::sort(buffer.begin(), buffer.end());
191 ASSERT(buffer.size() == number_of_hits,
"Test std::sort.");
192 ASSERT(is_sorted(buffer.begin(), buffer.end()),
"Test std::sort.");
201 for (
int i = 0; i != NUMBER_OF_MODULES; ++i) {
207 for (
int i = 0; i != NUMBER_OF_PMTS; ++i) {
208 n += input[i].size();
213 JBuffer1D::iterator out = buffer.begin();
217 for (
int i = 0; i != NUMBER_OF_PMTS; ++i) {
219 out = copy(input[i].begin(), input[i].end(), out);
221 delimiter[i+1] =
distance(buffer.begin(), out);
224 inplace_merge(buffer.begin(), delimiter.size(), delimiter.data());
228 ASSERT(buffer.size() == number_of_hits,
"Test std::inplace_merge.");
229 ASSERT(is_sorted(buffer.begin(), buffer.end()),
"Test std::inplace_merge.");
237 for (
int i = 0; i != NUMBER_OF_MODULES; ++i) {
243 merge(input, buffer);
247 ASSERT(buffer.size() == number_of_hits,
"Test JMergeSort.");
248 ASSERT(is_sorted(buffer.begin(), buffer.end()),
"Test JMergeSort.");
253 for (
int i = 0; i !=
sizeof(timer)/
sizeof(timer[0]); ++i) {
254 timer[i].
print(cout,
false);
std::vector< T >::difference_type distance(typename std::vector< T >::const_iterator first, typename PhysicsEvent::const_iterator< T > second)
Specialisation of STL distance.