/* Test vsub works correctly.  */
/* { dg-do run } */
/* { dg-options "--save-temps" } */

#include <arm_neon.h>

#define FLT_EPSILON __FLT_EPSILON__
#define DBL_EPSILON __DBL_EPSILON__

#define TESTA0 1
#define TESTA1 0.2223
#define TESTA2 0
#define TESTA3 -0.76544
/* 2^54, double has 53 significand bits
   according to Double-precision floating-point format.  */
#define TESTA4 18014398509481984
#define TESTA5 2.0

#define TESTB0 0.66667
#define TESTB1 2
#define TESTB2 0
#define TESTB3 -2
#define TESTB4 1.0
#define TESTB5 (1.0 / TESTA4)

#define ANSW0 0.33333
#define ANSW1 -1.7777
#define ANSW2 0
#define ANSW3 1.23456
#define ANSW4 TESTA4
#define ANSW5 2.0

extern void abort (void);

#define EPSILON __DBL_EPSILON__
#define ISNAN(a) __builtin_isnan (a)
/* FP_equals is implemented like this to execute subtraction
   exectly once during a single test run.  */
#define FP_equals(a, b, epsilon)		\
(						\
 ((a) == (b))					\
 || (ISNAN (a) && ISNAN (b))			\
 || (((a > b) && (a < (b + epsilon)))		\
     || ((b > a) && (b < (a + epsilon))))	\
)

int
test_vsub_f64 ()
{
  float64x1_t a;
  float64x1_t b;
  float64x1_t c;

  a = TESTA0;
  b = TESTB0;
  c = ANSW0;

  a = vsub_f64 (a, b);
  if (!FP_equals (a, c, EPSILON))
    return 1;

  a = TESTA1;
  b = TESTB1;
  c = ANSW1;

  a = vsub_f64 (a, b);
  if (!FP_equals (a, c, EPSILON))
    return 1;

  a = TESTA2;
  b = TESTB2;
  c = ANSW2;

  a = vsub_f64 (a, b);
  if (!FP_equals (a, c, EPSILON))
    return 1;

  a = TESTA3;
  b = TESTB3;
  c = ANSW3;

  a = vsub_f64 (a, b);
  if (!FP_equals (a, c, EPSILON))
    return 1;

  a = TESTA4;
  b = TESTB4;
  c = ANSW4;

  a = vsub_f64 (a, b);
  if (!FP_equals (a, c, EPSILON))
    return 1;

  a = TESTA5;
  b = TESTB5;
  c = ANSW5;

  a = vsub_f64 (a, b);
  if (!FP_equals (a, c, EPSILON))
    return 1;

  return 0;
}

/* { dg-final { scan-assembler-times "fsub\\td\[0-9\]+, d\[0-9\]+, d\[0-9\]+" 6 } } */

int
main (int argc, char **argv)
{
  if (test_vsub_f64 ())
    abort ();

  return 0;
}

/* { dg-final { cleanup-saved-temps } } */