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