2 /* { dg-options "-mavx512f -O2 -std=c99" } */
3 /* { dg-require-effective-target avx512f } */
4 /* { dg-require-effective-target c99_runtime } */
6 #include "avx512f-check.h"
11 long long int int_val
;
16 get_norm_mant (double source
, int signctrl
, int interv
)
18 long long src
, sign
, exp
, fraction
;
20 union fp_int_t bin_conv
;
22 bin_conv
.fp_val
= source
;
23 src
= bin_conv
.int_val
;
24 sign
= (signctrl
& 0x1) ? 0 : (src
>> 63);
25 exp
= (src
& 0x7ff0000000000000) >> 52;
26 fraction
= (src
& 0xfffffffffffff);
29 return signbit (source
) ? -NAN
: NAN
;
30 if (source
== 0.0 || source
== -0.0 || isinf (source
))
31 return sign
? -1.0 : 1.0;
32 if (signbit (source
) && (signctrl
& 0x2))
34 if (!isnormal (source
))
36 src
= (src
& 0xfff7ffffffffffff);
38 while (!(src
& 0x8000000000000))
40 src
+= fraction
& 0x8000000000000;
41 fraction
= fraction
<< 1;
52 exp
= ((exp
- 0x3ff) & 0x1) ? 0x3fe : 0x3ff;
58 exp
= (fraction
& 0x8000000000000) ? 0x3fe : 0x3ff;
64 bin_conv
.int_val
= (sign
<< 63) | (exp
<< 52) | fraction
;
65 return bin_conv
.fp_val
;
69 compute_vgetmantsd (double *r
, double *s1
, double *s2
, int interv
,
72 r
[0] = get_norm_mant (s2
[0], signctrl
, interv
);
79 union128d res1
, src1
, src2
;
81 int interv
= _MM_MANT_NORM_p5_1
;
82 int signctrl
= _MM_MANT_SIGN_src
;
84 src1
.x
= _mm_set_pd (-3.0, 111.111);
85 src2
.x
= _mm_set_pd (222.222, -2.0);
87 res1
.x
= _mm_getmant_sd (src1
.x
, src2
.x
, interv
, signctrl
);
89 compute_vgetmantsd (res_ref
, src1
.a
, src2
.a
, interv
, signctrl
);
91 if (check_union128d (res1
, res_ref
))