diff options
author | Andrew Turner <andrew@FreeBSD.org> | 2025-01-10 10:34:52 +0000 |
---|---|---|
committer | Andrew Turner <andrew@FreeBSD.org> | 2025-01-10 10:39:34 +0000 |
commit | 9d1de25930735261c16ed874a933b4c1f1d9041e (patch) | |
tree | b0cac1c933cc1ecb885c7e757b89ffbf13f1f012 /pl/math/sv_powif.c | |
parent | edc5c0de794f521eb620d2b6cbaee2434442a8f3 (diff) |
Update the Arm Optimized Routinesvendor/arm-optimized-routines/v25.01vendor/arm-optimized-routines
Import the v25.01 release of the Arm Optimized Routines [1].
[1] https://github.com/ARM-software/optimized-routines/tree/v25.01
Sponsored by: Arm Ltd
Diffstat (limited to 'pl/math/sv_powif.c')
-rw-r--r-- | pl/math/sv_powif.c | 48 |
1 files changed, 0 insertions, 48 deletions
diff --git a/pl/math/sv_powif.c b/pl/math/sv_powif.c deleted file mode 100644 index 7e032fd86a20..000000000000 --- a/pl/math/sv_powif.c +++ /dev/null @@ -1,48 +0,0 @@ -/* - * Single-precision SVE powi(x, n) function. - * - * Copyright (c) 2020-2023, Arm Limited. - * SPDX-License-Identifier: MIT OR Apache-2.0 WITH LLVM-exception - */ - -#include "sv_math.h" - -/* Optimized single-precision vector powi (float base, integer power). - powi is developed for environments in which accuracy is of much less - importance than performance, hence we provide no estimate for worst-case - error. */ -svfloat32_t -_ZGVsMxvv_powi (svfloat32_t as, svint32_t ns, svbool_t p) -{ - /* Compute powi by successive squaring, right to left. */ - svfloat32_t acc = sv_f32 (1.f); - svbool_t want_recip = svcmplt (p, ns, 0); - svuint32_t ns_abs = svreinterpret_u32 (svabs_x (p, ns)); - - /* We use a max to avoid needing to check whether any lane != 0 on each - iteration. */ - uint32_t max_n = svmaxv (p, ns_abs); - - svfloat32_t c = as; - /* Successively square c, and use merging predication (_m) to determine - whether or not to perform the multiplication or keep the previous - iteration. */ - while (true) - { - svbool_t px = svcmpeq (p, svand_x (p, ns_abs, 1), 1); - acc = svmul_m (px, acc, c); - max_n >>= 1; - if (max_n == 0) - break; - - ns_abs = svlsr_x (p, ns_abs, 1); - c = svmul_x (p, c, c); - } - - /* Negative powers are handled by computing the abs(n) version and then - taking the reciprocal. */ - if (svptest_any (want_recip, want_recip)) - acc = svdivr_m (want_recip, acc, 1.0f); - - return acc; -} |