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
11void
12check_mm256_fmsub_pd (__m256d __A, __m256d __B, __m256d __C)
13{
14  union256d a, b, c, e;
15  a.x = __A;
16  b.x = __B;
17  c.x = __C;
18  double d[4];
19  int i;
20  e.x = _mm256_fmsub_pd (__A, __B, __C);
21  for (i = 0; i < 4; i++)
22    {
23      d[i] = a.a[i] * b.a[i] - c.a[i];
24    }
25  if (check_union256d (e, d))
26    abort ();
27}
28
29void
30check_mm256_fmsub_ps (__m256 __A, __m256 __B, __m256 __C)
31{
32  union256 a, b, c, e;
33  a.x = __A;
34  b.x = __B;
35  c.x = __C;
36  float d[8];
37  int i;
38  e.x = _mm256_fmsub_ps (__A, __B, __C);
39  for (i = 0; i < 8; i++)
40    {
41      d[i] = a.a[i] * b.a[i] - c.a[i];
42    }
43  if (check_union256 (e, d))
44    abort ();
45}
46
47static void
48fma_test (void)
49{
50  union256 c[3];
51  union256d d[3];
52  int i, j;
53  for (i = 0; i < 3; i++)
54    {
55      for (j = 0; j < 8; j++)
56	c[i].a[j] = i * j + 3.5;
57      for (j = 0; j < 4; j++)
58	d[i].a[j] = i * j + 3.5;
59    }
60  check_mm256_fmsub_pd (d[0].x, d[1].x, d[2].x);
61  check_mm256_fmsub_ps (c[0].x, c[1].x, c[2].x);
62}
63