1  /* { dg-do run { target avx512fp16 } } */
       2  /* { dg-options "-O2 -mavx512fp16 -mfpmath=sse" } */
       3  
       4  #include <string.h>
       5  #include <stdlib.h>
       6  static void do_test (void);
       7  
       8  #define DO_TEST do_test
       9  #define AVX512FP16
      10  #include "avx512-check.h"
      11  #include "avx512fp16-11a.c"
      12  
      13  /* Get random float16 between -50.x to 50.x.  */
      14  _Float16
      15  get_float16_noround()
      16  {
      17    return ((int) (100.0 * rand ()/ (RAND_MAX + 1.0)) - 50)
      18      + 0.1f * (int) (10 * rand() / (RAND_MAX + 1.0));
      19  }
      20  
      21  static void
      22  do_test (void)
      23  {
      24    _Float16 x[32];
      25    _Float16 y[32];
      26    _Float16 res_add[32];
      27    _Float16 res_sub[32];
      28    _Float16 res_mul[32];
      29    _Float16 res_div[32];
      30    for (int i = 0 ; i != 32; i++)
      31      {
      32        x[i] = get_float16_noround ();
      33        y[i] = get_float16_noround ();
      34        if (y[i] == 0)
      35  	y[i] = 1.0f;
      36        res_add[i] = x[i] + y[i];
      37        res_sub[i] = x[i] - y[i];
      38        res_mul[i] = x[i] * y[i];
      39        res_div[i] = x[i] / y[i];
      40  
      41      }
      42  
      43    union512h u512 = { x[0], x[1], x[2], x[3], x[4], x[5], x[6], x[7],
      44        x[8], x[9], x[10], x[11], x[12], x[13], x[14], x[15],
      45        x[16], x[17], x[18], x[19], x[20], x[21], x[22], x[23],
      46        x[24], x[25], x[26], x[27], x[28], x[29], x[30], x[31] };
      47    union512h u512_1 = {y[0], y[1], y[2], y[3], y[4], y[5], y[6], y[7],
      48        y[8], y[9], y[10], y[11], y[12], y[13], y[14], y[15],
      49        y[16], y[17], y[18], y[19], y[20], y[21], y[22], y[23],
      50        y[24], y[25], y[26], y[27], y[28], y[29], y[30], y[31] };
      51  
      52    __m512h v512;
      53    union512h a512;
      54  
      55    memset (&v512, -1, sizeof (v512));
      56    v512 = vadd512 (u512.x, u512_1.x);
      57    a512.x = v512;
      58    if (check_union512h (a512, res_add))
      59      abort ();
      60    memset (&v512, -1, sizeof (v512));
      61    v512 = vsub512 (u512.x, u512_1.x);
      62    a512.x = v512;
      63    if (check_union512h (a512, res_sub))
      64      abort ();
      65    memset (&v512, -1, sizeof (v512));
      66    v512 = vmul512 (u512.x, u512_1.x);
      67    a512.x = v512;
      68    if (check_union512h (a512, res_mul))
      69      abort ();
      70    memset (&v512, -1, sizeof (v512));
      71    v512 = vdiv512 (u512.x, u512_1.x);
      72    a512.x = v512;
      73    if (check_union512h (a512, res_div))
      74      abort ();
      75  }