AK: Remove SIMDMath.h

This header is unused in ladybird. The associated test code is also
removed.
This commit is contained in:
R-Goc
2025-12-17 17:16:41 +01:00
committed by Gregory Bertilson
parent 0c37b093e9
commit 1665f4e980
3 changed files with 0 additions and 135 deletions

View File

@@ -1,98 +0,0 @@
/*
* Copyright (c) 2021, Stephan Unverwerth <s.unverwerth@serenityos.org>
*
* SPDX-License-Identifier: BSD-2-Clause
*/
#pragma once
#include <AK/Math.h>
#include <AK/SIMD.h>
#include <AK/SIMDExtras.h>
#include <math.h>
// Functions returning vectors or accepting vector arguments have different calling conventions
// depending on whether the target architecture supports SSE or not. GCC generates warning "psabi"
// when compiling for non-SSE architectures. We disable this warning because these functions
// are static and should never be visible from outside the translation unit that includes this header.
#pragma GCC diagnostic push
#pragma GCC diagnostic ignored "-Wpsabi"
namespace AK::SIMD {
// Functions ending in "_int_range" only accept arguments within range [INT_MIN, INT_MAX].
// Other inputs will generate unexpected results.
ALWAYS_INLINE static f32x4 truncate_int_range(f32x4 v)
{
return to_f32x4(to_i32x4(v));
}
ALWAYS_INLINE static f32x4 floor_int_range(f32x4 v)
{
auto t = truncate_int_range(v);
return t > v ? t - 1.0f : t;
}
ALWAYS_INLINE static f32x4 ceil_int_range(f32x4 v)
{
auto t = truncate_int_range(v);
return t < v ? t + 1.0f : t;
}
ALWAYS_INLINE static f32x4 frac_int_range(f32x4 v)
{
return v - floor_int_range(v);
}
ALWAYS_INLINE static f32x4 clamp(f32x4 v, f32x4 min, f32x4 max)
{
return v < min ? min : (v > max ? max : v);
}
ALWAYS_INLINE static f32x4 clamp(f32x4 v, float min, float max)
{
return v < min ? min : (v > max ? max : v);
}
ALWAYS_INLINE static f32x4 exp(f32x4 v)
{
// FIXME: This should be replaced with a vectorized algorithm instead of calling the scalar expf 4 times
return f32x4 {
expf(v[0]),
expf(v[1]),
expf(v[2]),
expf(v[3]),
};
}
ALWAYS_INLINE static f32x4 exp_approximate(f32x4 v)
{
static constexpr int number_of_iterations = 10;
auto result = 1.f + v / (1 << number_of_iterations);
for (int i = 0; i < number_of_iterations; ++i)
result *= result;
return result;
}
ALWAYS_INLINE static f32x4 sqrt(f32x4 v)
{
#if ARCH(X86_64)
# if defined(AK_COMPILER_CLANG)
return __builtin_elementwise_sqrt(v);
# else
return __builtin_ia32_sqrtps(v);
# endif
#else
return f32x4 {
AK::sqrt(v[0]),
AK::sqrt(v[1]),
AK::sqrt(v[2]),
AK::sqrt(v[3]),
};
#endif
}
}
#pragma GCC diagnostic pop

View File

@@ -58,7 +58,6 @@ set(AK_TEST_SOURCES
TestRedBlackTree.cpp
TestRefPtr.cpp
TestSegmentedVector.cpp
TestSIMD.cpp
TestSinglyLinkedList.cpp
TestSourceGenerator.cpp
TestSourceLocation.cpp

View File

@@ -1,36 +0,0 @@
/*
* Copyright (c) 2023, Jelle Raaijmakers <jelle@ladybird.org>
*
* SPDX-License-Identifier: BSD-2-Clause
*/
#include <LibTest/TestCase.h>
#include <AK/SIMD.h>
#include <AK/SIMDMath.h>
// See the comment in <AK/SIMDMath.h>
#pragma GCC diagnostic ignored "-Wpsabi"
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);
}