| //===-- Common utils for exp function ---------------------------*- C++ -*-===// | |
| // | |
| // Part of the LLVM Project, under the Apache License v2.0 with LLVM Exceptions. | |
| // See https://llvm.org/LICENSE.txt for license information. | |
| // SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception | |
| // | |
| //===----------------------------------------------------------------------===// | |
| #ifndef LLVM_LIBC_SRC___SUPPORT_MATHVEC_EXP_UTILS_H | |
| #define LLVM_LIBC_SRC___SUPPORT_MATHVEC_EXP_UTILS_H | |
| #include "src/__support/CPP/simd.h" | |
| #include "src/__support/FPUtil/FPBits.h" | |
| #include "src/__support/mathvec/common_constants.h" | |
| namespace LIBC_NAMESPACE_DECL { | |
| namespace mathvec { | |
| template <size_t N> | |
| LIBC_INLINE static cpp::simd<double, N> exp_lookup(cpp::simd<uint64_t, N> u) { | |
| cpp::simd<uint64_t, N> index = u & cpp::simd<uint64_t, N>(0x3f); | |
| cpp::simd<uint64_t, N> mantissa = | |
| cpp::gather<cpp::simd<uint64_t, N>>(true, index, EXP_MANTISSA); | |
| cpp::simd<uint64_t, N> exponent = (u >> 6) << 52; | |
| cpp::simd<uint64_t, N> result = mantissa | exponent; | |
| return cpp::bit_cast<cpp::simd<double, N>>(result); | |
| } | |
| } // namespace mathvec | |
| } // namespace LIBC_NAMESPACE_DECL | |
| #endif // LLVM_LIBC_SRC___SUPPORT_MATHVEC_EXP_UTILS_H |