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
      11  check_mm_fmsub_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_fmsub_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    if (check_union128d (e, d))
      25      abort ();
      26  }
      27  
      28  void
      29  check_mm_fmsub_ps (__m128 __A, __m128 __B, __m128 __C)
      30  {
      31    union128 a, b, c, e;
      32    a.x = __A;
      33    b.x = __B;
      34    c.x = __C;
      35    float d[4];
      36    int i;
      37    e.x = _mm_fmsub_ps (__A, __B, __C);
      38    for (i = 0; i < 4; i++)
      39      {
      40        d[i] = a.a[i] * b.a[i] - c.a[i];
      41      }
      42    if (check_union128 (e, d))
      43      abort ();
      44  }
      45  
      46  void
      47  check_mm_fmsub_sd (__m128d __A, __m128d __B, __m128d __C)
      48  {
      49    union128d a, b, c, e;
      50    a.x = __A;
      51    b.x = __B;
      52    c.x = __C;
      53    double d[2];
      54    int i;
      55    e.x = _mm_fmsub_sd (__A, __B, __C);
      56    for (i = 1; i < 2; i++)
      57      {
      58        d[i] = a.a[i];
      59      }
      60    d[0] = a.a[0] * b.a[0] - c.a[0];
      61    if (check_union128d (e, d))
      62      abort ();
      63  }
      64  
      65  void
      66  check_mm_fmsub_ss (__m128 __A, __m128 __B, __m128 __C)
      67  {
      68    union128 a, b, c, e;
      69    a.x = __A;
      70    b.x = __B;
      71    c.x = __C;
      72    float d[4];
      73    int i;
      74    e.x = _mm_fmsub_ss (__A, __B, __C);
      75    for (i = 1; i < 4; i++)
      76      {
      77        d[i] = a.a[i];
      78      }
      79    d[0] = a.a[0] * b.a[0] - c.a[0];
      80    if (check_union128 (e, d))
      81      abort ();
      82  }
      83  
      84  static void
      85  fma_test (void)
      86  {
      87    union128 a[3];
      88    union128d b[3];
      89    int i, j;
      90    for (i = 0; i < 3; i++)
      91      {
      92        for (j = 0; j < 4; j++)
      93  	a[i].a[j] = i * j + 3.5;
      94        for (j = 0; j < 2; j++)
      95  	b[i].a[j] = i * j + 3.5;
      96      }
      97    check_mm_fmsub_pd (b[0].x, b[1].x, b[2].x);
      98    check_mm_fmsub_sd (b[0].x, b[1].x, b[2].x);
      99    check_mm_fmsub_ps (a[0].x, a[1].x, a[2].x);
     100    check_mm_fmsub_ss (a[0].x, a[1].x, a[2].x);
     101  }