diff --git a/benchmarks/BM_divide.cpp b/benchmarks/BM_divide.cpp
new file mode 100644
index 00000000..824fa120
--- /dev/null
+++ b/benchmarks/BM_divide.cpp
@@ -0,0 +1,96 @@
+// Copyright (c) 2019, Paul Ferrand
+// All rights reserved.
+
+// Redistribution and use in source and binary forms, with or without
+// modification, are permitted provided that the following conditions are met:
+
+// 1. Redistributions of source code must retain the above copyright notice, this
+//    list of conditions and the following disclaimer.
+// 2. Redistributions in binary form must reproduce the above copyright notice,
+//    this list of conditions and the following disclaimer in the documentation
+//    and/or other materials provided with the distribution.
+
+// THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND
+// ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED
+// WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE
+// DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR
+// ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES
+// (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
+// LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND
+// ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
+// (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
+// SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
+
+#include <benchmark/benchmark.h>
+#include <random>
+#include <numeric>
+#include <vector>
+#include <cmath>
+#include <iostream>
+#include "SIMDHelpers.h"
+
+
+class Divide : public benchmark::Fixture {
+public:
+  void SetUp(const ::benchmark::State& state) {
+    std::random_device rd { };
+    std::mt19937 gen { rd() };
+    std::uniform_real_distribution<float> dist { 0, 1 };
+    input = std::vector<float>(state.range(0));
+    output = std::vector<float>(state.range(0));
+    divisor = std::vector<float>(state.range(0));
+    std::generate(divisor.begin(), divisor.end(), [&]() { return dist(gen); });
+    std::generate(input.begin(), input.end(), [&]() { return dist(gen); });
+  }
+
+  void TearDown(const ::benchmark::State& state [[maybe_unused]]) {
+
+  }
+
+  std::vector<float> divisor;
+  std::vector<float> input;
+  std::vector<float> output;
+};
+
+BENCHMARK_DEFINE_F(Divide, Straight)(benchmark::State& state) {
+    for (auto _ : state)
+    {
+        for (int i = 0; i < state.range(0); ++i)
+            output[i] = input[i] / divisor[i];
+    }
+}
+
+BENCHMARK_DEFINE_F(Divide, Scalar)(benchmark::State& state) {
+    for (auto _ : state)
+    {
+        sfz::divide<float, false>(input, divisor, absl::MakeSpan(output));
+    }
+}
+
+BENCHMARK_DEFINE_F(Divide, SIMD)(benchmark::State& state) {
+    for (auto _ : state)
+    {
+        sfz::divide<float, true>(input, divisor, absl::MakeSpan(output));
+    }
+}
+
+BENCHMARK_DEFINE_F(Divide, Scalar_Unaligned)(benchmark::State& state) {
+    for (auto _ : state)
+    {
+        sfz::divide<float, false>(absl::MakeSpan(input).subspan(1), absl::MakeSpan(divisor).subspan(1), absl::MakeSpan(output).subspan(1));
+    }
+}
+
+BENCHMARK_DEFINE_F(Divide, SIMD_Unaligned)(benchmark::State& state) {
+    for (auto _ : state)
+    {
+        sfz::divide<float, true>(absl::MakeSpan(input).subspan(1), absl::MakeSpan(divisor).subspan(1), absl::MakeSpan(output).subspan(1));
+    }
+}
+
+BENCHMARK_REGISTER_F(Divide, Straight)->RangeMultiplier(4)->Range(1 << 2, 1 << 12);
+BENCHMARK_REGISTER_F(Divide, Scalar)->RangeMultiplier(4)->Range(1 << 2, 1 << 12);
+BENCHMARK_REGISTER_F(Divide, SIMD)->RangeMultiplier(4)->Range(1 << 2, 1 << 12);
+BENCHMARK_REGISTER_F(Divide, Scalar_Unaligned)->RangeMultiplier(4)->Range(1 << 2, 1 << 12);
+BENCHMARK_REGISTER_F(Divide, SIMD_Unaligned)->RangeMultiplier(4)->Range(1 << 2, 1 << 12);
+BENCHMARK_MAIN();
diff --git a/benchmarks/CMakeLists.txt b/benchmarks/CMakeLists.txt
index cb8f105d..0294a83c 100644
--- a/benchmarks/CMakeLists.txt
+++ b/benchmarks/CMakeLists.txt
@@ -45,6 +45,10 @@ add_executable(bm_gain BM_gain.cpp ${BENCHMARK_SIMD_SOURCES})
 target_link_libraries(bm_gain benchmark absl::span)
 target_include_directories(bm_gain PRIVATE ../src/sfizz ../src/external)
 
+add_executable(bm_divide BM_divide.cpp ${BENCHMARK_SIMD_SOURCES})
+target_link_libraries(bm_divide benchmark absl::span)
+target_include_directories(bm_divide PRIVATE ../src/sfizz ../src/external)
+
 add_executable(bm_looping BM_looping.cpp ${BENCHMARK_SIMD_SOURCES})
 target_link_libraries(bm_looping benchmark absl::span absl::algorithm)
 target_include_directories(bm_looping PRIVATE ../src/sfizz ../src/external)
@@ -123,6 +127,7 @@ add_dependencies(sfizz_benchmarks
 	bm_interpolationCast
 	bm_mathfuns
 	bm_gain
+	bm_divide
 	bm_looping
 	bm_saturating
 	bm_ramp