10#ifndef THIRD_PARTY_EIGEN3_EIGEN_SRC_CORE_ARCH_AVX512_MATHFUNCTIONS_H_
11#define THIRD_PARTY_EIGEN3_EIGEN_SRC_CORE_ARCH_AVX512_MATHFUNCTIONS_H_
18#if EIGEN_GNUC_AT_LEAST(5, 3) || EIGEN_COMP_CLANG || EIGEN_COMP_MSVC >= 1923
20#define _EIGEN_DECLARE_CONST_Packet16f(NAME, X) \
21 const Packet16f p16f_##NAME = pset1<Packet16f>(X)
23#define _EIGEN_DECLARE_CONST_Packet16f_FROM_INT(NAME, X) \
24 const Packet16f p16f_##NAME = preinterpret<Packet16f,Packet16i>(pset1<Packet16i>(X))
26#define _EIGEN_DECLARE_CONST_Packet8d(NAME, X) \
27 const Packet8d p8d_##NAME = pset1<Packet8d>(X)
29#define _EIGEN_DECLARE_CONST_Packet8d_FROM_INT64(NAME, X) \
30 const Packet8d p8d_##NAME = _mm512_castsi512_pd(_mm512_set1_epi64(X))
32#define _EIGEN_DECLARE_CONST_Packet16bf(NAME, X) \
33 const Packet16bf p16bf_##NAME = pset1<Packet16bf>(X)
35#define _EIGEN_DECLARE_CONST_Packet16bf_FROM_INT(NAME, X) \
36 const Packet16bf p16bf_##NAME = preinterpret<Packet16bf,Packet16i>(pset1<Packet16i>(X))
39EIGEN_DEFINE_FUNCTION_ALLOWING_MULTIPLE_DEFINITIONS EIGEN_UNUSED Packet16f
40plog<Packet16f>(
const Packet16f& _x) {
41 return plog_float(_x);
45EIGEN_DEFINE_FUNCTION_ALLOWING_MULTIPLE_DEFINITIONS EIGEN_UNUSED Packet8d
46plog<Packet8d>(
const Packet8d& _x) {
47 return plog_double(_x);
50F16_PACKET_FUNCTION(Packet16f, Packet16h, plog)
51BF16_PACKET_FUNCTION(Packet16f, Packet16bf, plog)
54EIGEN_DEFINE_FUNCTION_ALLOWING_MULTIPLE_DEFINITIONS EIGEN_UNUSED Packet16f
55plog2<Packet16f>(
const Packet16f& _x) {
56 return plog2_float(_x);
60EIGEN_DEFINE_FUNCTION_ALLOWING_MULTIPLE_DEFINITIONS EIGEN_UNUSED Packet8d
61plog2<Packet8d>(
const Packet8d& _x) {
62 return plog2_double(_x);
65F16_PACKET_FUNCTION(Packet16f, Packet16h, plog2)
66BF16_PACKET_FUNCTION(Packet16f, Packet16bf, plog2)
72EIGEN_DEFINE_FUNCTION_ALLOWING_MULTIPLE_DEFINITIONS EIGEN_UNUSED Packet16f
73pexp<Packet16f>(
const Packet16f& _x) {
74 _EIGEN_DECLARE_CONST_Packet16f(1, 1.0f);
75 _EIGEN_DECLARE_CONST_Packet16f(half, 0.5f);
76 _EIGEN_DECLARE_CONST_Packet16f(127, 127.0f);
78 _EIGEN_DECLARE_CONST_Packet16f(exp_hi, 88.3762626647950f);
79 _EIGEN_DECLARE_CONST_Packet16f(exp_lo, -88.3762626647949f);
81 _EIGEN_DECLARE_CONST_Packet16f(cephes_LOG2EF, 1.44269504088896341f);
83 _EIGEN_DECLARE_CONST_Packet16f(cephes_exp_p0, 1.9875691500E-4f);
84 _EIGEN_DECLARE_CONST_Packet16f(cephes_exp_p1, 1.3981999507E-3f);
85 _EIGEN_DECLARE_CONST_Packet16f(cephes_exp_p2, 8.3334519073E-3f);
86 _EIGEN_DECLARE_CONST_Packet16f(cephes_exp_p3, 4.1665795894E-2f);
87 _EIGEN_DECLARE_CONST_Packet16f(cephes_exp_p4, 1.6666665459E-1f);
88 _EIGEN_DECLARE_CONST_Packet16f(cephes_exp_p5, 5.0000001201E-1f);
91 Packet16f x = pmax(pmin(_x, p16f_exp_hi), p16f_exp_lo);
95 Packet16f m = _mm512_floor_ps(pmadd(x, p16f_cephes_LOG2EF, p16f_half));
99 _EIGEN_DECLARE_CONST_Packet16f(nln2, -0.6931471805599453f);
100 Packet16f r = _mm512_fmadd_ps(m, p16f_nln2, x);
101 Packet16f r2 = pmul(r, r);
102 Packet16f r3 = pmul(r2, r);
106 y = pmadd(p16f_cephes_exp_p0, r, p16f_cephes_exp_p1);
107 y1 = pmadd(p16f_cephes_exp_p3, r, p16f_cephes_exp_p4);
108 y2 = padd(r, p16f_1);
109 y = pmadd(y, r, p16f_cephes_exp_p2);
110 y1 = pmadd(y1, r, p16f_cephes_exp_p5);
111 y = pmadd(y, r3, y1);
112 y = pmadd(y, r2, y2);
115 Packet16i emm0 = _mm512_cvttps_epi32(padd(m, p16f_127));
116 emm0 = _mm512_slli_epi32(emm0, 23);
119 return pmax(pmul(y, _mm512_castsi512_ps(emm0)), _x);
123EIGEN_DEFINE_FUNCTION_ALLOWING_MULTIPLE_DEFINITIONS EIGEN_UNUSED Packet8d
124pexp<Packet8d>(
const Packet8d& _x) {
125 return pexp_double(_x);
128F16_PACKET_FUNCTION(Packet16f, Packet16h, pexp)
129BF16_PACKET_FUNCTION(Packet16f, Packet16bf, pexp)
132EIGEN_STRONG_INLINE Packet16h pfrexp(
const Packet16h& a, Packet16h& exponent) {
134 const Packet16h out = float2half(pfrexp<Packet16f>(half2float(a), fexponent));
135 exponent = float2half(fexponent);
140EIGEN_STRONG_INLINE Packet16h pldexp(
const Packet16h& a,
const Packet16h& exponent) {
141 return float2half(pldexp<Packet16f>(half2float(a), half2float(exponent)));
145EIGEN_STRONG_INLINE Packet16bf pfrexp(
const Packet16bf& a, Packet16bf& exponent) {
147 const Packet16bf out = F32ToBf16(pfrexp<Packet16f>(Bf16ToF32(a), fexponent));
148 exponent = F32ToBf16(fexponent);
153EIGEN_STRONG_INLINE Packet16bf pldexp(
const Packet16bf& a,
const Packet16bf& exponent) {
154 return F32ToBf16(pldexp<Packet16f>(Bf16ToF32(a), Bf16ToF32(exponent)));
165EIGEN_DEFINE_FUNCTION_ALLOWING_MULTIPLE_DEFINITIONS EIGEN_UNUSED Packet16f
166psqrt<Packet16f>(
const Packet16f& _x) {
167 Packet16f neg_half = pmul(_x, pset1<Packet16f>(-.5f));
168 __mmask16 denormal_mask = _mm512_kand(
169 _mm512_cmp_ps_mask(_x, pset1<Packet16f>((std::numeric_limits<float>::min)()),
171 _mm512_cmp_ps_mask(_x, _mm512_setzero_ps(), _CMP_GE_OQ));
173 Packet16f x = _mm512_rsqrt14_ps(_x);
176 x = pmul(x, pmadd(neg_half, pmul(x, x), pset1<Packet16f>(1.5f)));
179 return _mm512_mask_blend_ps(denormal_mask, pmul(_x,x), _mm512_setzero_ps());
183EIGEN_DEFINE_FUNCTION_ALLOWING_MULTIPLE_DEFINITIONS EIGEN_UNUSED Packet8d
184psqrt<Packet8d>(
const Packet8d& _x) {
185 Packet8d neg_half = pmul(_x, pset1<Packet8d>(-.5));
186 __mmask16 denormal_mask = _mm512_kand(
187 _mm512_cmp_pd_mask(_x, pset1<Packet8d>((std::numeric_limits<double>::min)()),
189 _mm512_cmp_pd_mask(_x, _mm512_setzero_pd(), _CMP_GE_OQ));
191 Packet8d x = _mm512_rsqrt14_pd(_x);
194 x = pmul(x, pmadd(neg_half, pmul(x, x), pset1<Packet8d>(1.5)));
197 x = pmul(x, pmadd(neg_half, pmul(x, x), pset1<Packet8d>(1.5)));
199 return _mm512_mask_blend_pd(denormal_mask, pmul(_x,x), _mm512_setzero_pd());
203EIGEN_STRONG_INLINE Packet16f psqrt<Packet16f>(
const Packet16f& x) {
204 return _mm512_sqrt_ps(x);
208EIGEN_STRONG_INLINE Packet8d psqrt<Packet8d>(
const Packet8d& x) {
209 return _mm512_sqrt_pd(x);
213F16_PACKET_FUNCTION(Packet16f, Packet16h, psqrt)
214BF16_PACKET_FUNCTION(Packet16f, Packet16bf, psqrt)
217#if defined(EIGEN_VECTORIZE_AVX512ER)
220EIGEN_STRONG_INLINE Packet16f prsqrt<Packet16f>(
const Packet16f& x) {
221 return _mm512_rsqrt28_ps(x);
226EIGEN_DEFINE_FUNCTION_ALLOWING_MULTIPLE_DEFINITIONS EIGEN_UNUSED Packet16f
227prsqrt<Packet16f>(
const Packet16f& _x) {
228 _EIGEN_DECLARE_CONST_Packet16f_FROM_INT(inf, 0x7f800000);
229 _EIGEN_DECLARE_CONST_Packet16f(one_point_five, 1.5f);
230 _EIGEN_DECLARE_CONST_Packet16f(minus_half, -0.5f);
232 Packet16f neg_half = pmul(_x, p16f_minus_half);
235 __mmask16 inf_mask = _mm512_cmp_ps_mask(_x, p16f_inf, _CMP_EQ_OQ);
236 __mmask16 not_pos_mask = _mm512_cmp_ps_mask(_x, _mm512_setzero_ps(), _CMP_LE_OQ);
237 __mmask16 not_finite_pos_mask = not_pos_mask | inf_mask;
241 Packet16f y_approx = _mm512_rsqrt14_ps(_x);
247 Packet16f y_newton = pmul(y_approx, pmadd(y_approx, pmul(neg_half, y_approx), p16f_one_point_five));
252 return _mm512_mask_blend_ps(not_finite_pos_mask, y_newton, y_approx);
257EIGEN_STRONG_INLINE Packet16f prsqrt<Packet16f>(
const Packet16f& x) {
258 _EIGEN_DECLARE_CONST_Packet16f(one, 1.0f);
259 return _mm512_div_ps(p16f_one, _mm512_sqrt_ps(x));
263F16_PACKET_FUNCTION(Packet16f, Packet16h, prsqrt)
264BF16_PACKET_FUNCTION(Packet16f, Packet16bf, prsqrt)
269EIGEN_DEFINE_FUNCTION_ALLOWING_MULTIPLE_DEFINITIONS EIGEN_UNUSED Packet8d
270prsqrt<Packet8d>(
const Packet8d& _x) {
271 _EIGEN_DECLARE_CONST_Packet8d(one_point_five, 1.5);
272 _EIGEN_DECLARE_CONST_Packet8d(minus_half, -0.5);
273 _EIGEN_DECLARE_CONST_Packet8d_FROM_INT64(inf, 0x7ff0000000000000LL);
275 Packet8d neg_half = pmul(_x, p8d_minus_half);
278 __mmask8 inf_mask = _mm512_cmp_pd_mask(_x, p8d_inf, _CMP_EQ_OQ);
279 __mmask8 not_pos_mask = _mm512_cmp_pd_mask(_x, _mm512_setzero_pd(), _CMP_LE_OQ);
280 __mmask8 not_finite_pos_mask = not_pos_mask | inf_mask;
284#if defined(EIGEN_VECTORIZE_AVX512ER)
285 Packet8d y_approx = _mm512_rsqrt28_pd(_x);
287 Packet8d y_approx = _mm512_rsqrt14_pd(_x);
296 Packet8d y_newton = pmul(y_approx, pmadd(neg_half, pmul(y_approx, y_approx), p8d_one_point_five));
297#if !defined(EIGEN_VECTORIZE_AVX512ER)
298 y_newton = pmul(y_newton, pmadd(y_newton, pmul(neg_half, y_newton), p8d_one_point_five));
303 return _mm512_mask_blend_pd(not_finite_pos_mask, y_newton, y_approx);
307EIGEN_STRONG_INLINE Packet8d prsqrt<Packet8d>(
const Packet8d& x) {
308 _EIGEN_DECLARE_CONST_Packet8d(one, 1.0f);
309 return _mm512_div_pd(p8d_one, _mm512_sqrt_pd(x));
313template<> EIGEN_DEFINE_FUNCTION_ALLOWING_MULTIPLE_DEFINITIONS EIGEN_UNUSED
314Packet16f plog1p<Packet16f>(
const Packet16f& _x) {
315 return generic_plog1p(_x);
318F16_PACKET_FUNCTION(Packet16f, Packet16h, plog1p)
319BF16_PACKET_FUNCTION(Packet16f, Packet16bf, plog1p)
321template<> EIGEN_DEFINE_FUNCTION_ALLOWING_MULTIPLE_DEFINITIONS EIGEN_UNUSED
322Packet16f pexpm1<Packet16f>(
const Packet16f& _x) {
323 return generic_expm1(_x);
326F16_PACKET_FUNCTION(Packet16f, Packet16h, pexpm1)
327BF16_PACKET_FUNCTION(Packet16f, Packet16bf, pexpm1)
333EIGEN_DEFINE_FUNCTION_ALLOWING_MULTIPLE_DEFINITIONS EIGEN_UNUSED Packet16f
334psin<Packet16f>(
const Packet16f& _x) {
335 return psin_float(_x);
339EIGEN_DEFINE_FUNCTION_ALLOWING_MULTIPLE_DEFINITIONS EIGEN_UNUSED Packet16f
340pcos<Packet16f>(
const Packet16f& _x) {
341 return pcos_float(_x);
345EIGEN_DEFINE_FUNCTION_ALLOWING_MULTIPLE_DEFINITIONS EIGEN_UNUSED Packet16f
346ptanh<Packet16f>(
const Packet16f& _x) {
347 return internal::generic_fast_tanh_float(_x);
350F16_PACKET_FUNCTION(Packet16f, Packet16h, psin)
351F16_PACKET_FUNCTION(Packet16f, Packet16h, pcos)
352F16_PACKET_FUNCTION(Packet16f, Packet16h, ptanh)
354BF16_PACKET_FUNCTION(Packet16f, Packet16bf, psin)
355BF16_PACKET_FUNCTION(Packet16f, Packet16bf, pcos)
356BF16_PACKET_FUNCTION(Packet16f, Packet16bf, ptanh)
Namespace containing all symbols from the Eigen library.
Definition: Core:141