Bench for exps

Author: ProExpertProg

.idea/misc.xml

@@ -25,6 +25,9 @@ target_link_libraries(bench-naive src-lib benchmark::benchmark)
 add_executable(bench-hl-filter hl-filter.cpp)
 target_link_libraries(bench-hl-filter src-lib benchmark::benchmark_main)
 
+add_executable(bench-exps exps.cpp)
+target_link_libraries(bench-exps src-lib benchmark::benchmark_main)
+
 if (ENABLE_CILK)
     add_executable(bench-fftw-cilk fftw-cilk.cpp)
     target_link_libraries(bench-fftw-cilk fftw-cpp benchmark::benchmark fftw3_threads)

cpp/bench/CMakeLists.txt

@@ -0,0 +1,31 @@
+#include "Exponentials.hpp"
+
+#include "benchmark-util.hpp"
+#include <iostream>
+
+using namespace ahr;
+template <typename Exp> static void BM_ExpKXKY(benchmark::State &state) {
+  Dim const X = state.range(0);
+  Dim const Y = state.range(1);
+
+  Grid grid{1, X, Y};
+  Real dt = 1.0;
+  HyperCoefficients hyper = HyperCoefficients::calculate(dt, grid);
+
+  Exp exp{grid};
+  exp.update(hyper, dt);
+
+  for (auto _ : state) {
+    grid.for_each_kxky([&](Dim kx, Dim ky) { benchmark::DoNotOptimize(exp(kx, ky)); });
+  }
+}
+
+BENCHMARK_WMIN(BM_ExpKXKY<exp::Eta>)
+    ->ArgsProduct({{2048, 4096, 8192}, {2048, 4096, 8192}})
+    ->Unit(benchmark::kMillisecond);
+BENCHMARK_WMIN(BM_ExpKXKY<exp::Nu>)
+    ->ArgsProduct({{2048, 4096, 8192}, {2048, 4096, 8192}})
+    ->Unit(benchmark::kMillisecond);
+BENCHMARK_WMIN(BM_ExpKXKY<exp::NuG>)
+    ->ArgsProduct({{2048, 4096, 8192}, {2048, 4096, 8192}})
+    ->Unit(benchmark::kMillisecond);