1 /* { dg-do run } */
2 /* { dg-options "-mavx512f -O2" } */
3 /* { dg-require-effective-target avx512f } */
4
5 #include <math.h>
6 #include "avx512f-check.h"
7 #include "avx512f-helper.h"
8
9 static void
compute_vrsqrt14sd(double * s1,double * s2,double * r)10 compute_vrsqrt14sd (double *s1, double *s2, double *r)
11 {
12 r[0] = 1.0 / sqrt (s2[0]);
13 r[1] = s1[1];
14 }
15
16 static void
avx512f_test(void)17 avx512f_test (void)
18 {
19 union128d s1, s2, res1, res2, res3;
20 __mmask8 m = 0;
21 double res_ref[2];
22
23 s1.x = _mm_set_pd (-3.0, 111.111);
24 s2.x = _mm_set_pd (222.222, 4.0);
25
26 res1.x = _mm_rsqrt14_sd (s1.x, s2.x);
27
28 compute_vrsqrt14sd (s1.a, s2.a, res_ref);
29
30 if (check_fp_union128d (res1, res_ref))
31 abort ();
32
33 res2.x = _mm_set_pd (-4.0, DEFAULT_VALUE);
34 res2.x = _mm_mask_rsqrt14_sd(res2.x, m, s1.x, s2.x);
35
36 MASK_MERGE (d) (res_ref, m, 1);
37 if (checkVd (res2.a, res_ref, 2))
38 abort();
39
40 res3.x = _mm_maskz_rsqrt14_sd(m, s1.x, s2.x);
41
42 MASK_ZERO (d) (res_ref, m, 1);
43 if (checkVd (res3.a, res_ref, 2))
44 abort();
45 }
46