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

#include <arm_neon.h>

#define TESTA0 0.33333
#define TESTA1 -1.7777
#define TESTA2 0
#define TESTA3 1.23456
/* 2^54, double has 53 significand bits
   according to Double-precision floating-point format.  */
#define TESTA4 18014398509481984
#define TESTA5 (1.0 / TESTA4)

#define TESTB0 0.66667
#define TESTB1 2
#define TESTB2 0
#define TESTB3 -2
#define TESTB4 1.0
#define TESTB5 2.0

#define ANSW0 1
#define ANSW1 0.2223
#define ANSW2 0
#define ANSW3 -0.76544
#define ANSW4 TESTA4
#define ANSW5 2.0

extern void abort (void);

#define EPSILON __DBL_EPSILON__
#define ABS(a) __builtin_fabs (a)
#define ISNAN(a) __builtin_isnan (a)
#define FP_equals(a, b, epsilon)			\
  (							\
   ((a) == (b))						\
    || (ISNAN (a) && ISNAN (b))				\
    || (ABS (a - b) < epsilon)				\
   )

#define TEST(N)					\
int						\
test_vadd_f64_##N ()				\
{						\
  float64x1_t a = { TESTA##N };			\
  float64x1_t b = { TESTB##N };			\
  float64x1_t c = { ANSW##N };			\
						\
  a = vadd_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 "fadd\\td\[0-9\]+, d\[0-9\]+, d\[0-9\]+" 6 } } */

int
main (int argc, char **argv)
{
  if (test_vadd_f64_0 ())
    abort ();
  if (test_vadd_f64_1 ())
    abort ();
  if (test_vadd_f64_2 ())
    abort ();
  if (test_vadd_f64_3 ())
    abort ();
  if (test_vadd_f64_4 ())
    abort ();
  if (test_vadd_f64_5 ())
    abort ();

  return 0;
}