benchmarkHyperRectDomain-catch.cpp

 
#include <iostream>
#include <numeric>
#include <chrono>
 
#include "DGtal/base/Common.h"
#include "DGtal/kernel/SpaceND.h"
#include "DGtal/kernel/PointVector.h"
#include "DGtal/kernel/domains/HyperRectDomain.h"
 
#include "DGtalCatch.h"
 
using namespace DGtal;
using namespace std;
 
// Context for each benchmark
struct BenchDomain
{
  static constexpr std::size_t dim = 3;
  static constexpr signed long long int size = 200;
 
  using Space = DGtal::SpaceND<dim>;
  using Point = Space::Point;
  using Domain = DGtal::HyperRectDomain<Space>;
 
  BenchDomain()
    : a(Point::diagonal(0))
    , b(Point::diagonal(size))
    , domain(Domain(a, b))
    , dimensions(Point::dimension)
    {
      std::iota(dimensions.begin(), dimensions.end(), Dimension(0));
    }
 
  Point a, b;
  Domain domain;
  std::vector<Point::Dimension> dimensions;
};
 
 
TEST_CASE_METHOD(BenchDomain, "Benchmarking HyperRectDomain iterators using Catch2", "[catch]")
{
  BENCHMARK("Domain traversal")
    {
      Point check;
 
      for (auto const& pt : domain)
        check += pt;
 
      return check;
    };
 
  BENCHMARK("Domain reverse traversal")
    {
      Point check;
 
      for (auto it = domain.rbegin(), it_end = domain.rend(); it != it_end; ++it)
        check += *it;
 
      return check;
    };
 
  BENCHMARK("Domain traversal using subRange")
    {
      Point check;
 
      for (auto const& pt : domain.subRange(dimensions))
        check += pt;
 
      return check;
    };
 
  BENCHMARK("Domain reverse traversal using subRange")
    {
      Point check;
      const auto range = domain.subRange(dimensions);
 
      for (auto it = range.rbegin(), it_end = range.rend(); it != it_end; ++it)
        check += *it;
 
      return check;
    };
}