-
Notifications
You must be signed in to change notification settings - Fork 9
Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
Setup initial work for pow implementation
- Loading branch information
Showing
19 changed files
with
555 additions
and
432 deletions.
There are no files selected for viewing
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
32 changes: 32 additions & 0 deletions
32
include/ccmath/internal/math/generic/func/power/pow_impl.hpp
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,32 @@ | ||
| /* | ||
| * Copyright (c) Ian Pike | ||
| * Copyright (c) CCMath contributors | ||
| * | ||
| * CCMath is provided under the Apache-2.0 License WITH LLVM-exception. | ||
| * See LICENSE for more information. | ||
| * | ||
| * SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception | ||
| */ | ||
|
|
||
| #pragma once | ||
|
|
||
| #include <type_traits> | ||
|
|
||
| namespace ccm::gen::impl | ||
| { | ||
| namespace internal::impl | ||
| { | ||
| template <typename T> | ||
| constexpr T pow_impl(T base, T exp) noexcept | ||
| { | ||
| return 0; | ||
| } | ||
| } | ||
|
|
||
| template <typename T> | ||
| constexpr T pow_impl(T base, T exp) noexcept | ||
| { | ||
| return internal::impl::pow_impl(base, exp); | ||
| } | ||
|
|
||
| } // namespace ccm |
32 changes: 32 additions & 0 deletions
32
include/ccmath/internal/math/generic/func/power/powf_impl.hpp
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,32 @@ | ||
| /* | ||
| * Copyright (c) Ian Pike | ||
| * Copyright (c) CCMath contributors | ||
| * | ||
| * CCMath is provided under the Apache-2.0 License WITH LLVM-exception. | ||
| * See LICENSE for more information. | ||
| * | ||
| * SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception | ||
| */ | ||
|
|
||
| #pragma once | ||
|
|
||
| #include <type_traits> | ||
|
|
||
| namespace ccm::gen::impl | ||
| { | ||
| namespace internal::impl | ||
| { | ||
| template <typename T> | ||
| constexpr T powf_impl(T base, T exp) noexcept | ||
| { | ||
| return 0; | ||
| } | ||
| } | ||
|
|
||
| template <typename T> | ||
| constexpr T powf_impl(T base, T exp) noexcept | ||
| { | ||
| return internal::impl::powf_impl(base, exp); | ||
| } | ||
|
|
||
| } // namespace ccm |
73 changes: 73 additions & 0 deletions
73
include/ccmath/internal/math/runtime/func/power/pow_rt.hpp
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,73 @@ | ||
| /* | ||
| * Copyright (c) Ian Pike | ||
| * Copyright (c) CCMath contributors | ||
| * | ||
| * CCMath is provided under the Apache-2.0 License WITH LLVM-exception. | ||
| * See LICENSE for more information. | ||
| * | ||
| * SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception | ||
| */ | ||
|
|
||
| #pragma once | ||
|
|
||
| #include "ccmath/internal/config/type_support.hpp" | ||
| #include "ccmath/internal/math/generic/func/power/pow_gen.hpp" | ||
| #include "ccmath/internal/math/runtime/simd/func/sqrt.hpp" | ||
| #include "ccmath/internal/support/always_false.hpp" | ||
| #include "ccmath/internal/support/bits.hpp" | ||
| #include "ccmath/internal/support/fenv/rounding_mode.hpp" | ||
| #include "ccmath/internal/support/fp/fp_bits.hpp" | ||
| #include "ccmath/internal/support/is_constant_evaluated.hpp" | ||
|
|
||
| #include <type_traits> | ||
|
|
||
| namespace ccm::rt::simd_impl | ||
| { | ||
| #if defined(CCM_TYPES_LONG_DOUBLE_IS_FLOAT64) | ||
| template <typename T, std::enable_if_t<std::is_floating_point_v<T>, bool> = true> | ||
| #else | ||
| template <typename T, std::enable_if_t<std::is_floating_point_v<T> && !std::is_same_v<T, long double>, bool> = true> | ||
| #endif | ||
| [[nodiscard]] inline T pow_simd_impl(T num) noexcept; | ||
|
|
||
| #ifdef CCMATH_HAS_SIMD | ||
| #if defined(CCM_TYPES_LONG_DOUBLE_IS_FLOAT64) | ||
| template <typename T, std::enable_if_t<std::is_floating_point_v<T>, bool>> | ||
| #else | ||
| template <typename T, std::enable_if_t<std::is_floating_point_v<T> && !std::is_same_v<T, long double>, bool>> | ||
| #endif | ||
| [[nodiscard]] inline T pow_simd_impl(T num) noexcept | ||
| { | ||
| intrin::simd<T, intrin::abi::native> const num_m(num); | ||
| intrin::simd<T, intrin::abi::native> const sqrt_m = intrin::sqrt(num_m); | ||
| return sqrt_m.convert(); | ||
| } | ||
| #endif | ||
| } // namespace ccm::rt::simd_impl | ||
|
|
||
| namespace ccm::rt | ||
| { | ||
| template <typename T, std::enable_if_t<std::is_floating_point_v<T>, bool> = true> | ||
| T pow_rt(T num) | ||
| { | ||
| #if CCM_HAS_BUILTIN(__builtin_pow) || defined(__builtin_pow) // Prefer the builtins if available. | ||
| if constexpr (std::is_same_v<T, float>) { return __builtin_powf(num); } | ||
| else if constexpr (std::is_same_v<T, double>) { return __builtin_pow(num); } | ||
| else if constexpr (std::is_same_v<T, long double>) { return __builtin_powl(num); } | ||
| else { return static_cast<T>(__builtin_powl(static_cast<long double>(num))); } | ||
| #elif defined(CCMATH_HAS_SIMD) | ||
| // In the unlikely event, the rounding mode is not the default, use the runtime implementation instead. | ||
| if (CCM_UNLIKELY(ccm::support::fenv::get_rounding_mode() != FE_TONEAREST)) { return gen::pow_gen<T>(num); } | ||
| #if !defined(CCM_TYPES_LONG_DOUBLE_IS_FLOAT64) // If long double is different from double, use the generic implementation instead. | ||
| if constexpr (std::is_same_v<T, float> || std::is_same_v<T, double>) { return simd_impl::pow_simd_impl(num); } | ||
| else { return gen::pow_gen<T>(num); } | ||
| #else // If long double is the same as double, we can use the SIMD implementation instead. | ||
| if constexpr (std::is_same_v<T, float> || std::is_same_v<T, double>) { return simd_impl::sqrt_simd_impl(num); } | ||
| else if constexpr (std::is_same_v<T, long double>) { return static_cast<long double>(simd_impl::sqrt_simd_impl(static_cast<double>(num))); } | ||
| else { return ccm::gen::sqrt_gen<T>(num); } | ||
| #endif | ||
| #else // If we don't have a builtin or SIMD, use the generic implementation. | ||
| return ccm::gen::pow_gen<T>(num); | ||
| #endif | ||
| } | ||
| } // namespace ccm::rt |
Oops, something went wrong.