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

#include <arm_neon.h>

#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))))	\
)

#define TEST(N)					\
int						\
test_vsub_f64_##N ()				\
{						\
  float64x1_t a = { TESTA##N };			\
  float64x1_t b = { TESTB##N };			\
  float64x1_t c = { ANSW##N };			\
						\
  a = vsub_f64 (a, b);				\
  return !FP_equals (a[0], c[0], EPSILON);	\
}

TEST (0)
TEST (1)
TEST (2)
TEST (3)
TEST (4)
TEST (5)

/* { 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_0 ())
    abort ();
  if (test_vsub_f64_1 ())
    abort ();
  if (test_vsub_f64_2 ())
    abort ();
  if (test_vsub_f64_3 ())
    abort ();
  if (test_vsub_f64_4 ())
    abort ();
  if (test_vsub_f64_5 ())
    abort ();

  return 0;
}