1  /* { dg-do run } */
       2  /* { dg-options "-O3 -mpower8-vector -Wno-psabi" } */
       3  /* { dg-require-effective-target p8vector_hw } */
       4  
       5  #ifndef CHECK_H
       6  #define CHECK_H "sse3-check.h"
       7  #endif
       8  
       9  #include CHECK_H
      10  
      11  #ifndef TEST
      12  #define TEST sse3_test_addsubpd_1
      13  #endif
      14  
      15  #define NO_WARN_X86_INTRINSICS 1
      16  #include <pmmintrin.h>
      17  
      18  static void
      19  sse3_test_addsubpd (double *i1, double *i2, double *r)
      20  {
      21    __m128d t1 = _mm_loadu_pd (i1);
      22    __m128d t2 = _mm_loadu_pd (i2);
      23  
      24    t1 = _mm_addsub_pd (t1, t2);
      25  
      26    _mm_storeu_pd (r, t1);
      27  }
      28  
      29  static void
      30  sse3_test_addsubpd_subsume (double *i1, double *i2, double *r)
      31  {
      32    __m128d t1 = _mm_load_pd (i1);
      33    __m128d t2 = _mm_load_pd (i2);
      34  
      35    t1 = _mm_addsub_pd (t1, t2);
      36  
      37    _mm_storeu_pd (r, t1);
      38  }
      39  
      40  static int
      41  chk_pd (double *v1, double *v2)
      42  {
      43    int i;
      44    int n_fails = 0;
      45  
      46    for (i = 0; i < 2; i++)
      47      if (v1[i] != v2[i])
      48        n_fails += 1;
      49  
      50    return n_fails;
      51  }
      52  
      53  static double p1[2] __attribute__ ((aligned(16)));
      54  static double p2[2] __attribute__ ((aligned(16)));
      55  static double p3[2];
      56  static double ck[2];
      57  
      58  double vals[] =
      59    {
      60      100.0,  200.0, 300.0, 400.0, 5.0, -1.0, .345, -21.5,
      61      1100.0, 0.235, 321.3, 53.40, 0.3, 10.0, 42.0, 32.52,
      62      32.6,   123.3, 1.234, 2.156, 0.1, 3.25, 4.75, 32.44,
      63      12.16,  52.34, 64.12, 71.13, -.1, 2.30, 5.12, 3.785,
      64      541.3,  321.4, 231.4, 531.4, 71., 321., 231., -531.,
      65      23.45,  23.45, 23.45, 23.45, 23.45, 23.45, 23.45, 23.45,
      66      23.45,  -1.43, -6.74, 6.345, -20.1, -20.1, -40.1, -40.1,
      67      1.234,  2.345, 3.456, 4.567, 5.678, 6.789, 7.891, 8.912,
      68      -9.32,  -8.41, -7.50, -6.59, -5.68, -4.77, -3.86, -2.95,
      69      9.32,  8.41, 7.50, 6.59, -5.68, -4.77, -3.86, -2.95
      70    };
      71  
      72  static
      73  void
      74  TEST (void)
      75  {
      76    int i;
      77    int fail = 0;
      78  
      79    for (i = 0; i < sizeof (vals) / sizeof (vals[0]); i += 4)
      80      {
      81        p1[0] = vals[i+0];
      82        p1[1] = vals[i+1];
      83  
      84        p2[0] = vals[i+2];
      85        p2[1] = vals[i+3];
      86  
      87        ck[0] = p1[0] - p2[0];
      88        ck[1] = p1[1] + p2[1];
      89  
      90        sse3_test_addsubpd (p1, p2, p3);
      91  
      92        fail += chk_pd (ck, p3);
      93  
      94        sse3_test_addsubpd_subsume (p1, p2, p3);
      95  
      96        fail += chk_pd (ck, p3);
      97      }
      98  
      99    if (fail != 0)
     100      abort ();
     101  }