1 /* { dg-do run } */
2 /* { dg-require-effective-target fma } */
3 /* { dg-options "-O2 -mfma" } */
4
5 #include "fma-check.h"
6
7 #include <x86intrin.h>
8 #include "m256-check.h"
9
10 void
check_mm_fmadd_pd(__m128d __A,__m128d __B,__m128d __C)11 check_mm_fmadd_pd (__m128d __A, __m128d __B, __m128d __C)
12 {
13 union128d a, b, c, e;
14 a.x = __A;
15 b.x = __B;
16 c.x = __C;
17 double d[2];
18 int i;
19 e.x = _mm_fmadd_pd (__A, __B, __C);
20 for (i = 0; i < 2; i++)
21 {
22 d[i] = a.a[i] * b.a[i] + c.a[i];
23 }
24
25 if (check_union128d (e, d))
26 abort ();
27 }
28
29 void
check_mm_fmadd_ps(__m128 __A,__m128 __B,__m128 __C)30 check_mm_fmadd_ps (__m128 __A, __m128 __B, __m128 __C)
31 {
32 union128 a, b, c, e;
33 a.x = __A;
34 b.x = __B;
35 c.x = __C;
36 float d[4];
37 int i;
38 e.x = _mm_fmadd_ps (__A, __B, __C);
39 for (i = 0; i < 4; i++)
40 {
41 d[i] = a.a[i] * b.a[i] + c.a[i];
42 }
43 if (check_union128 (e, d))
44 abort ();
45 }
46
47 void
check_mm_fmadd_sd(__m128d __A,__m128d __B,__m128d __C)48 check_mm_fmadd_sd (__m128d __A, __m128d __B, __m128d __C)
49 {
50 union128d a, b, c, e;
51 a.x = __A;
52 b.x = __B;
53 c.x = __C;
54 double d[2];
55 int i;
56 e.x = _mm_fmadd_sd (__A, __B, __C);
57 for (i = 1; i < 2; i++)
58 {
59 d[i] = a.a[i];
60 }
61 d[0] = a.a[0] * b.a[0] + c.a[0];
62 if (check_union128d (e, d))
63 abort ();
64 }
65
66 void
check_mm_fmadd_ss(__m128 __A,__m128 __B,__m128 __C)67 check_mm_fmadd_ss (__m128 __A, __m128 __B, __m128 __C)
68 {
69 union128 a, b, c, e;
70 a.x = __A;
71 b.x = __B;
72 c.x = __C;
73 float d[4];
74 int i;
75 e.x = _mm_fmadd_ss (__A, __B, __C);
76 for (i = 1; i < 4; i++)
77 {
78 d[i] = a.a[i];
79 }
80 d[0] = a.a[0] * b.a[0] + c.a[0];
81 if (check_union128 (e, d))
82 abort ();
83 }
84
85 static void
fma_test(void)86 fma_test (void)
87 {
88 union128 a[3];
89 union128d b[3];
90 int i, j;
91 for (i = 0; i < 3; i++)
92 {
93 for (j = 0; j < 4; j++)
94 a[i].a[j] = i * j + 3.5;
95 for (j = 0; j < 2; j++)
96 b[i].a[j] = i * j + 3.5;
97 }
98 check_mm_fmadd_pd (b[0].x, b[1].x, b[2].x);
99 check_mm_fmadd_sd (b[0].x, b[1].x, b[2].x);
100 check_mm_fmadd_ps (a[0].x, a[1].x, a[2].x);
101 check_mm_fmadd_ss (a[0].x, a[1].x, a[2].x);
102 }
103