gcc/testsuite/gcc.target/i386/avx512f-vgetmantps-2.c

   1 /* { dg-do run } */
   2 /* { dg-options "-O2 -mavx512f -std=c99 -mfpmath=sse" } */
   3 /* { dg-require-effective-target avx512f } */
   4 /* { dg-require-effective-target c99_runtime } */
   5
   6 #define AVX512F
   7
   8 #include "avx512f-helper.h"
   9
  10 #define SIZE (AVX512F_LEN / 32)
  11 #include "avx512f-mask-type.h"
  12 #include <math.h>
  13
  14 #ifndef GET_NORM_MANT
  15 #define GET_NORM_MANT
  16
  17 union fp_int_t
  18 {
  19   int int_val;
  20   float fp_val;
  21 };
  22
  23 float
  24 get_norm_mant (float source, int signctrl, int interv)
  25 {
  26   int src, sign, exp, fraction;
  27   union fp_int_t bin_conv;
  28
  29   bin_conv.fp_val = source;
  30   src = bin_conv.int_val;
  31   sign = (signctrl & 0x1) ? 0 : (src >> 31);
  32   exp = (src & 0x7f800000) >> 23;
  33   fraction = (src & 0x7fffff);
  34
  35   if (isnan (source))
  36     return signbit (source) ? -NAN : NAN;
  37   if (source == 0.0 || source == -0.0 || isinf (source))
  38     return sign ? -1.0 : 1.0;
  39   if (signbit (source) && (signctrl & 0x2))
  40     return -NAN;
  41   if (!isnormal (source))
  42     {
  43       src = (src & 0xffbfffff);
  44       exp = 0x7f;
  45       while (!(src & 0x400000))
  46         {
  47           src += fraction & 0x400000;
  48           fraction = fraction << 1;
  49           exp--;
  50         }
  51     }
  52
  53   switch (interv)
  54     {
  55     case 0:
  56       exp = 0x7f;
  57       break;
  58     case 1:
  59       exp = ((exp - 0x7f) & 0x1) ? 0x7e : 0x7f;
  60       break;
  61     case 2:
  62       exp = 0x7e;
  63       break;
  64     case 3:
  65       exp = (fraction & 0x400000) ? 0x7e : 0x7f;
  66       break;
  67     default:
  68       abort ();
  69     }
  70
  71   bin_conv.int_val = (sign << 31) | (exp << 23) | fraction;
  72
  73   return bin_conv.fp_val;
  74 }
  75 #endif
  76
  77 void static
  78 CALC (float *r, float *s, int interv, int signctrl)
  79 {
  80   int i;
  81   for (i = 0; i < SIZE; i++)
  82     {
  83       r[i] = get_norm_mant (s[i], signctrl, interv);
  84     }
  85 }
  86
  87 void
  88 TEST (void)
  89 {
  90   int i, sign;
  91   UNION_TYPE (AVX512F_LEN,) res1, res2, res3, src;
  92   MASK_TYPE mask = MASK_VALUE;
  93   float res_ref[SIZE];
  94   int interv = _MM_MANT_NORM_p5_1;
  95   int signctrl = _MM_MANT_SIGN_src;
  96
  97   sign = -1;
  98   for (i = 0; i < SIZE; i++)
  99     {
 100       src.a[i] = 34.67 * i * sign;
 101       sign = sign * -1;
 102     }
 103   for (i = 0; i < SIZE; i++)
 104     res2.a[i] = DEFAULT_VALUE;
 105
 106   res1.x = INTRINSIC (_getmant_ps) (src.x, interv, signctrl);
 107   res2.x =
 108     INTRINSIC (_mask_getmant_ps) (res2.x, mask, src.x, interv,
 109                                   signctrl);
 110   res3.x =
 111     INTRINSIC (_maskz_getmant_ps) (mask, src.x, interv, signctrl);
 112
 113   CALC (res_ref, src.a, interv, signctrl);
 114
 115   if (UNION_CHECK (AVX512F_LEN,) (res1, res_ref))
 116     abort ();
 117
 118   MASK_MERGE ()(res_ref, mask, SIZE);
 119   if (UNION_CHECK (AVX512F_LEN,) (res2, res_ref))
 120     abort ();
 121
 122   MASK_ZERO ()(res_ref, mask, SIZE);
 123   if (UNION_CHECK (AVX512F_LEN,) (res3, res_ref))
 124     abort ();
 125 }