2 * This file is part of the GROMACS molecular simulation package.
4 * Copyright (c) 2014,2015,2016,2017,2018, by the GROMACS development team, led by
5 * Mark Abraham, David van der Spoel, Berk Hess, and Erik Lindahl,
6 * and including many others, as listed in the AUTHORS file in the
7 * top-level source directory and at http://www.gromacs.org.
9 * GROMACS is free software; you can redistribute it and/or
10 * modify it under the terms of the GNU Lesser General Public License
11 * as published by the Free Software Foundation; either version 2.1
12 * of the License, or (at your option) any later version.
14 * GROMACS is distributed in the hope that it will be useful,
15 * but WITHOUT ANY WARRANTY; without even the implied warranty of
16 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
17 * Lesser General Public License for more details.
19 * You should have received a copy of the GNU Lesser General Public
20 * License along with GROMACS; if not, see
21 * http://www.gnu.org/licenses, or write to the Free Software Foundation,
22 * Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.
24 * If you want to redistribute modifications to GROMACS, please
25 * consider that scientific software is very special. Version
26 * control is crucial - bugs must be traceable. We will be happy to
27 * consider code for inclusion in the official distribution, but
28 * derived work must not be called official GROMACS. Details are found
29 * in the README & COPYING files - if they are missing, get the
30 * official version at http://www.gromacs.org.
32 * To help us fund GROMACS development, we humbly ask that you cite
33 * the research papers on the package. Check out http://www.gromacs.org.
36 #ifndef GMX_SIMD_IMPL_X86_AVX_256_SIMD_FLOAT_H
37 #define GMX_SIMD_IMPL_X86_AVX_256_SIMD_FLOAT_H
45 #include <immintrin.h>
47 #include "gromacs/math/utilities.h"
57 SimdFloat(float f
) : simdInternal_(_mm256_set1_ps(f
)) {}
59 // Internal utility constructor to simplify return statements
60 SimdFloat(__m256 simd
) : simdInternal_(simd
) {}
70 SimdFInt32(std::int32_t i
) : simdInternal_(_mm256_set1_epi32(i
)) {}
72 // Internal utility constructor to simplify return statements
73 SimdFInt32(__m256i simd
) : simdInternal_(simd
) {}
75 __m256i simdInternal_
;
83 SimdFBool(bool b
) : simdInternal_(_mm256_castsi256_ps(_mm256_set1_epi32( b
? 0xFFFFFFFF : 0))) {}
85 // Internal utility constructor to simplify return statements
86 SimdFBool(__m256 simd
) : simdInternal_(simd
) {}
91 static inline SimdFloat gmx_simdcall
92 simdLoad(const float *m
, SimdFloatTag
/*unused*/ = {})
94 assert(std::size_t(m
) % 32 == 0);
100 static inline void gmx_simdcall
101 store(float *m
, SimdFloat a
)
103 assert(std::size_t(m
) % 32 == 0);
104 _mm256_store_ps(m
, a
.simdInternal_
);
107 static inline SimdFloat gmx_simdcall
108 simdLoadU(const float *m
, SimdFloatTag
/*unused*/ = {})
115 static inline void gmx_simdcall
116 storeU(float *m
, SimdFloat a
)
118 _mm256_storeu_ps(m
, a
.simdInternal_
);
121 static inline SimdFloat gmx_simdcall
129 static inline SimdFInt32 gmx_simdcall
130 simdLoad(const std::int32_t * m
, SimdFInt32Tag
/*unused*/)
132 assert(std::size_t(m
) % 32 == 0);
134 _mm256_load_si256(reinterpret_cast<const __m256i
*>(m
))
138 static inline void gmx_simdcall
139 store(std::int32_t * m
, SimdFInt32 a
)
141 assert(std::size_t(m
) % 32 == 0);
142 _mm256_store_si256(reinterpret_cast<__m256i
*>(m
), a
.simdInternal_
);
145 static inline SimdFInt32 gmx_simdcall
146 simdLoadU(const std::int32_t *m
, SimdFInt32Tag
/*unused*/)
149 _mm256_loadu_si256(reinterpret_cast<const __m256i
*>(m
))
153 static inline void gmx_simdcall
154 storeU(std::int32_t * m
, SimdFInt32 a
)
156 _mm256_storeu_si256(reinterpret_cast<__m256i
*>(m
), a
.simdInternal_
);
159 static inline SimdFInt32 gmx_simdcall
163 _mm256_setzero_si256()
168 static inline std::int32_t gmx_simdcall
169 extract(SimdFInt32 a
)
171 return _mm_extract_epi32(_mm256_extractf128_si256(a
.simdInternal_
, index
>>2), index
& 0x3);
174 static inline SimdFloat gmx_simdcall
175 operator&(SimdFloat a
, SimdFloat b
)
178 _mm256_and_ps(a
.simdInternal_
, b
.simdInternal_
)
182 static inline SimdFloat gmx_simdcall
183 andNot(SimdFloat a
, SimdFloat b
)
186 _mm256_andnot_ps(a
.simdInternal_
, b
.simdInternal_
)
190 static inline SimdFloat gmx_simdcall
191 operator|(SimdFloat a
, SimdFloat b
)
194 _mm256_or_ps(a
.simdInternal_
, b
.simdInternal_
)
198 static inline SimdFloat gmx_simdcall
199 operator^(SimdFloat a
, SimdFloat b
)
202 _mm256_xor_ps(a
.simdInternal_
, b
.simdInternal_
)
206 static inline SimdFloat gmx_simdcall
207 operator+(SimdFloat a
, SimdFloat b
)
210 _mm256_add_ps(a
.simdInternal_
, b
.simdInternal_
)
214 static inline SimdFloat gmx_simdcall
215 operator-(SimdFloat a
, SimdFloat b
)
218 _mm256_sub_ps(a
.simdInternal_
, b
.simdInternal_
)
222 static inline SimdFloat gmx_simdcall
223 operator-(SimdFloat x
)
226 _mm256_xor_ps(x
.simdInternal_
, _mm256_set1_ps(GMX_FLOAT_NEGZERO
))
230 static inline SimdFloat gmx_simdcall
231 operator*(SimdFloat a
, SimdFloat b
)
234 _mm256_mul_ps(a
.simdInternal_
, b
.simdInternal_
)
238 // Override for AVX2 and higher
239 #if GMX_SIMD_X86_AVX_256
240 static inline SimdFloat gmx_simdcall
241 fma(SimdFloat a
, SimdFloat b
, SimdFloat c
)
244 _mm256_add_ps(_mm256_mul_ps(a
.simdInternal_
, b
.simdInternal_
), c
.simdInternal_
)
248 static inline SimdFloat gmx_simdcall
249 fms(SimdFloat a
, SimdFloat b
, SimdFloat c
)
252 _mm256_sub_ps(_mm256_mul_ps(a
.simdInternal_
, b
.simdInternal_
), c
.simdInternal_
)
256 static inline SimdFloat gmx_simdcall
257 fnma(SimdFloat a
, SimdFloat b
, SimdFloat c
)
260 _mm256_sub_ps(c
.simdInternal_
, _mm256_mul_ps(a
.simdInternal_
, b
.simdInternal_
))
264 static inline SimdFloat gmx_simdcall
265 fnms(SimdFloat a
, SimdFloat b
, SimdFloat c
)
268 _mm256_sub_ps(_mm256_setzero_ps(), _mm256_add_ps(_mm256_mul_ps(a
.simdInternal_
, b
.simdInternal_
), c
.simdInternal_
))
273 static inline SimdFloat gmx_simdcall
277 _mm256_rsqrt_ps(x
.simdInternal_
)
281 static inline SimdFloat gmx_simdcall
285 _mm256_rcp_ps(x
.simdInternal_
)
289 static inline SimdFloat gmx_simdcall
290 maskAdd(SimdFloat a
, SimdFloat b
, SimdFBool m
)
293 _mm256_add_ps(a
.simdInternal_
, _mm256_and_ps(b
.simdInternal_
, m
.simdInternal_
))
297 static inline SimdFloat gmx_simdcall
298 maskzMul(SimdFloat a
, SimdFloat b
, SimdFBool m
)
301 _mm256_and_ps(_mm256_mul_ps(a
.simdInternal_
, b
.simdInternal_
), m
.simdInternal_
)
305 static inline SimdFloat
306 maskzFma(SimdFloat a
, SimdFloat b
, SimdFloat c
, SimdFBool m
)
309 _mm256_and_ps(_mm256_add_ps(_mm256_mul_ps(a
.simdInternal_
, b
.simdInternal_
), c
.simdInternal_
), m
.simdInternal_
)
313 static inline SimdFloat
314 maskzRsqrt(SimdFloat x
, SimdFBool m
)
317 x
.simdInternal_
= _mm256_blendv_ps(_mm256_set1_ps(1.0f
), x
.simdInternal_
, m
.simdInternal_
);
320 _mm256_and_ps(_mm256_rsqrt_ps(x
.simdInternal_
), m
.simdInternal_
)
324 static inline SimdFloat
325 maskzRcp(SimdFloat x
, SimdFBool m
)
328 x
.simdInternal_
= _mm256_blendv_ps(_mm256_set1_ps(1.0f
), x
.simdInternal_
, m
.simdInternal_
);
331 _mm256_and_ps(_mm256_rcp_ps(x
.simdInternal_
), m
.simdInternal_
)
335 static inline SimdFloat gmx_simdcall
339 _mm256_andnot_ps( _mm256_set1_ps(GMX_FLOAT_NEGZERO
), x
.simdInternal_
)
343 static inline SimdFloat gmx_simdcall
344 max(SimdFloat a
, SimdFloat b
)
347 _mm256_max_ps(a
.simdInternal_
, b
.simdInternal_
)
351 static inline SimdFloat gmx_simdcall
352 min(SimdFloat a
, SimdFloat b
)
355 _mm256_min_ps(a
.simdInternal_
, b
.simdInternal_
)
359 static inline SimdFloat gmx_simdcall
363 _mm256_round_ps(x
.simdInternal_
, _MM_FROUND_NINT
)
367 static inline SimdFloat gmx_simdcall
371 _mm256_round_ps(x
.simdInternal_
, _MM_FROUND_TRUNC
)
375 // Override for AVX2 and higher
376 #if GMX_SIMD_X86_AVX_256
377 static inline SimdFloat gmx_simdcall
378 frexp(SimdFloat value
, SimdFInt32
* exponent
)
380 const __m256 exponentMask
= _mm256_castsi256_ps(_mm256_set1_epi32(0x7F800000));
381 const __m256 mantissaMask
= _mm256_castsi256_ps(_mm256_set1_epi32(0x807FFFFF));
382 const __m256 half
= _mm256_set1_ps(0.5);
383 const __m128i exponentBias
= _mm_set1_epi32(126); // add 1 to make our definition identical to frexp()
385 __m128i iExponentLow
, iExponentHigh
;
387 iExponent
= _mm256_castps_si256(_mm256_and_ps(value
.simdInternal_
, exponentMask
));
388 iExponentHigh
= _mm256_extractf128_si256(iExponent
, 0x1);
389 iExponentLow
= _mm256_castsi256_si128(iExponent
);
390 iExponentLow
= _mm_srli_epi32(iExponentLow
, 23);
391 iExponentHigh
= _mm_srli_epi32(iExponentHigh
, 23);
392 iExponentLow
= _mm_sub_epi32(iExponentLow
, exponentBias
);
393 iExponentHigh
= _mm_sub_epi32(iExponentHigh
, exponentBias
);
394 iExponent
= _mm256_castsi128_si256(iExponentLow
);
395 exponent
->simdInternal_
= _mm256_insertf128_si256(iExponent
, iExponentHigh
, 0x1);
398 _mm256_or_ps(_mm256_and_ps(value
.simdInternal_
, mantissaMask
), half
)
403 template <MathOptimization opt
= MathOptimization::Safe
>
404 static inline SimdFloat gmx_simdcall
405 ldexp(SimdFloat value
, SimdFInt32 exponent
)
407 const __m128i exponentBias
= _mm_set1_epi32(127);
409 __m128i iExponentLow
, iExponentHigh
;
411 iExponentHigh
= _mm256_extractf128_si256(exponent
.simdInternal_
, 0x1);
412 iExponentLow
= _mm256_castsi256_si128(exponent
.simdInternal_
);
414 iExponentLow
= _mm_add_epi32(iExponentLow
, exponentBias
);
415 iExponentHigh
= _mm_add_epi32(iExponentHigh
, exponentBias
);
417 if (opt
== MathOptimization::Safe
)
419 // Make sure biased argument is not negative
420 iExponentLow
= _mm_max_epi32(iExponentLow
, _mm_setzero_si128());
421 iExponentHigh
= _mm_max_epi32(iExponentHigh
, _mm_setzero_si128());
424 iExponentLow
= _mm_slli_epi32(iExponentLow
, 23);
425 iExponentHigh
= _mm_slli_epi32(iExponentHigh
, 23);
426 iExponent
= _mm256_castsi128_si256(iExponentLow
);
427 iExponent
= _mm256_insertf128_si256(iExponent
, iExponentHigh
, 0x1);
429 _mm256_mul_ps(value
.simdInternal_
, _mm256_castsi256_ps(iExponent
))
434 static inline float gmx_simdcall
438 t0
= _mm_add_ps(_mm256_castps256_ps128(a
.simdInternal_
), _mm256_extractf128_ps(a
.simdInternal_
, 0x1));
439 t0
= _mm_add_ps(t0
, _mm_permute_ps(t0
, _MM_SHUFFLE(1, 0, 3, 2)));
440 t0
= _mm_add_ss(t0
, _mm_permute_ps(t0
, _MM_SHUFFLE(0, 3, 2, 1)));
441 return *reinterpret_cast<float *>(&t0
);
444 static inline SimdFBool gmx_simdcall
445 operator==(SimdFloat a
, SimdFloat b
)
448 _mm256_cmp_ps(a
.simdInternal_
, b
.simdInternal_
, _CMP_EQ_OQ
)
452 static inline SimdFBool gmx_simdcall
453 operator!=(SimdFloat a
, SimdFloat b
)
456 _mm256_cmp_ps(a
.simdInternal_
, b
.simdInternal_
, _CMP_NEQ_OQ
)
460 static inline SimdFBool gmx_simdcall
461 operator<(SimdFloat a
, SimdFloat b
)
464 _mm256_cmp_ps(a
.simdInternal_
, b
.simdInternal_
, _CMP_LT_OQ
)
468 static inline SimdFBool gmx_simdcall
469 operator<=(SimdFloat a
, SimdFloat b
)
472 _mm256_cmp_ps(a
.simdInternal_
, b
.simdInternal_
, _CMP_LE_OQ
)
476 // Override for AVX2 and higher
477 #if GMX_SIMD_X86_AVX_256
478 static inline SimdFBool gmx_simdcall
479 testBits(SimdFloat a
)
481 __m256 tst
= _mm256_cvtepi32_ps(_mm256_castps_si256(a
.simdInternal_
));
484 _mm256_cmp_ps(tst
, _mm256_setzero_ps(), _CMP_NEQ_OQ
)
489 static inline SimdFBool gmx_simdcall
490 operator&&(SimdFBool a
, SimdFBool b
)
493 _mm256_and_ps(a
.simdInternal_
, b
.simdInternal_
)
497 static inline SimdFBool gmx_simdcall
498 operator||(SimdFBool a
, SimdFBool b
)
501 _mm256_or_ps(a
.simdInternal_
, b
.simdInternal_
)
505 static inline bool gmx_simdcall
506 anyTrue(SimdFBool a
) { return _mm256_movemask_ps(a
.simdInternal_
) != 0; }
508 static inline SimdFloat gmx_simdcall
509 selectByMask(SimdFloat a
, SimdFBool mask
)
512 _mm256_and_ps(a
.simdInternal_
, mask
.simdInternal_
)
516 static inline SimdFloat gmx_simdcall
517 selectByNotMask(SimdFloat a
, SimdFBool mask
)
520 _mm256_andnot_ps(mask
.simdInternal_
, a
.simdInternal_
)
524 static inline SimdFloat gmx_simdcall
525 blend(SimdFloat a
, SimdFloat b
, SimdFBool sel
)
528 _mm256_blendv_ps(a
.simdInternal_
, b
.simdInternal_
, sel
.simdInternal_
)
532 static inline SimdFInt32 gmx_simdcall
536 _mm256_cvtps_epi32(a
.simdInternal_
)
540 static inline SimdFInt32 gmx_simdcall
544 _mm256_cvttps_epi32(a
.simdInternal_
)
548 static inline SimdFloat gmx_simdcall
552 _mm256_cvtepi32_ps(a
.simdInternal_
)
558 #endif // GMX_SIMD_IMPL_X86_AVX_256_SIMD_FLOAT_H