blob: c54d8cd249947d568be507d2920a20ea8ceaa8e4 [file] [log] [blame]
/*
* Copyright (c) 2015 Cisco and/or its affiliates.
* Licensed under the Apache License, Version 2.0 (the "License");
* you may not use this file except in compliance with the License.
* You may obtain a copy of the License at:
*
* http://www.apache.org/licenses/LICENSE-2.0
*
* Unless required by applicable law or agreed to in writing, software
* distributed under the License is distributed on an "AS IS" BASIS,
* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
* See the License for the specific language governing permissions and
* limitations under the License.
*/
#ifndef included_vector_avx512_h
#define included_vector_avx512_h
#include <vppinfra/clib.h>
#include <x86intrin.h>
/* *INDENT-OFF* */
#define foreach_avx512_vec512i \
_(i,8,64,epi8) _(i,16,32,epi16) _(i,32,16,epi32) _(i,64,8,epi64)
#define foreach_avx512_vec512u \
_(u,8,64,epi8) _(u,16,32,epi16) _(u,32,16,epi32) _(u,64,8,epi64)
#define foreach_avx512_vec512f \
_(f,32,8,ps) _(f,64,4,pd)
/* splat, load_unaligned, store_unaligned, is_all_zero, is_equal,
is_all_equal, is_zero_mask */
#define _(t, s, c, i) \
static_always_inline t##s##x##c \
t##s##x##c##_splat (t##s x) \
{ return (t##s##x##c) _mm512_set1_##i (x); } \
\
static_always_inline t##s##x##c \
t##s##x##c##_load_unaligned (void *p) \
{ return (t##s##x##c) _mm512_loadu_si512 (p); } \
\
static_always_inline void \
t##s##x##c##_store_unaligned (t##s##x##c v, void *p) \
{ _mm512_storeu_si512 ((__m512i *) p, (__m512i) v); } \
\
static_always_inline int \
t##s##x##c##_is_all_zero (t##s##x##c v) \
{ return (_mm512_test_epi64_mask ((__m512i) v, (__m512i) v) == 0); } \
\
static_always_inline int \
t##s##x##c##_is_equal (t##s##x##c a, t##s##x##c b) \
{ return t##s##x##c##_is_all_zero (a ^ b); } \
\
static_always_inline int \
t##s##x##c##_is_all_equal (t##s##x##c v, t##s x) \
{ return t##s##x##c##_is_equal (v, t##s##x##c##_splat (x)); } \
\
static_always_inline u##c \
t##s##x##c##_is_zero_mask (t##s##x##c v) \
{ return _mm512_test_##i##_mask ((__m512i) v, (__m512i) v); } \
\
static_always_inline t##s##x##c \
t##s##x##c##_interleave_lo (t##s##x##c a, t##s##x##c b) \
{ return (t##s##x##c) _mm512_unpacklo_##i ((__m512i) a, (__m512i) b); } \
\
static_always_inline t##s##x##c \
t##s##x##c##_interleave_hi (t##s##x##c a, t##s##x##c b) \
{ return (t##s##x##c) _mm512_unpackhi_##i ((__m512i) a, (__m512i) b); } \
foreach_avx512_vec512i foreach_avx512_vec512u
#undef _
/* *INDENT-ON* */
static_always_inline u32
u16x32_msb_mask (u16x32 v)
{
return (u32) _mm512_movepi16_mask ((__m512i) v);
}
static_always_inline u32x16
u32x16_byte_swap (u32x16 v)
{
u8x64 swap = {
3, 2, 1, 0, 7, 6, 5, 4, 11, 10, 9, 8, 15, 14, 13, 12,
3, 2, 1, 0, 7, 6, 5, 4, 11, 10, 9, 8, 15, 14, 13, 12,
3, 2, 1, 0, 7, 6, 5, 4, 11, 10, 9, 8, 15, 14, 13, 12,
3, 2, 1, 0, 7, 6, 5, 4, 11, 10, 9, 8, 15, 14, 13, 12
};
return (u32x16) _mm512_shuffle_epi8 ((__m512i) v, (__m512i) swap);
}
static_always_inline u16x32
u16x32_byte_swap (u16x32 v)
{
u8x64 swap = {
1, 0, 3, 2, 5, 4, 7, 6, 9, 8, 11, 10, 13, 12, 15, 14,
1, 0, 3, 2, 5, 4, 7, 6, 9, 8, 11, 10, 13, 12, 15, 14,
1, 0, 3, 2, 5, 4, 7, 6, 9, 8, 11, 10, 13, 12, 15, 14,
1, 0, 3, 2, 5, 4, 7, 6, 9, 8, 11, 10, 13, 12, 15, 14
};
return (u16x32) _mm512_shuffle_epi8 ((__m512i) v, (__m512i) swap);
}
static_always_inline u32x8
u32x16_extract_lo (u32x16 v)
{
return (u32x8) _mm512_extracti64x4_epi64 ((__m512i) v, 0);
}
static_always_inline u32x8
u32x16_extract_hi (u32x16 v)
{
return (u32x8) _mm512_extracti64x4_epi64 ((__m512i) v, 1);
}
static_always_inline u32
u32x16_min_scalar (u32x16 v)
{
return u32x8_min_scalar (u32x8_min (u32x16_extract_lo (v),
u32x16_extract_hi (v)));
}
static_always_inline u32x16
u32x16_insert_lo (u32x16 r, u32x8 v)
{
return (u32x16) _mm512_inserti64x4 ((__m512i) r, (__m256i) v, 0);
}
static_always_inline u32x16
u32x16_insert_hi (u32x16 r, u32x8 v)
{
return (u32x16) _mm512_inserti64x4 ((__m512i) r, (__m256i) v, 1);
}
static_always_inline u64x8
u64x8_permute (u64x8 a, u64x8 b, u64x8 mask)
{
return (u64x8) _mm512_permutex2var_epi64 ((__m512i) a, (__m512i) mask,
(__m512i) b);
}
#define u32x16_ternary_logic(a, b, c, d) \
(u32x16) _mm512_ternarylogic_epi32 ((__m512i) a, (__m512i) b, (__m512i) c, d)
static_always_inline void
u32x16_transpose (u32x16 m[16])
{
__m512i r[16], a, b, c, d, x, y;
/* *INDENT-OFF* */
__m512i pm1 = (__m512i) (u64x8) { 0, 1, 8, 9, 4, 5, 12, 13};
__m512i pm2 = (__m512i) (u64x8) { 2, 3, 10, 11, 6, 7, 14, 15};
__m512i pm3 = (__m512i) (u64x8) { 0, 1, 2, 3, 8, 9, 10, 11};
__m512i pm4 = (__m512i) (u64x8) { 4, 5, 6, 7, 12, 13, 14, 15};
/* *INDENT-ON* */
r[0] = _mm512_unpacklo_epi32 ((__m512i) m[0], (__m512i) m[1]);
r[1] = _mm512_unpacklo_epi32 ((__m512i) m[2], (__m512i) m[3]);
r[2] = _mm512_unpacklo_epi32 ((__m512i) m[4], (__m512i) m[5]);
r[3] = _mm512_unpacklo_epi32 ((__m512i) m[6], (__m512i) m[7]);
r[4] = _mm512_unpacklo_epi32 ((__m512i) m[8], (__m512i) m[9]);
r[5] = _mm512_unpacklo_epi32 ((__m512i) m[10], (__m512i) m[11]);
r[6] = _mm512_unpacklo_epi32 ((__m512i) m[12], (__m512i) m[13]);
r[7] = _mm512_unpacklo_epi32 ((__m512i) m[14], (__m512i) m[15]);
r[8] = _mm512_unpackhi_epi32 ((__m512i) m[0], (__m512i) m[1]);
r[9] = _mm512_unpackhi_epi32 ((__m512i) m[2], (__m512i) m[3]);
r[10] = _mm512_unpackhi_epi32 ((__m512i) m[4], (__m512i) m[5]);
r[11] = _mm512_unpackhi_epi32 ((__m512i) m[6], (__m512i) m[7]);
r[12] = _mm512_unpackhi_epi32 ((__m512i) m[8], (__m512i) m[9]);
r[13] = _mm512_unpackhi_epi32 ((__m512i) m[10], (__m512i) m[11]);
r[14] = _mm512_unpackhi_epi32 ((__m512i) m[12], (__m512i) m[13]);
r[15] = _mm512_unpackhi_epi32 ((__m512i) m[14], (__m512i) m[15]);
a = _mm512_unpacklo_epi64 (r[0], r[1]);
b = _mm512_unpacklo_epi64 (r[2], r[3]);
c = _mm512_unpacklo_epi64 (r[4], r[5]);
d = _mm512_unpacklo_epi64 (r[6], r[7]);
x = _mm512_permutex2var_epi64 (a, pm1, b);
y = _mm512_permutex2var_epi64 (c, pm1, d);
m[0] = (u32x16) _mm512_permutex2var_epi64 (x, pm3, y);
m[8] = (u32x16) _mm512_permutex2var_epi64 (x, pm4, y);
x = _mm512_permutex2var_epi64 (a, pm2, b);
y = _mm512_permutex2var_epi64 (c, pm2, d);
m[4] = (u32x16) _mm512_permutex2var_epi64 (x, pm3, y);
m[12] = (u32x16) _mm512_permutex2var_epi64 (x, pm4, y);
a = _mm512_unpacklo_epi64 (r[8], r[9]);
b = _mm512_unpacklo_epi64 (r[10], r[11]);
c = _mm512_unpacklo_epi64 (r[12], r[13]);
d = _mm512_unpacklo_epi64 (r[14], r[15]);
x = _mm512_permutex2var_epi64 (a, pm1, b);
y = _mm512_permutex2var_epi64 (c, pm1, d);
m[2] = (u32x16) _mm512_permutex2var_epi64 (x, pm3, y);
m[10] = (u32x16) _mm512_permutex2var_epi64 (x, pm4, y);
x = _mm512_permutex2var_epi64 (a, pm2, b);
y = _mm512_permutex2var_epi64 (c, pm2, d);
m[6] = (u32x16) _mm512_permutex2var_epi64 (x, pm3, y);
m[14] = (u32x16) _mm512_permutex2var_epi64 (x, pm4, y);
a = _mm512_unpackhi_epi64 (r[0], r[1]);
b = _mm512_unpackhi_epi64 (r[2], r[3]);
c = _mm512_unpackhi_epi64 (r[4], r[5]);
d = _mm512_unpackhi_epi64 (r[6], r[7]);
x = _mm512_permutex2var_epi64 (a, pm1, b);
y = _mm512_permutex2var_epi64 (c, pm1, d);
m[1] = (u32x16) _mm512_permutex2var_epi64 (x, pm3, y);
m[9] = (u32x16) _mm512_permutex2var_epi64 (x, pm4, y);
x = _mm512_permutex2var_epi64 (a, pm2, b);
y = _mm512_permutex2var_epi64 (c, pm2, d);
m[5] = (u32x16) _mm512_permutex2var_epi64 (x, pm3, y);
m[13] = (u32x16) _mm512_permutex2var_epi64 (x, pm4, y);
a = _mm512_unpackhi_epi64 (r[8], r[9]);
b = _mm512_unpackhi_epi64 (r[10], r[11]);
c = _mm512_unpackhi_epi64 (r[12], r[13]);
d = _mm512_unpackhi_epi64 (r[14], r[15]);
x = _mm512_permutex2var_epi64 (a, pm1, b);
y = _mm512_permutex2var_epi64 (c, pm1, d);
m[3] = (u32x16) _mm512_permutex2var_epi64 (x, pm3, y);
m[11] = (u32x16) _mm512_permutex2var_epi64 (x, pm4, y);
x = _mm512_permutex2var_epi64 (a, pm2, b);
y = _mm512_permutex2var_epi64 (c, pm2, d);
m[7] = (u32x16) _mm512_permutex2var_epi64 (x, pm3, y);
m[15] = (u32x16) _mm512_permutex2var_epi64 (x, pm4, y);
}
static_always_inline void
u64x8_transpose (u64x8 m[8])
{
__m512i r[8], x, y;
/* *INDENT-OFF* */
__m512i pm1 = (__m512i) (u64x8) { 0, 1, 8, 9, 4, 5, 12, 13};
__m512i pm2 = (__m512i) (u64x8) { 2, 3, 10, 11, 6, 7, 14, 15};
__m512i pm3 = (__m512i) (u64x8) { 0, 1, 2, 3, 8, 9, 10, 11};
__m512i pm4 = (__m512i) (u64x8) { 4, 5, 6, 7, 12, 13, 14, 15};
/* *INDENT-ON* */
r[0] = _mm512_unpacklo_epi64 ((__m512i) m[0], (__m512i) m[1]);
r[1] = _mm512_unpacklo_epi64 ((__m512i) m[2], (__m512i) m[3]);
r[2] = _mm512_unpacklo_epi64 ((__m512i) m[4], (__m512i) m[5]);
r[3] = _mm512_unpacklo_epi64 ((__m512i) m[6], (__m512i) m[7]);
r[4] = _mm512_unpackhi_epi64 ((__m512i) m[0], (__m512i) m[1]);
r[5] = _mm512_unpackhi_epi64 ((__m512i) m[2], (__m512i) m[3]);
r[6] = _mm512_unpackhi_epi64 ((__m512i) m[4], (__m512i) m[5]);
r[7] = _mm512_unpackhi_epi64 ((__m512i) m[6], (__m512i) m[7]);
x = _mm512_permutex2var_epi64 (r[0], pm1, r[1]);
y = _mm512_permutex2var_epi64 (r[2], pm1, r[3]);
m[0] = (u64x8) _mm512_permutex2var_epi64 (x, pm3, y);
m[4] = (u64x8) _mm512_permutex2var_epi64 (x, pm4, y);
x = _mm512_permutex2var_epi64 (r[0], pm2, r[1]);
y = _mm512_permutex2var_epi64 (r[2], pm2, r[3]);
m[2] = (u64x8) _mm512_permutex2var_epi64 (x, pm3, y);
m[6] = (u64x8) _mm512_permutex2var_epi64 (x, pm4, y);
x = _mm512_permutex2var_epi64 (r[4], pm1, r[5]);
y = _mm512_permutex2var_epi64 (r[6], pm1, r[7]);
m[1] = (u64x8) _mm512_permutex2var_epi64 (x, pm3, y);
m[5] = (u64x8) _mm512_permutex2var_epi64 (x, pm4, y);
x = _mm512_permutex2var_epi64 (r[4], pm2, r[5]);
y = _mm512_permutex2var_epi64 (r[6], pm2, r[7]);
m[3] = (u64x8) _mm512_permutex2var_epi64 (x, pm3, y);
m[7] = (u64x8) _mm512_permutex2var_epi64 (x, pm4, y);
}
#endif /* included_vector_avx512_h */
/*
* fd.io coding-style-patch-verification: ON
*
* Local Variables:
* eval: (c-set-style "gnu")
* End:
*/