/* { dg-do run } */
|
/* { dg-do run } */
|
/* { dg-options "-O2 -ffast-math -msse -mfpmath=sse -mrecip" } */
|
/* { dg-options "-O2 -ffast-math -msse -mfpmath=sse -mrecip" } */
|
/* { dg-require-effective-target sse } */
|
/* { dg-require-effective-target sse } */
|
|
|
#include "sse-check.h"
|
#include "sse-check.h"
|
|
|
extern float sqrtf (float);
|
extern float sqrtf (float);
|
extern float fabsf (float);
|
extern float fabsf (float);
|
|
|
#define N 8
|
#define N 8
|
|
|
static void
|
static void
|
sse_test (void)
|
sse_test (void)
|
{
|
{
|
float a[N] = { 0.f, 18.f, 108.f, 324.f, 720.f, 1944.f, 3087.f, 5832.f };
|
float a[N] = { 0.f, 18.f, 108.f, 324.f, 720.f, 1944.f, 3087.f, 5832.f };
|
float b[N] = { 1.f, 2.f, 3.f, 4.f, 5.f, 6.f, 7.f, 8.f };
|
float b[N] = { 1.f, 2.f, 3.f, 4.f, 5.f, 6.f, 7.f, 8.f };
|
float r[N];
|
float r[N];
|
|
|
float rc[N] = { 0.f, 3.f, 6.f, 9.f, 12.f, 18.f, 21.f, 27.f };
|
float rc[N] = { 0.f, 3.f, 6.f, 9.f, 12.f, 18.f, 21.f, 27.f };
|
|
|
int i;
|
int i;
|
|
|
for (i = 0; i < N; i++)
|
for (i = 0; i < N; i++)
|
{
|
{
|
r[i] = sqrtf (a[i] / b[i]);
|
r[i] = sqrtf (a[i] / b[i]);
|
}
|
}
|
|
|
/* check results: */
|
/* check results: */
|
for (i = 0; i < N; i++)
|
for (i = 0; i < N; i++)
|
{
|
{
|
if (fabsf (r[i] - rc[i]) > 0.0001)
|
if (fabsf (r[i] - rc[i]) > 0.0001)
|
abort();
|
abort();
|
}
|
}
|
}
|
}
|
|
|