1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
|
/* { dg-do run } */
/* { dg-require-effective-target avx512er } */
/* { dg-options "-O2 -mavx512er" } */
#include "avx512er-check.h"
#include "avx512f-mask-type.h"
#include "avx512f-helper.h"
#include <math.h>
void static
compute_vrsqrt28ps (float *s, float *r)
{
int i;
for (i = 0; i < 16; i++)
r[i] = 1.0 / sqrt (s[i]);
}
void static
avx512er_test (void)
{
union512 src, res1, res2, res3;
__mmask16 mask = MASK_VALUE;
float res_ref[16];
int i;
for (i = 0; i < 16; i++)
{
src.a[i] = 179.345 - 6.5645 * i;
res2.a[i] = DEFAULT_VALUE;
}
res1.x = _mm512_rsqrt28_ps (src.x);
res2.x = _mm512_mask_rsqrt28_ps (res2.x, mask, src.x);
res3.x = _mm512_maskz_rsqrt28_ps (mask, src.x);
compute_vrsqrt28ps (src.a, res_ref);
if (check_rough_union512 (res1, res_ref, 0.0001))
abort ();
MASK_MERGE ()(res_ref, mask, 16);
if (check_rough_union512 (res2, res_ref, 0.0001))
abort ();
MASK_ZERO ()(res_ref, mask, 16);
if (check_rough_union512 (res3, res_ref, 0.0001))
abort ();
}
|
