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  }