Logo ROOT  

Reference Guide
histspeedtest.cxx
Go to the documentation of this file.
1/// \file histspeedtest.cxx
2///
3/// \warning This is part of the ROOT 7 prototype! It will change without notice. It might trigger earthquakes. Feedback
4/// is welcome!
5/// \author Axel Naumann <axel@cern.ch>
6
7#include "TRandom3.h"
8#include <vector>
9#include <chrono>
10#include <iostream>
11
12#include "TH1.h"
13#include "TH2.h"
14#include "TH3.h"
15
16#include "ROOT/RHist.hxx"
17#include "ROOT/RHistBufferedFill.hxx"
18
19using namespace ROOT;
20using namespace std;
21
22constexpr unsigned short gRepeat = 2;
23
24/* DataTypes:
25
26 THistDataContent
27 THistDataUncertainty
28 THistDataMomentUncert
29 THistDataRuntime
30and
31 MyTHistDataNoStat
32 MyTHistDataContent
33 MyTHistDataMomentUncert
34
35 /opt/build/root_builds/rootcling.cmake/include/ROOT/THistBinIter.h:53:50: error: no member named 'GetUncertainty' in
36'ROOT::Experimental::THistDataContent<2, double, ROOT::Experimental::THistDataDefaultStorage>::TBinStat<double>' auto
37GetUncertainty() const { return GetStat().GetUncertainty(); }
38 ~~~~~~~~~ ^
39
40 new ones (STATCLASSES)
41
42 THistStatContent
43 THistStatUncertainty
44 THistStatTotalSumOfWeights
45 THistStatTotalSumOfSquaredWeights
46 THistDataMomentUncert
47
48
49 */
50
51#ifndef STATCLASSES
52#define STATCLASSES Experimental::RHistStatContent, Experimental::RHistStatUncertainty
53#endif
54
55struct Timer {
56 using TimePoint_t = decltype(std::chrono::high_resolution_clock::now());
57
58 const char *fTitle;
59 size_t fCount;
60 TimePoint_t fStart;
61
62 Timer(const char *title, size_t count)
63 : fTitle(title), fCount(count), fStart(std::chrono::high_resolution_clock::now())
64 {}
65
66 ~Timer()
67 {
68 using namespace std::chrono;
69 auto end = high_resolution_clock::now();
70 duration<double> time_span = duration_cast<duration<double>>(end - fStart);
71 // std::cout << fCount << " * " << fTitle << ": " << time_span.count() << " seconds \n";
72 // std::cout << fCount << " * " << fTitle << ": " << fCount / (1e6) / time_span.count()
73 // << " millions per seconds\n";
74 std::cout << fCount << " * " << fTitle << ": " << time_span.count() << " seconds, \t";
75 std::cout << fCount / (1e6) / time_span.count() << " millions per seconds \n";
76 }
77};
78
79constexpr UInt_t gStride = 32; // TH1::GetDefaultBufferSize()
80
81struct BinEdges {
82 static constexpr size_t fNBinsX = 4;
83 static constexpr size_t fNBinsY = 5;
84 double fXBins[4];
85 double fYBins[5];
86
87 BinEdges(double minValue, double maxValue)
88 {
89 if (maxValue < minValue)
90 swap(minValue, maxValue);
91 double range = maxValue - minValue;
92
93 double x[fNBinsX] = {0., 0.1, 0.3, 1.};
94 double y[fNBinsY] = {0., 0.1, 0.2, 0.3, 1.};
95
96 for (size_t i = 0; i < fNBinsX; ++i)
97 fXBins[i] = minValue + range * x[i];
98 for (size_t i = 0; i < fNBinsY; ++i)
99 fYBins[i] = minValue + range * y[i];
100 }
101
102 using AConf_t = Experimental::RAxisConfig;
103
104 AConf_t GetConfigX() const { return AConf_t(std::span<const double>(fXBins).to_vector()); }
105 AConf_t GetConfigY() const { return AConf_t(std::span<const double>(fYBins).to_vector()); }
106};
107
108template <typename T>
109void GenerateInput(std::vector<T> &numbers, double minVal, double maxVal, UInt_t seed)
110{
111 Timer t("GenerateInput", numbers.size());
112 if (minVal > maxVal) {
113 std::swap(minVal, maxVal);
114 }
115 T range = maxVal - minVal;
116 TRandom3 r(seed);
117 size_t len = numbers.size();
118 for (auto c = numbers.begin(); c != numbers.end(); ++c) {
119 *c = minVal + range * r.Rndm();
120 }
121}
122
123std::string
124MakeTitle(std::string_view version, std::string_view histname, std::string_view title, std::string_view axis)
125{
126 std::string result =
127 std::string(version) + " " + std::string(histname) + " " + std::string(title) + " [" + std::string(axis) + "]";
128 return result;
129}
130
131template <int dim, typename type>
132const char *GetHist();
133
134template <>
135const char *GetHist<2, double>()
136{
137 return "2D";
138};
139template <>
140const char *GetHist<2, float>()
141{
142 return "2F";
143};
144
145template <>
146const char *GetHist<1, double>()
147{
148 return "1D";
149};
150template <>
151const char *GetHist<1, float>()
152{
153 return "1F";
154};
155
156namespace R7 {
157const char *gVersion = "R7";
158
159template <typename T, unsigned short kNDim>
160struct Dim;
161
162template <typename T>
163struct Dim<T, 2> {
164
165 constexpr static unsigned short kNDim = 2;
166 using ExpTH2 = Experimental::RHist<kNDim, T, STATCLASSES>;
167
168 using FillFunc_t = std::add_pointer_t<long(ExpTH2 &hist, std::vector<double> &input, std::string_view type)>;
169
170 struct EE {
171 static constexpr const char *const gType = "regular bin size ";
172
173 template <FillFunc_t filler>
174 static long Execute(std::vector<double> &input, double minVal, double maxVal)
175 {
176
177 ExpTH2 hist({100, minVal, maxVal}, {5, minVal, maxVal});
178 return filler(hist, input, gType);
179 }
180 };
181
182 struct II {
183 static constexpr const char *const gType = "irregular bin size";
184
185 template <FillFunc_t filler>
186 static long Execute(std::vector<double> &input, double minVal, double maxVal)
187 {
188
189 BinEdges edges(minVal, maxVal);
190 ExpTH2 hist(edges.GetConfigX(), edges.GetConfigY());
191 return filler(hist, input, gType);
192 }
193 };
194
195 inline static long fillN(ExpTH2 &hist, std::vector<double> &input, std::string_view gType)
196 {
197
198 using array_t = Experimental::Hist::RCoordArray<2>;
199 array_t *values = (array_t *)(&input[0]);
200 constexpr size_t stride = gStride;
201
202 std::string title = MakeTitle(gVersion, GetHist<kNDim, T>(), "fills N (stride 32)", gType);
203 {
204 Timer t(title.c_str(), input.size() / 2);
205 for (size_t i = 0; i < (input.size() - (stride * 2 - 1)); i += (stride * 2), values += 32) {
206 std::span<array_t> coords(values, 32);
207 hist.FillN(coords);
208 }
209 }
210 return hist.GetNDim();
211 }
212
213 inline static long fillBuffered(ExpTH2 &hist, std::vector<double> &input, std::string_view gType)
214 {
215 Experimental::RHistBufferedFill<ExpTH2> filler(hist);
216 std::string title = MakeTitle(gVersion, GetHist<kNDim, T>(), "fills (buffered) ", gType);
217 {
218 Timer t(title.c_str(), input.size() / 2);
219 for (size_t i = 0; i < input.size() - 1; i += 2)
220 filler.Fill({input[i], input[i + 1]});
221 }
222 return hist.GetNDim();
223 }
224
225 inline static long fill(ExpTH2 &hist, std::vector<double> &input, std::string_view gType)
226 {
227 std::string title = MakeTitle(gVersion, GetHist<kNDim, T>(), "fills ", gType);
228 {
229 Timer t(title.c_str(), input.size() / 2);
230 for (size_t i = 0; i < input.size() - 1; i += 2)
231 hist.Fill({input[i], input[i + 1]});
232 }
233 return hist.GetNDim();
234 }
235}; // DimD
236
237template <typename T>
238struct Dim<T, 1> {
239
240 constexpr static unsigned short kNDim = 1;
241 using ExpTH1 = Experimental::RHist<kNDim, T, STATCLASSES>;
242
243 using FillFunc_t = std::add_pointer_t<long(ExpTH1 &hist, std::vector<double> &input, std::string_view type)>;
244
245 struct EE {
246 static constexpr const char *const gType = "regular bin size ";
247
248 template <FillFunc_t filler>
249 static long Execute(std::vector<double> &input, double minVal, double maxVal)
250 {
251 long result = 0;
252 for (unsigned short i = 0; i < gRepeat; ++i) {
253 ExpTH1 hist({100, minVal, maxVal});
254 result += filler(hist, input, gType);
255 }
256 return result;
257 }
258 };
259
260 struct II {
261 static constexpr const char *const gType = "irregular bin size";
262
263 template <FillFunc_t filler>
264 static long Execute(std::vector<double> &input, double minVal, double maxVal)
265 {
266
267 long result = 0;
268 BinEdges edges(minVal, maxVal);
269 for (unsigned short i = 0; i < gRepeat; ++i) {
270 ExpTH1 hist(edges.GetConfigX());
271 result += filler(hist, input, gType);
272 }
273 return result;
274 }
275 };
276
277 inline static long fillN(ExpTH1 &hist, std::vector<double> &input, std::string_view gType)
278 {
279
280 using array_t = Experimental::Hist::RCoordArray<1>;
281 array_t *values = (array_t *)(&input[0]);
282 constexpr size_t stride = gStride;
283
284 std::string title = MakeTitle(gVersion, GetHist<kNDim, T>(), "fills N (stride 32)", gType);
285 {
286 Timer t(title.c_str(), input.size());
287 for (size_t i = 0; i < (input.size() - (stride - 1)); i += (stride), values += 32) {
288 std::span<array_t> coords(values, 32);
289 hist.FillN(coords);
290 }
291 }
292 return hist.GetNDim();
293 }
294
295 inline static long fillBuffered(ExpTH1 &hist, std::vector<double> &input, std::string_view gType)
296 {
297 Experimental::RHistBufferedFill<ExpTH1> filler(hist);
298 std::string title = MakeTitle(gVersion, GetHist<kNDim, T>(), "fills (buffered) ", gType);
299 {
300 Timer t(title.c_str(), input.size());
301 for (size_t i = 0; i < input.size(); ++i)
302 filler.Fill({input[i]});
303 }
304 return hist.GetNDim();
305 }
306
307 inline static long fill(ExpTH1 &hist, std::vector<double> &input, std::string_view gType)
308 {
309 std::string title = MakeTitle(gVersion, GetHist<kNDim, T>(), "fills ", gType);
310 {
311 Timer t(title.c_str(), input.size());
312 for (size_t i = 0; i < input.size(); ++i)
313 hist.Fill({input[i]});
314 }
315 return hist.GetNDim();
316 }
317}; // Dim1
318
319} // namespace R7
320
321namespace R6 {
322const char *gVersion = "R6";
323
324template <int ndim, typename T>
325struct Redirect;
326template <>
327struct Redirect<2, float> {
328 using HistType_t = TH2F;
329};
330template <>
331struct Redirect<2, double> {
332 using HistType_t = TH2D;
333};
334template <>
335struct Redirect<1, float> {
336 using HistType_t = TH1F;
337};
338template <>
339struct Redirect<1, double> {
340 using HistType_t = TH1D;
341};
342
343template <typename T, int kNDim>
344struct Dim;
345
346template <typename T>
347struct Dim<T, 2> {
348
349 constexpr static unsigned short kNDim = 2;
350 using HistType_t = typename Redirect<kNDim, T>::HistType_t;
351
352 using FillFunc_t = std::add_pointer_t<long(HistType_t &hist, std::vector<double> &input, std::string_view type)>;
353
354 struct EE {
355
356 static constexpr const char *const gType = "regular bin size ";
357
358 template <FillFunc_t filler>
359 static long Execute(std::vector<double> &input, double minVal, double maxVal)
360 {
361
362 long result = 0;
363 for (unsigned short i = 0; i < gRepeat; ++i) {
364 HistType_t hist("a", "a hist", 100, minVal, maxVal, 5, minVal, maxVal);
365 result += filler(hist, input, gType);
366 }
367 return result;
368 }
369 };
370
371 struct II {
372
373 static constexpr const char *const gType = "irregular bin size";
374
375 template <FillFunc_t filler>
376 static long Execute(std::vector<double> &input, double minVal, double maxVal)
377 {
378 long result = 0;
379 BinEdges edges(minVal, maxVal);
380 for (unsigned short i = 0; i < gRepeat; ++i) {
381 HistType_t hist("a", "a hist", edges.fNBinsX - 1, edges.fXBins, edges.fNBinsY - 1, edges.fYBins);
382 result += filler(hist, input, gType);
383 }
384 return result;
385 }
386 };
387
388 static long fillBuffered(HistType_t &hist, std::vector<double> &input, std::string_view gType)
389 {
390 std::string title = MakeTitle(gVersion, GetHist<kNDim, T>(), "fills (buffered) ", gType);
391 hist.SetBuffer(TH1::GetDefaultBufferSize());
392 {
393 // Timer t("R6 2D fills [regular bins size]",input.size()/2);
394 Timer t(title.c_str(), input.size() / 2);
395 for (size_t i = 0; i < input.size() - 1; i += 2)
396 hist.Fill(input[i], input[i + 1]);
397 }
398 return (long)hist.GetEntries();
399 }
400
401 static long fillN(HistType_t &hist, std::vector<double> &input, std::string_view gType)
402 {
403 std::string title = MakeTitle(gVersion, GetHist<kNDim, T>(), "fills N (stride 32)", gType);
404 constexpr size_t stride = gStride;
405 {
406 // Timer t("R6 2D fills [regular bins size]",input.size()/2);
407 Timer t(title.c_str(), input.size() / 2);
408 for (size_t i = 0; i < (input.size() - (stride * 2 - 1)); i += (stride * 2))
409 hist.FillN(gStride, &(input[i]), &(input[i + gStride]), nullptr);
410 }
411 return (long)hist.GetEntries();
412 }
413
414 static long fill(HistType_t &hist, std::vector<double> &input, std::string_view gType)
415 {
416 std::string title = MakeTitle(gVersion, GetHist<kNDim, T>(), "fills ", gType);
417 {
418 // Timer t("R6 2D fills [regular bins size]",input.size()/2);
419 Timer t(title.c_str(), input.size() / 2);
420 for (size_t i = 0; i < input.size() - 1; i += 2)
421 hist.Fill(input[i], input[i + 1]);
422 }
423 return (long)hist.GetEntries();
424 }
425}; // Dim
426
427template <typename T>
428struct Dim<T, 1> {
429
430 constexpr static unsigned short kNDim = 1;
431 using HistType_t = typename Redirect<kNDim, T>::HistType_t;
432
433 using FillFunc_t = std::add_pointer_t<long(HistType_t &hist, std::vector<double> &input, std::string_view type)>;
434
435 struct EE {
436
437 static constexpr const char *const gType = "regular bin size ";
438
439 template <FillFunc_t filler>
440 static long Execute(std::vector<double> &input, double minVal, double maxVal)
441 {
442
443 long result = 0;
444 for (unsigned short i = 0; i < gRepeat; ++i) {
445 HistType_t hist("a", "a hist", 100, minVal, maxVal);
446 result += filler(hist, input, gType);
447 }
448 return result;
449 }
450 };
451
452 struct II {
453
454 static constexpr const char *const gType = "irregular bin size";
455
456 template <FillFunc_t filler>
457 static long Execute(std::vector<double> &input, double minVal, double maxVal)
458 {
459 long result = 0;
460 BinEdges edges(minVal, maxVal);
461 for (unsigned short i = 0; i < gRepeat; ++i) {
462 HistType_t hist("a", "a hist", edges.fNBinsX - 1, edges.fXBins);
463 result += filler(hist, input, gType);
464 }
465 return result;
466 }
467 };
468
469 static long fillBuffered(HistType_t &hist, std::vector<double> &input, std::string_view gType)
470 {
471 std::string title = MakeTitle(gVersion, GetHist<kNDim, T>(), "fills (buffered) ", gType);
472 hist.SetBuffer(TH1::GetDefaultBufferSize());
473 {
474 // Timer t("R6 2D fills [regular bins size]",input.size()/2);
475 Timer t(title.c_str(), input.size());
476 for (size_t i = 0; i < input.size() - 1; ++i)
477 hist.Fill(input[i]);
478 }
479 return (long)hist.GetEntries();
480 }
481
482 static long fillN(HistType_t &hist, std::vector<double> &input, std::string_view gType)
483 {
484 std::string title = MakeTitle(gVersion, GetHist<kNDim, T>(), "fills N (stride 32)", gType);
485 constexpr size_t stride = gStride;
486 {
487 // Timer t("R6 2D fills [regular bins size]",input.size()/2);
488 Timer t(title.c_str(), input.size());
489 for (size_t i = 0; i < (input.size() - (stride - 1)); i += (stride))
490 hist.FillN(gStride, &(input[i]), nullptr);
491 }
492 return (long)hist.GetEntries();
493 }
494
495 static long fill(HistType_t &hist, std::vector<double> &input, std::string_view gType)
496 {
497 std::string title = MakeTitle(gVersion, GetHist<kNDim, T>(), "fills ", gType);
498 {
499 // Timer t("R6 2D fills [regular bins size]",input.size()/2);
500 Timer t(title.c_str(), input.size());
501 for (size_t i = 0; i < input.size(); ++i)
502 hist.Fill(input[i]);
503 }
504 return (long)hist.GetEntries();
505 }
506}; // Dim1
507} // namespace R6
508
509template <typename T, unsigned short kNDim>
510void speedtest(size_t count = (size_t)(1e6));
511
512template <>
513void speedtest<double, 2>(size_t count)
514{
515 using DataType_t = double;
516 static constexpr unsigned short kNDim = 2;
517
518 TH1::AddDirectory(kFALSE);
519
520 std::vector<double> input; // (count);
521 input.resize(count);
522
523 double minVal = -5.0;
524 double maxVal = +5.0;
525 GenerateInput(input, minVal, maxVal, 0);
526
527 // Make sure we have some overflow.
528 minVal *= 0.9;
529 maxVal *= 0.9;
530
531 cout << '\n';
532
533 R7::Dim<DataType_t, kNDim>::EE::Execute<R7::Dim<DataType_t, kNDim>::fillBuffered>(input, minVal, maxVal);
534 R7::Dim<DataType_t, kNDim>::EE::Execute<R7::Dim<DataType_t, kNDim>::fillBuffered>(input, minVal, maxVal);
535 R7::Dim<DataType_t, kNDim>::EE::Execute<R7::Dim<DataType_t, kNDim>::fillN>(input, minVal, maxVal);
536 R7::Dim<DataType_t, kNDim>::EE::Execute<R7::Dim<DataType_t, kNDim>::fillN>(input, minVal, maxVal);
537 R7::Dim<DataType_t, kNDim>::EE::Execute<R7::Dim<DataType_t, kNDim>::fill>(input, minVal, maxVal);
538 R7::Dim<DataType_t, kNDim>::EE::Execute<R7::Dim<DataType_t, kNDim>::fill>(input, minVal, maxVal);
539
540 cout << '\n';
541
542 R7::Dim<DataType_t, kNDim>::II::Execute<R7::Dim<DataType_t, kNDim>::fillBuffered>(input, minVal, maxVal);
543 R7::Dim<DataType_t, kNDim>::II::Execute<R7::Dim<DataType_t, kNDim>::fillBuffered>(input, minVal, maxVal);
544 R7::Dim<DataType_t, kNDim>::II::Execute<R7::Dim<DataType_t, kNDim>::fillN>(input, minVal, maxVal);
545 R7::Dim<DataType_t, kNDim>::II::Execute<R7::Dim<DataType_t, kNDim>::fillN>(input, minVal, maxVal);
546 R7::Dim<DataType_t, kNDim>::II::Execute<R7::Dim<DataType_t, kNDim>::fill>(input, minVal, maxVal);
547 R7::Dim<DataType_t, kNDim>::II::Execute<R7::Dim<DataType_t, kNDim>::fill>(input, minVal, maxVal);
548
549 cout << '\n';
550
551 R6::Dim<DataType_t, kNDim>::EE::Execute<R6::Dim<DataType_t, kNDim>::fillBuffered>(input, minVal, maxVal);
552 R6::Dim<DataType_t, kNDim>::EE::Execute<R6::Dim<DataType_t, kNDim>::fillBuffered>(input, minVal, maxVal);
553 R6::Dim<DataType_t, kNDim>::EE::Execute<R6::Dim<DataType_t, kNDim>::fillN>(input, minVal, maxVal);
554 R6::Dim<DataType_t, kNDim>::EE::Execute<R6::Dim<DataType_t, kNDim>::fillN>(input, minVal, maxVal);
555 R6::Dim<DataType_t, kNDim>::EE::Execute<R6::Dim<DataType_t, kNDim>::fill>(input, minVal, maxVal);
556 R6::Dim<DataType_t, kNDim>::EE::Execute<R6::Dim<DataType_t, kNDim>::fill>(input, minVal, maxVal);
557
558 cout << '\n';
559
560 R6::Dim<DataType_t, kNDim>::II::Execute<R6::Dim<DataType_t, kNDim>::fillBuffered>(input, minVal, maxVal);
561 R6::Dim<DataType_t, kNDim>::II::Execute<R6::Dim<DataType_t, kNDim>::fillBuffered>(input, minVal, maxVal);
562 R6::Dim<DataType_t, kNDim>::II::Execute<R6::Dim<DataType_t, kNDim>::fillN>(input, minVal, maxVal);
563 R6::Dim<DataType_t, kNDim>::II::Execute<R6::Dim<DataType_t, kNDim>::fillN>(input, minVal, maxVal);
564 R6::Dim<DataType_t, kNDim>::II::Execute<R6::Dim<DataType_t, kNDim>::fill>(input, minVal, maxVal);
565 R6::Dim<DataType_t, kNDim>::II::Execute<R6::Dim<DataType_t, kNDim>::fill>(input, minVal, maxVal);
566
567 cout << '\n';
568}
569
570// these are copy/paste to work around failure to properly instantiate the template :(
571
572template <>
573void speedtest<float, 2>(size_t count)
574{
575 using DataType_t = float;
576 constexpr unsigned short kNDim = 2;
577
578 TH1::AddDirectory(kFALSE);
579
580 std::vector<double> input; // (count);
581 input.resize(count);
582
583 double minVal = -5.0;
584 double maxVal = +5.0;
585 GenerateInput(input, minVal, maxVal, 0);
586
587 // Make sure we have some overflow.
588 minVal *= 0.9;
589 maxVal *= 0.9;
590
591 cout << '\n';
592
593 R7::Dim<DataType_t, kNDim>::EE::Execute<R7::Dim<DataType_t, kNDim>::fillBuffered>(input, minVal, maxVal);
594 R7::Dim<DataType_t, kNDim>::EE::Execute<R7::Dim<DataType_t, kNDim>::fillBuffered>(input, minVal, maxVal);
595 R7::Dim<DataType_t, kNDim>::EE::Execute<R7::Dim<DataType_t, kNDim>::fillN>(input, minVal, maxVal);
596 R7::Dim<DataType_t, kNDim>::EE::Execute<R7::Dim<DataType_t, kNDim>::fillN>(input, minVal, maxVal);
597 R7::Dim<DataType_t, kNDim>::EE::Execute<R7::Dim<DataType_t, kNDim>::fill>(input, minVal, maxVal);
598 R7::Dim<DataType_t, kNDim>::EE::Execute<R7::Dim<DataType_t, kNDim>::fill>(input, minVal, maxVal);
599
600 cout << '\n';
601
602 R7::Dim<DataType_t, kNDim>::II::Execute<R7::Dim<DataType_t, kNDim>::fillBuffered>(input, minVal, maxVal);
603 R7::Dim<DataType_t, kNDim>::II::Execute<R7::Dim<DataType_t, kNDim>::fillBuffered>(input, minVal, maxVal);
604 R7::Dim<DataType_t, kNDim>::II::Execute<R7::Dim<DataType_t, kNDim>::fillN>(input, minVal, maxVal);
605 R7::Dim<DataType_t, kNDim>::II::Execute<R7::Dim<DataType_t, kNDim>::fillN>(input, minVal, maxVal);
606 R7::Dim<DataType_t, kNDim>::II::Execute<R7::Dim<DataType_t, kNDim>::fill>(input, minVal, maxVal);
607 R7::Dim<DataType_t, kNDim>::II::Execute<R7::Dim<DataType_t, kNDim>::fill>(input, minVal, maxVal);
608
609 cout << '\n';
610
611 R6::Dim<DataType_t, kNDim>::EE::Execute<R6::Dim<DataType_t, kNDim>::fillBuffered>(input, minVal, maxVal);
612 R6::Dim<DataType_t, kNDim>::EE::Execute<R6::Dim<DataType_t, kNDim>::fillBuffered>(input, minVal, maxVal);
613 R6::Dim<DataType_t, kNDim>::EE::Execute<R6::Dim<DataType_t, kNDim>::fillN>(input, minVal, maxVal);
614 R6::Dim<DataType_t, kNDim>::EE::Execute<R6::Dim<DataType_t, kNDim>::fillN>(input, minVal, maxVal);
615 R6::Dim<DataType_t, kNDim>::EE::Execute<R6::Dim<DataType_t, kNDim>::fill>(input, minVal, maxVal);
616 R6::Dim<DataType_t, kNDim>::EE::Execute<R6::Dim<DataType_t, kNDim>::fill>(input, minVal, maxVal);
617
618 cout << '\n';
619
620 R6::Dim<DataType_t, kNDim>::II::Execute<R6::Dim<DataType_t, kNDim>::fillBuffered>(input, minVal, maxVal);
621 R6::Dim<DataType_t, kNDim>::II::Execute<R6::Dim<DataType_t, kNDim>::fillBuffered>(input, minVal, maxVal);
622 R6::Dim<DataType_t, kNDim>::II::Execute<R6::Dim<DataType_t, kNDim>::fillN>(input, minVal, maxVal);
623 R6::Dim<DataType_t, kNDim>::II::Execute<R6::Dim<DataType_t, kNDim>::fillN>(input, minVal, maxVal);
624 R6::Dim<DataType_t, kNDim>::II::Execute<R6::Dim<DataType_t, kNDim>::fill>(input, minVal, maxVal);
625 R6::Dim<DataType_t, kNDim>::II::Execute<R6::Dim<DataType_t, kNDim>::fill>(input, minVal, maxVal);
626}
627
628template <>
629void speedtest<double, 1>(size_t count)
630{
631 using DataType_t = double;
632 static constexpr unsigned short kNDim = 1;
633
634 TH1::AddDirectory(kFALSE);
635
636 std::vector<double> input; // (count);
637 input.resize(count);
638
639 double minVal = -5.0;
640 double maxVal = +5.0;
641 GenerateInput(input, minVal, maxVal, 0);
642
643 // Make sure we have some overflow.
644 minVal *= 0.9;
645 maxVal *= 0.9;
646
647 cout << '\n';
648
649 R7::Dim<DataType_t, kNDim>::EE::Execute<R7::Dim<DataType_t, kNDim>::fillBuffered>(input, minVal, maxVal);
650 R7::Dim<DataType_t, kNDim>::EE::Execute<R7::Dim<DataType_t, kNDim>::fillBuffered>(input, minVal, maxVal);
651 R7::Dim<DataType_t, kNDim>::EE::Execute<R7::Dim<DataType_t, kNDim>::fillN>(input, minVal, maxVal);
652 R7::Dim<DataType_t, kNDim>::EE::Execute<R7::Dim<DataType_t, kNDim>::fillN>(input, minVal, maxVal);
653 R7::Dim<DataType_t, kNDim>::EE::Execute<R7::Dim<DataType_t, kNDim>::fill>(input, minVal, maxVal);
654 R7::Dim<DataType_t, kNDim>::EE::Execute<R7::Dim<DataType_t, kNDim>::fill>(input, minVal, maxVal);
655
656 cout << '\n';
657
658 R7::Dim<DataType_t, kNDim>::II::Execute<R7::Dim<DataType_t, kNDim>::fillBuffered>(input, minVal, maxVal);
659 R7::Dim<DataType_t, kNDim>::II::Execute<R7::Dim<DataType_t, kNDim>::fillBuffered>(input, minVal, maxVal);
660 R7::Dim<DataType_t, kNDim>::II::Execute<R7::Dim<DataType_t, kNDim>::fillN>(input, minVal, maxVal);
661 R7::Dim<DataType_t, kNDim>::II::Execute<R7::Dim<DataType_t, kNDim>::fillN>(input, minVal, maxVal);
662 R7::Dim<DataType_t, kNDim>::II::Execute<R7::Dim<DataType_t, kNDim>::fill>(input, minVal, maxVal);
663 R7::Dim<DataType_t, kNDim>::II::Execute<R7::Dim<DataType_t, kNDim>::fill>(input, minVal, maxVal);
664
665 cout << '\n';
666
667 R6::Dim<DataType_t, kNDim>::EE::Execute<R6::Dim<DataType_t, kNDim>::fillBuffered>(input, minVal, maxVal);
668 R6::Dim<DataType_t, kNDim>::EE::Execute<R6::Dim<DataType_t, kNDim>::fillBuffered>(input, minVal, maxVal);
669 R6::Dim<DataType_t, kNDim>::EE::Execute<R6::Dim<DataType_t, kNDim>::fillN>(input, minVal, maxVal);
670 R6::Dim<DataType_t, kNDim>::EE::Execute<R6::Dim<DataType_t, kNDim>::fillN>(input, minVal, maxVal);
671 R6::Dim<DataType_t, kNDim>::EE::Execute<R6::Dim<DataType_t, kNDim>::fill>(input, minVal, maxVal);
672 R6::Dim<DataType_t, kNDim>::EE::Execute<R6::Dim<DataType_t, kNDim>::fill>(input, minVal, maxVal);
673
674 cout << '\n';
675
676 R6::Dim<DataType_t, kNDim>::II::Execute<R6::Dim<DataType_t, kNDim>::fillBuffered>(input, minVal, maxVal);
677 R6::Dim<DataType_t, kNDim>::II::Execute<R6::Dim<DataType_t, kNDim>::fillBuffered>(input, minVal, maxVal);
678 R6::Dim<DataType_t, kNDim>::II::Execute<R6::Dim<DataType_t, kNDim>::fillN>(input, minVal, maxVal);
679 R6::Dim<DataType_t, kNDim>::II::Execute<R6::Dim<DataType_t, kNDim>::fillN>(input, minVal, maxVal);
680 R6::Dim<DataType_t, kNDim>::II::Execute<R6::Dim<DataType_t, kNDim>::fill>(input, minVal, maxVal);
681 R6::Dim<DataType_t, kNDim>::II::Execute<R6::Dim<DataType_t, kNDim>::fill>(input, minVal, maxVal);
682}
683
684template <>
685void speedtest<float, 1>(size_t count)
686{
687 using DataType_t = float;
688 static constexpr unsigned short kNDim = 1;
689
690 TH1::AddDirectory(kFALSE);
691
692 std::vector<double> input; // (count);
693 input.resize(count);
694
695 double minVal = -5.0;
696 double maxVal = +5.0;
697 GenerateInput(input, minVal, maxVal, 0);
698
699 // Make sure we have some overflow.
700 minVal *= 0.9;
701 maxVal *= 0.9;
702
703 cout << '\n';
704
705 R7::Dim<DataType_t, kNDim>::EE::Execute<R7::Dim<DataType_t, kNDim>::fillBuffered>(input, minVal, maxVal);
706 R7::Dim<DataType_t, kNDim>::EE::Execute<R7::Dim<DataType_t, kNDim>::fillBuffered>(input, minVal, maxVal);
707 R7::Dim<DataType_t, kNDim>::EE::Execute<R7::Dim<DataType_t, kNDim>::fillN>(input, minVal, maxVal);
708 R7::Dim<DataType_t, kNDim>::EE::Execute<R7::Dim<DataType_t, kNDim>::fillN>(input, minVal, maxVal);
709 R7::Dim<DataType_t, kNDim>::EE::Execute<R7::Dim<DataType_t, kNDim>::fill>(input, minVal, maxVal);
710 R7::Dim<DataType_t, kNDim>::EE::Execute<R7::Dim<DataType_t, kNDim>::fill>(input, minVal, maxVal);
711
712 cout << '\n';
713
714 R7::Dim<DataType_t, kNDim>::II::Execute<R7::Dim<DataType_t, kNDim>::fillBuffered>(input, minVal, maxVal);
715 R7::Dim<DataType_t, kNDim>::II::Execute<R7::Dim<DataType_t, kNDim>::fillBuffered>(input, minVal, maxVal);
716 R7::Dim<DataType_t, kNDim>::II::Execute<R7::Dim<DataType_t, kNDim>::fillN>(input, minVal, maxVal);
717 R7::Dim<DataType_t, kNDim>::II::Execute<R7::Dim<DataType_t, kNDim>::fillN>(input, minVal, maxVal);
718 R7::Dim<DataType_t, kNDim>::II::Execute<R7::Dim<DataType_t, kNDim>::fill>(input, minVal, maxVal);
719 R7::Dim<DataType_t, kNDim>::II::Execute<R7::Dim<DataType_t, kNDim>::fill>(input, minVal, maxVal);
720
721 cout << '\n';
722
723 R6::Dim<DataType_t, kNDim>::EE::Execute<R6::Dim<DataType_t, kNDim>::fillBuffered>(input, minVal, maxVal);
724 R6::Dim<DataType_t, kNDim>::EE::Execute<R6::Dim<DataType_t, kNDim>::fillBuffered>(input, minVal, maxVal);
725 R6::Dim<DataType_t, kNDim>::EE::Execute<R6::Dim<DataType_t, kNDim>::fillN>(input, minVal, maxVal);
726 R6::Dim<DataType_t, kNDim>::EE::Execute<R6::Dim<DataType_t, kNDim>::fillN>(input, minVal, maxVal);
727 R6::Dim<DataType_t, kNDim>::EE::Execute<R6::Dim<DataType_t, kNDim>::fill>(input, minVal, maxVal);
728 R6::Dim<DataType_t, kNDim>::EE::Execute<R6::Dim<DataType_t, kNDim>::fill>(input, minVal, maxVal);
729
730 cout << '\n';
731
732 R6::Dim<DataType_t, kNDim>::II::Execute<R6::Dim<DataType_t, kNDim>::fillBuffered>(input, minVal, maxVal);
733 R6::Dim<DataType_t, kNDim>::II::Execute<R6::Dim<DataType_t, kNDim>::fillBuffered>(input, minVal, maxVal);
734 R6::Dim<DataType_t, kNDim>::II::Execute<R6::Dim<DataType_t, kNDim>::fillN>(input, minVal, maxVal);
735 R6::Dim<DataType_t, kNDim>::II::Execute<R6::Dim<DataType_t, kNDim>::fillN>(input, minVal, maxVal);
736 R6::Dim<DataType_t, kNDim>::II::Execute<R6::Dim<DataType_t, kNDim>::fill>(input, minVal, maxVal);
737 R6::Dim<DataType_t, kNDim>::II::Execute<R6::Dim<DataType_t, kNDim>::fill>(input, minVal, maxVal);
738}
739
740void histspeedtest(size_t iter = 1e6, int what = 255)
741{
742 if (what & 1)
743 speedtest<double, 2>(iter);
744 if (what & 2)
745 speedtest<float, 2>(iter);
746 if (what & 4)
747 speedtest<double, 1>(iter);
748 if (what & 8)
749 speedtest<float, 1>(iter);
750}
751
752int main(int argc, char **argv)
753{
754
755 size_t iter = 1e7;
756 int what = 1 | 2 | 4 | 8;
757 if (argc > 1)
758 iter = atof(argv[1]);
759 if (argc > 2)
760 what = atoi(argv[2]);
761
762 histspeedtest(iter, what);
763}
r
ROOT::R::TRInterface & r
Definition: Object.C:4
c
#define c(i)
Definition: RSha256.hxx:101
UInt_t
unsigned int UInt_t
Definition: RtypesCore.h:42
kFALSE
const Bool_t kFALSE
Definition: RtypesCore.h:88
type
int type
Definition: TGX11.cxx:120
ROOT::Experimental::RAxisConfig
Objects used to configure the different axis types.
Definition: RAxis.hxx:300
ROOT::Experimental::RHist
Histogram class for histograms with DIMENSIONS dimensions, where each bin count is stored by a value ...
Definition: RHist.hxx:51
TH1D
1-D histogram with a double per channel (see TH1 documentation)}
Definition: TH1.h:614
TH1F
1-D histogram with a float per channel (see TH1 documentation)}
Definition: TH1.h:571
TH1::GetDefaultBufferSize
static Int_t GetDefaultBufferSize()
Static function return the default buffer size for automatic histograms the parameter fgBufferSize ma...
Definition: TH1.cxx:4277
TH1::AddDirectory
static void AddDirectory(Bool_t add=kTRUE)
Sets the flag controlling the automatic add of histograms in memory.
Definition: TH1.cxx:1226
TH2D
2-D histogram with a double per channel (see TH1 documentation)}
Definition: TH2.h:292
TH2F
2-D histogram with a float per channel (see TH1 documentation)}
Definition: TH2.h:251
TRandom3
Random number generator class based on M.
Definition: TRandom3.h:27
GetHist< 1, double >
const char * GetHist< 1, double >()
Definition: histspeedtest.cxx:146
gStride
constexpr UInt_t gStride
Definition: histspeedtest.cxx:79
speedtest
void speedtest(size_t count=(size_t)(1e6))
main
int main(int argc, char **argv)
Definition: histspeedtest.cxx:752
GetHist< 2, double >
const char * GetHist< 2, double >()
Definition: histspeedtest.cxx:135
speedtest< double, 1 >
void speedtest< double, 1 >(size_t count)
Definition: histspeedtest.cxx:629
gRepeat
constexpr unsigned short gRepeat
Definition: histspeedtest.cxx:22
speedtest< float, 1 >
void speedtest< float, 1 >(size_t count)
Definition: histspeedtest.cxx:685
GetHist< 2, float >
const char * GetHist< 2, float >()
Definition: histspeedtest.cxx:140
GetHist< 1, float >
const char * GetHist< 1, float >()
Definition: histspeedtest.cxx:151
histspeedtest
void histspeedtest(size_t iter=1e6, int what=255)
Definition: histspeedtest.cxx:740
GenerateInput
void GenerateInput(std::vector< T > &numbers, double minVal, double maxVal, UInt_t seed)
Definition: histspeedtest.cxx:109
GetHist
const char * GetHist()
MakeTitle
std::string MakeTitle(std::string_view version, std::string_view histname, std::string_view title, std::string_view axis)
Definition: histspeedtest.cxx:124
speedtest< double, 2 >
void speedtest< double, 2 >(size_t count)
Definition: histspeedtest.cxx:513
speedtest< float, 2 >
void speedtest< float, 2 >(size_t count)
Definition: histspeedtest.cxx:573
y
Double_t y[n]
Definition: legend1.C:17
x
Double_t x[n]
Definition: legend1.C:17
string_view
basic_string_view< char > string_view
Definition: libcpp_string_view.h:786
R6::gVersion
const char * gVersion
Definition: histspeedtest.cxx:322
R7::gVersion
const char * gVersion
Definition: histspeedtest.cxx:157
ROOT::Experimental::Internal::swap
void swap(RDirectoryEntry &e1, RDirectoryEntry &e2) noexcept
Definition: RDirectoryEntry.hxx:94
ROOT::Math::Chebyshev::T
double T(double x)
Definition: ChebyshevPol.h:34
ROOT
VSD Structures.
Definition: StringConv.hxx:21
RooFit::Timer
RooCmdArg Timer(Bool_t flag=kTRUE)
Definition: RooGlobalFunc.cxx:187
fit1_py.fill
fill
Definition: fit1_py.py:6
gType
static int gType
Definition: proofexecv.cxx:43