Function compare_2d

test/speed_root.cpp:46–58

}

void compare_2d(unsigned n, int distrib) {
  auto r = random_array(n, distrib);

  double best_root = std::numeric_limits<double>::max();
  for (unsigned k = 0; k < 20; ++k) {
    TH2I hroot("", "", 100, 0, 1, 100, 0, 1);
    auto t = clock();
    for (unsigned i = 0; i < n / 2; ++i) hroot.Fill(r[2 * i], r[2 * i + 1]);
    t = clock() - t;
    best_root = std::min(best_root, double(t) / CLOCKS_PER_SEC);
  }
  printf("root %.3f\n", best_root);
}

void compare_3d(unsigned n, int distrib) {
  auto r = random_array(n, distrib);