Viewing file:  vsub_f64.c (1.64 KB) -rw-r--r--
Select action/file-type:
 (+) |  (+) |  (+) | Code (+) | Session (+) |  (+) | SDB (+) |  (+) |  (+) |  (+) |  (+) |  (+) |
/* 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;
}
|
