AK: Add `AK::SIMD::exp_approximate`

This approximation tries to generate values within 0.1% of their actual
expected value. Microbenchmarks indicate that this iterative SIMD
version can be up to 60x faster than `AK::SIMD::exp`.

2 files changed, 34 insertions, 0 deletions

diff --git a/Tests/AK/CMakeLists.txt b/Tests/AK/CMakeLists.txt
index 77a288d2ea..e26a38e41f 100644
--- a/Tests/AK/CMakeLists.txt
+++ b/Tests/AK/CMakeLists.txt
@@ -61,6 +61,7 @@ set(AK_TEST_SOURCES
     TestQuickSort.cpp
     TestRedBlackTree.cpp
     TestRefPtr.cpp
+    TestSIMD.cpp
     TestSinglyLinkedList.cpp
     TestSourceGenerator.cpp
     TestSourceLocation.cpp
diff --git a/Tests/AK/TestSIMD.cpp b/Tests/AK/TestSIMD.cpp
new file mode 100644
index 0000000000..251904cbf2
--- /dev/null
+++ b/Tests/AK/TestSIMD.cpp
@@ -0,0 +1,33 @@
+/*
+ * Copyright (c) 2023, Jelle Raaijmakers <jelle@gmta.nl>
+ *
+ * SPDX-License-Identifier: BSD-2-Clause
+ */
+
+#include <LibTest/TestCase.h>
+
+#include <AK/SIMD.h>
+#include <AK/SIMDMath.h>
+
+TEST_CASE(exp)
+{
+    AK::SIMD::f32x4 v = { .2f, .4f, .6f, .8f };
+    auto result = AK::SIMD::exp(v);
+
+    EXPECT_APPROXIMATE(result[0], 1.22140276f);
+    EXPECT_APPROXIMATE(result[1], 1.49182470f);
+    EXPECT_APPROXIMATE(result[2], 1.82211880f);
+    EXPECT_APPROXIMATE(result[3], 2.22554093f);
+}
+
+TEST_CASE(exp_approximate)
+{
+    AK::SIMD::f32x4 v = { .2f, .4f, .6f, .8f };
+    auto result = AK::SIMD::exp_approximate(v);
+    constexpr float accuracy = .001f;
+
+    EXPECT(fabsf(result[0] - 1.22140276f) <= accuracy);
+    EXPECT(fabsf(result[1] - 1.49182470f) <= accuracy);
+    EXPECT(fabsf(result[2] - 1.82211880f) <= accuracy);
+    EXPECT(fabsf(result[3] - 2.22554093f) <= accuracy);
+}