blob: e067a81e5620b7203ba4ed0732f57294f573768f [file] [log] [blame]
/* { dg-do run } */
/* { dg-require-effective-target avx512er } */
/* { dg-options "-O2 -ffast-math -ftree-vectorize -mavx512er" } */
#include <math.h>
#include "avx512er-check.h"
#define MAX 1000
#define EPS 0.00001
__attribute__ ((noinline, optimize (1)))
void static
compute_sqrt_ref (float *a, float *r)
{
for (int i = 0; i < MAX; i++)
r[i] = sqrtf (a[i]);
}
__attribute__ ((noinline))
void static
compute_sqrt_exp (float *a, float *r)
{
for (int i = 0; i < MAX; i++)
r[i] = sqrtf (a[i]);
}
void static
avx512er_test (void)
{
float in[MAX];
float ref[MAX];
float exp[MAX];
for (int i = 0; i < MAX; i++)
in[i] = 8765.987 - 8.6756 * i;
compute_sqrt_ref (in, ref);
compute_sqrt_exp (in, exp);
for (int i = 0; i < MAX; i++)
{
float rel_err = (ref[i] - exp[i]) / ref[i];
rel_err = rel_err > 0.0 ? rel_err : -rel_err;
if (rel_err > EPS)
abort ();
}
}