| /* { dg-do run } */ |
| /* { dg-require-effective-target sse5 } */ |
| /* { dg-options "-O2 -msse5" } */ |
| |
| #include "sse5-check.h" |
| |
| #include <bmmintrin.h> |
| #include <string.h> |
| |
| #define NUM 20 |
| |
| union |
| { |
| __m128 x[NUM]; |
| float f[NUM * 4]; |
| __m128d y[NUM]; |
| double d[NUM * 2]; |
| } dst, res, src1, src2, src3; |
| |
| /* Note that in macc*,msub*,mnmacc* and mnsub* instructions, the intermdediate |
| product is not rounded, only the addition is rounded. */ |
| |
| static void |
| init_nmsubps () |
| { |
| int i; |
| for (i = 0; i < NUM * 4; i++) |
| { |
| src1.f[i] = i; |
| src2.f[i] = i + 10; |
| src3.f[i] = i + 20; |
| } |
| } |
| |
| static void |
| init_nmsubpd () |
| { |
| int i; |
| for (i = 0; i < NUM * 4; i++) |
| { |
| src1.d[i] = i; |
| src2.d[i] = i + 10; |
| src3.d[i] = i + 20; |
| } |
| } |
| |
| static int |
| check_nmsubps () |
| { |
| int i, j, check_fails = 0; |
| for (i = 0; i < NUM * 4; i = i + 4) |
| for (j = 0; j < 4; j++) |
| { |
| res.f[i + j] = - (src1.f[i + j] * src2.f[i + j]) - src3.f[i + j]; |
| if (dst.f[i + j] != res.f[i + j]) |
| check_fails++; |
| } |
| return check_fails++; |
| } |
| |
| static int |
| check_nmsubpd () |
| { |
| int i, j, check_fails = 0; |
| for (i = 0; i < NUM * 2; i = i + 2) |
| for (j = 0; j < 2; j++) |
| { |
| res.d[i + j] = - (src1.d[i + j] * src2.d[i + j]) - src3.d[i + j]; |
| if (dst.d[i + j] != res.d[i + j]) |
| check_fails++; |
| } |
| return check_fails++; |
| } |
| |
| |
| static int |
| check_nmsubss () |
| { |
| int i, j, check_fails = 0; |
| for (i = 0; i < NUM * 4; i = i + 4) |
| { |
| res.f[i] = - (src1.f[i] * src2.f[i]) - src3.f[i]; |
| if (dst.f[i] != res.f[i]) |
| check_fails++; |
| } |
| return check_fails++; |
| } |
| |
| static int |
| check_nmsubsd () |
| { |
| int i, j, check_fails = 0; |
| for (i = 0; i < NUM * 2; i = i + 2) |
| { |
| res.d[i] = - (src1.d[i] * src2.d[i]) - src3.d[i]; |
| if (dst.d[i] != res.d[i]) |
| check_fails++; |
| } |
| return check_fails++; |
| } |
| |
| static void |
| sse5_test (void) |
| { |
| int i; |
| |
| /* Check nmsubps */ |
| init_nmsubps (); |
| |
| for (i = 0; i < NUM; i++) |
| dst.x[i] = _mm_nmsub_ps (src1.x[i], src2.x[i], src3.x[i]); |
| |
| if (check_nmsubps (&dst.x[i], &src1.f[i * 4], &src2.f[i * 4], &src3.f[i * 4])) |
| abort (); |
| |
| /* check nmsubss */ |
| for (i = 0; i < NUM; i++) |
| dst.x[i] = _mm_nmsub_ss (src1.x[i], src2.x[i], src3.x[i]); |
| |
| if (check_nmsubss (&dst.x[i], &src1.f[i * 4], &src2.f[i * 4], &src3.f[i * 4])) |
| abort (); |
| |
| /* Check nmsubpd */ |
| init_nmsubpd (); |
| |
| for (i = 0; i < NUM; i++) |
| dst.y[i] = _mm_nmsub_pd (src1.y[i], src2.y[i], src3.y[i]); |
| |
| if (check_nmsubpd (&dst.y[i], &src1.d[i * 2], &src2.d[i * 2], &src3.d[i * 2])) |
| abort (); |
| |
| /* Check nmsubps */ |
| for (i = 0; i < NUM; i++) |
| dst.y[i] = _mm_nmsub_sd (src1.y[i], src2.y[i], src3.y[i]); |
| |
| if (check_nmsubsd (&dst.y[i], &src1.d[i * 2], &src2.d[i * 2], &src3.d[i * 2])) |
| abort (); |
| |
| } |