gcc/testsuite/gcc.target/i386/sse5-nmsubXX.c

   1 /* { dg-do run } */
   2 /* { dg-require-effective-target sse5 } */
   3 /* { dg-options "-O2 -msse5" } */
   4
   5 #include "sse5-check.h"
   6
   7 #include <bmmintrin.h>
   8 #include <string.h>
   9
  10 #define NUM 20
  11
  12 union
  13 {
  14   __m128 x[NUM];
  15   float f[NUM * 4];
  16   __m128d y[NUM];
  17   double d[NUM * 2];
  18 } dst, res, src1, src2, src3;
  19
  20 /* Note that in macc*,msub*,mnmacc* and mnsub* instructions, the intermdediate
  21    product is not rounded, only the addition is rounded. */
  22
  23 static void
  24 init_nmsubps ()
  25 {
  26   int i;
  27   for (i = 0; i < NUM * 4; i++)
  28     {
  29       src1.f[i] = i;
  30       src2.f[i] = i + 10;
  31       src3.f[i] = i + 20;
  32     }
  33 }
  34
  35 static void
  36 init_nmsubpd ()
  37 {
  38   int i;
  39   for (i = 0; i < NUM * 4; i++)
  40     {
  41       src1.d[i] = i;
  42       src2.d[i] = i + 10;
  43       src3.d[i] = i + 20;
  44     }
  45 }
  46
  47 static int
  48 check_nmsubps ()
  49 {
  50   int i, j, check_fails = 0;
  51   for (i = 0; i < NUM * 4; i = i + 4)
  52     for (j = 0; j < 4; j++)
  53       {
  54         res.f[i + j] = - (src1.f[i + j] * src2.f[i + j]) - src3.f[i + j];
  55         if (dst.f[i + j] != res.f[i + j])
  56           check_fails++;
  57       }
  58   return check_fails++;
  59 }
  60
  61 static int
  62 check_nmsubpd ()
  63 {
  64   int i, j, check_fails = 0;
  65   for (i = 0; i < NUM * 2; i = i + 2)
  66     for (j = 0; j < 2; j++)
  67       {
  68         res.d[i + j] = - (src1.d[i + j] * src2.d[i + j]) - src3.d[i + j];
  69         if (dst.d[i + j] != res.d[i + j])
  70           check_fails++;
  71       }
  72   return check_fails++;
  73 }
  74
  75
  76 static int
  77 check_nmsubss ()
  78 {
  79   int i, j, check_fails = 0;
  80   for (i = 0; i < NUM * 4; i = i + 4)
  81     {
  82       res.f[i] = - (src1.f[i] * src2.f[i]) - src3.f[i];
  83       if (dst.f[i] != res.f[i])
  84         check_fails++;
  85     }
  86   return check_fails++;
  87 }
  88
  89 static int
  90 check_nmsubsd ()
  91 {
  92   int i, j, check_fails = 0;
  93   for (i = 0; i < NUM * 2; i = i + 2)
  94     {
  95       res.d[i] = - (src1.d[i] * src2.d[i]) - src3.d[i];
  96       if (dst.d[i] != res.d[i])
  97         check_fails++;
  98     }
  99   return check_fails++;
 100 }
 101
 102 static void
 103 sse5_test (void)
 104 {
 105   int i;
 106
 107   /* Check nmsubps */
 108   init_nmsubps ();
 109
 110   for (i = 0; i < NUM; i++)
 111     dst.x[i] = _mm_nmsub_ps (src1.x[i], src2.x[i], src3.x[i]);
 112
 113   if (check_nmsubps (&dst.x[i], &src1.f[i * 4], &src2.f[i * 4], &src3.f[i * 4]))
 114     abort ();
 115
 116   /* check nmsubss */
 117   for (i = 0; i < NUM; i++)
 118     dst.x[i] = _mm_nmsub_ss (src1.x[i], src2.x[i], src3.x[i]);
 119
 120   if (check_nmsubss (&dst.x[i], &src1.f[i * 4], &src2.f[i * 4], &src3.f[i * 4]))
 121     abort ();
 122
 123   /* Check nmsubpd */
 124   init_nmsubpd ();
 125
 126   for (i = 0; i < NUM; i++)
 127     dst.y[i] = _mm_nmsub_pd (src1.y[i], src2.y[i], src3.y[i]);
 128
 129   if (check_nmsubpd (&dst.y[i], &src1.d[i * 2], &src2.d[i * 2], &src3.d[i * 2]))
 130     abort ();
 131
 132   /* Check nmsubps */
 133   for (i = 0; i < NUM; i++)
 134     dst.y[i] = _mm_nmsub_sd (src1.y[i], src2.y[i], src3.y[i]);
 135
 136   if (check_nmsubsd (&dst.y[i], &src1.d[i * 2], &src2.d[i * 2], &src3.d[i * 2]))
 137     abort ();
 138
 139 }