1 /* { dg-do run } */
2 /* { dg-require-effective-target avx512er } */
3 /* { dg-options "-O2 -ffast-math -ftree-vectorize -mavx512er" } */
4
5 #include "avx512er-check.h"
6
7 #define MAX 1000
8 #define EPS 0.00001
9
10 __attribute__ ((noinline, optimize (0)))
11 void static
compute_rcp_ref(float * a,float * b,float * r)12 compute_rcp_ref (float *a, float *b, float *r)
13 {
14 for (int i = 0; i < MAX; i++)
15 r[i] = a[i] / b[i];
16 }
17
18 __attribute__ ((noinline))
19 void static
compute_rcp_exp(float * a,float * b,float * r)20 compute_rcp_exp (float *a, float *b, float *r)
21 {
22 for (int i = 0; i < MAX; i++)
23 r[i] = a[i] / b[i];
24 }
25
26 void static
avx512er_test(void)27 avx512er_test (void)
28 {
29 float a[MAX];
30 float b[MAX];
31 float ref[MAX];
32 float exp[MAX];
33
34 for (int i = 0; i < MAX; i++)
35 {
36 a[i] = 179.345 - 6.5645 * i;
37 b[i] = 8765.987 - 8.6756 * i;
38 }
39
40 compute_rcp_ref (a, b, ref);
41 compute_rcp_exp (a, b, exp);
42
43 for (int i = 0; i < MAX; i++)
44 {
45 float rel_err = (ref[i] - exp[i]) / ref[i];
46 rel_err = rel_err > 0.0 ? rel_err : -rel_err;
47 if (rel_err > EPS)
48 abort ();
49 }
50 }
51