Marat Dukhan | 4a4a7fa | 2019-10-21 13:46:14 -0700 | [diff] [blame] | 1 | #include <algorithm> |
| 2 | #include <cfloat> |
| 3 | #include <cmath> |
| 4 | #include <functional> |
| 5 | #include <random> |
| 6 | #include <vector> |
| 7 | |
| 8 | #include "bench/utils.h" |
| 9 | #include <xnnpack/common.h> |
| 10 | #include <xnnpack/params.h> |
| 11 | #include <xnnpack/raddexpminusmax.h> |
| 12 | #include <xnnpack/raddstoreexpminusmax.h> |
| 13 | #include <xnnpack/rmax.h> |
Marat Dukhan | 05ac8e3 | 2019-10-21 15:39:33 -0700 | [diff] [blame^] | 14 | #include <xnnpack/vscale.h> |
Marat Dukhan | 4a4a7fa | 2019-10-21 13:46:14 -0700 | [diff] [blame] | 15 | #include <xnnpack/vscaleexpminusmax.h> |
| 16 | |
| 17 | #include <benchmark/benchmark.h> |
| 18 | |
| 19 | |
| 20 | static void ThreePassSoftargmaxWithRecomputing( |
| 21 | benchmark::State& state, |
| 22 | xnn_f32_rmax_ukernel_function rmax, |
| 23 | xnn_f32_raddexpminusmax_ukernel_function raddexpminusmax, |
| 24 | xnn_f32_vscaleexpminusmax_ukernel_function vscaleexpminusmax) |
| 25 | { |
| 26 | const size_t n = state.range(0); |
| 27 | const size_t cache_line_size_max = 128; |
| 28 | const size_t packed_n = benchmark::utils::roundUp(n, cache_line_size_max / sizeof(float)); |
| 29 | |
| 30 | std::random_device random_device; |
| 31 | auto rng = std::mt19937(random_device()); |
| 32 | auto f32rng = std::bind(std::uniform_real_distribution<float>(-1000.0f, 1000.0f), rng); |
| 33 | |
| 34 | const size_t num_buffers = 1 + |
| 35 | benchmark::utils::divideRoundUp<size_t>(benchmark::utils::GetMaxCacheSize(), packed_n * sizeof(float)); |
| 36 | std::vector<float> x(n); |
| 37 | std::vector<float> y(packed_n * num_buffers); |
| 38 | |
| 39 | std::generate(x.begin(), x.end(), std::ref(f32rng)); |
| 40 | |
| 41 | benchmark::utils::DisableDenormals(); |
| 42 | |
| 43 | size_t buffer_index = 0; |
| 44 | for (auto _ : state) { |
| 45 | benchmark::utils::prefetchToL1(x.data(), x.size() * sizeof(float)); |
| 46 | if (++buffer_index == num_buffers) { |
| 47 | buffer_index = 0; |
| 48 | } |
| 49 | |
| 50 | const auto start = std::chrono::high_resolution_clock::now(); |
| 51 | float x_max = nanf(""); |
| 52 | rmax(n * sizeof(float), x.data(), &x_max); |
| 53 | float y_sum = nanf(""); |
| 54 | raddexpminusmax(n * sizeof(float), x.data(), &y_sum, x_max); |
| 55 | vscaleexpminusmax(n * sizeof(float), x.data(), y.data() + packed_n * buffer_index, x_max, 1.0f / y_sum); |
| 56 | const auto end = std::chrono::high_resolution_clock::now(); |
| 57 | |
| 58 | const auto elapsed_seconds = |
| 59 | std::chrono::duration_cast<std::chrono::duration<double>>(end - start); |
| 60 | state.SetIterationTime(elapsed_seconds.count()); |
| 61 | } |
| 62 | |
| 63 | state.SetItemsProcessed(uint64_t(state.iterations()) * n); |
| 64 | state.SetBytesProcessed(uint64_t(state.iterations()) * 2 * sizeof(float) * n); |
| 65 | } |
| 66 | |
Marat Dukhan | 05ac8e3 | 2019-10-21 15:39:33 -0700 | [diff] [blame^] | 67 | static void ThreePassSoftargmaxWithReloading( |
| 68 | benchmark::State& state, |
| 69 | xnn_f32_rmax_ukernel_function rmax, |
| 70 | xnn_f32_raddstoreexpminusmax_ukernel_function raddstoreexpminusmax, |
| 71 | xnn_f32_vscale_ukernel_function vscale) |
| 72 | { |
| 73 | const size_t n = state.range(0); |
| 74 | const size_t cache_line_size_max = 128; |
| 75 | const size_t packed_n = benchmark::utils::roundUp(n, cache_line_size_max / sizeof(float)); |
| 76 | |
| 77 | std::random_device random_device; |
| 78 | auto rng = std::mt19937(random_device()); |
| 79 | auto f32rng = std::bind(std::uniform_real_distribution<float>(-1000.0f, 1000.0f), rng); |
| 80 | |
| 81 | const size_t num_buffers = 1 + |
| 82 | benchmark::utils::divideRoundUp<size_t>(benchmark::utils::GetMaxCacheSize(), packed_n * sizeof(float)); |
| 83 | std::vector<float> x(n); |
| 84 | std::vector<float> y(packed_n * num_buffers); |
| 85 | |
| 86 | std::generate(x.begin(), x.end(), std::ref(f32rng)); |
| 87 | |
| 88 | benchmark::utils::DisableDenormals(); |
| 89 | |
| 90 | size_t buffer_index = 0; |
| 91 | for (auto _ : state) { |
| 92 | benchmark::utils::prefetchToL1(x.data(), x.size() * sizeof(float)); |
| 93 | if (++buffer_index == num_buffers) { |
| 94 | buffer_index = 0; |
| 95 | } |
| 96 | |
| 97 | const auto start = std::chrono::high_resolution_clock::now(); |
| 98 | float x_max = nanf(""); |
| 99 | rmax(n * sizeof(float), x.data(), &x_max); |
| 100 | float y_sum = nanf(""); |
| 101 | raddstoreexpminusmax(n * sizeof(float), x.data(), y.data() + packed_n * buffer_index, &y_sum, x_max); |
| 102 | vscale(n * sizeof(float), y.data() + packed_n * buffer_index, y.data() + packed_n * buffer_index, 1.0f / y_sum); |
| 103 | const auto end = std::chrono::high_resolution_clock::now(); |
| 104 | |
| 105 | const auto elapsed_seconds = |
| 106 | std::chrono::duration_cast<std::chrono::duration<double>>(end - start); |
| 107 | state.SetIterationTime(elapsed_seconds.count()); |
| 108 | } |
| 109 | |
| 110 | state.SetItemsProcessed(uint64_t(state.iterations()) * n); |
| 111 | state.SetBytesProcessed(uint64_t(state.iterations()) * 2 * sizeof(float) * n); |
| 112 | } |
| 113 | |
Marat Dukhan | 4a4a7fa | 2019-10-21 13:46:14 -0700 | [diff] [blame] | 114 | static void CharacteristicArguments(benchmark::internal::Benchmark* b) { |
| 115 | for (int32_t n = 1000; n <= 10000000; n *= 10) { |
| 116 | b->Arg(n); |
| 117 | b->Arg(3 * n); |
| 118 | } |
| 119 | } |
| 120 | |
| 121 | #if XNN_ARCH_X86 || XNN_ARCH_X86_64 |
| 122 | BENCHMARK_CAPTURE(ThreePassSoftargmaxWithRecomputing, avx512f_p5_scalef_unroll128, |
| 123 | xnn_f32_rmax_ukernel__avx512f, xnn_f32_raddexpminusmax_ukernel__avx512f_p5_scalef_unroll128, xnn_f32_vscaleexpminusmax_ukernel__avx512f_p5_scalef_unroll128) |
| 124 | ->Apply(CharacteristicArguments)->UseManualTime(); |
Marat Dukhan | 05ac8e3 | 2019-10-21 15:39:33 -0700 | [diff] [blame^] | 125 | BENCHMARK_CAPTURE(ThreePassSoftargmaxWithReloading, avx512f_p5_scalef_unroll128, |
| 126 | xnn_f32_rmax_ukernel__avx512f, xnn_f32_raddstoreexpminusmax_ukernel__avx512f_p5_scalef_unroll128, xnn_f32_vscale_ukernel__avx512f_unroll64) |
| 127 | ->Apply(CharacteristicArguments)->UseManualTime(); |
Marat Dukhan | 4a4a7fa | 2019-10-21 13:46:14 -0700 | [diff] [blame] | 128 | |
| 129 | BENCHMARK_CAPTURE(ThreePassSoftargmaxWithRecomputing, avx2_p5_unroll64, |
| 130 | xnn_f32_rmax_ukernel__avx, xnn_f32_raddexpminusmax_ukernel__avx2_p5_unroll64, xnn_f32_vscaleexpminusmax_ukernel__avx2_p5_unroll64) |
| 131 | ->Apply(CharacteristicArguments)->UseManualTime(); |
Marat Dukhan | 05ac8e3 | 2019-10-21 15:39:33 -0700 | [diff] [blame^] | 132 | BENCHMARK_CAPTURE(ThreePassSoftargmaxWithReloading, avx2_p5_unroll64, |
| 133 | xnn_f32_rmax_ukernel__avx, xnn_f32_raddstoreexpminusmax_ukernel__avx2_p5_unroll64, xnn_f32_vscale_ukernel__avx_unroll32) |
| 134 | ->Apply(CharacteristicArguments)->UseManualTime(); |
Marat Dukhan | 4a4a7fa | 2019-10-21 13:46:14 -0700 | [diff] [blame] | 135 | #endif // XNN_ARCH_X86 || XNN_ARCH_X86_64 |
| 136 | |
| 137 | #ifndef XNNPACK_BENCHMARK_NO_MAIN |
| 138 | BENCHMARK_MAIN(); |
| 139 | #endif |