1  /* { dg-do run } */
       2  /* { dg-require-effective-target avx512er } */
       3  /* { dg-options "-O2 -ffast-math -ftree-vectorize -mavx512er" } */
       4  
       5  #include "avx512er-check.h"
       6  
       7  #define MAX 1000
       8  #define EPS 0.00001
       9  
      10  __attribute__ ((noinline, optimize (0)))
      11  void static
      12  compute_rcp_ref (float *a, float *b, float *r)
      13  {
      14    for (int i = 0; i < MAX; i++)
      15      r[i] = a[i] / b[i];
      16  }
      17  
      18  __attribute__ ((noinline))
      19  void static
      20  compute_rcp_exp (float *a, float *b, float *r)
      21  {
      22    for (int i = 0; i < MAX; i++)
      23      r[i] = a[i] / b[i];
      24  }
      25  
      26  void static
      27  avx512er_test (void)
      28  {
      29    float a[MAX];
      30    float b[MAX];
      31    float ref[MAX];
      32    float exp[MAX];
      33  
      34    for (int i = 0; i < MAX; i++)
      35      {
      36        a[i] = 179.345 - 6.5645 * i;
      37        b[i] = 8765.987 - 8.6756 * i;
      38      }
      39  
      40    compute_rcp_ref (a, b, ref);
      41    compute_rcp_exp (a, b, exp);
      42  
      43    for (int i = 0; i < MAX; i++)
      44      {
      45        float rel_err = (ref[i] - exp[i]) / ref[i];
      46        rel_err = rel_err > 0.0 ? rel_err : -rel_err;
      47        if (rel_err > EPS)
      48  	abort ();
      49      }
      50  }