gcc/testsuite/gcc.target/i386/sse4_1-dpps-2.c

   1 /* { dg-do run } */
   2 /* { dg-require-effective-target sse4 } */
   3 /* { dg-options "-O2 -msse4.1" } */
   4
   5 #ifndef CHECK_H
   6 #define CHECK_H "sse4_1-check.h"
   7 #endif
   8
   9 #ifndef TEST
  10 #define TEST sse4_1_test
  11 #endif
  12
  13 #include CHECK_H
  14
  15 #include <smmintrin.h>
  16
  17 #include <string.h>
  18
  19 #define lmskN  0x00
  20 #define lmsk0  0x01
  21 #define lmsk1  0x02
  22 #define lmsk2  0x04
  23 #define lmsk3  0x08
  24 #define lmsk01 0x03
  25 #define lmsk02 0x05
  26 #define lmsk03 0x09
  27 #define lmsk12 0x06
  28 #define lmsk13 0x0A
  29 #define lmsk23 0x0C
  30 #define lmskA  0x0F
  31
  32 #define hmskN  0x00
  33 #define hmskA  0xF0
  34 #define hmsk0  0x10
  35 #define hmsk1  0x20
  36 #define hmsk2  0x40
  37 #define hmsk3  0x80
  38 #define hmsk01 0x30
  39 #define hmsk02 0x50
  40 #define hmsk03 0x90
  41 #define hmsk12 0x60
  42 #define hmsk13 0xA0
  43 #define hmsk23 0xC0
  44
  45 #ifndef HIMASK
  46 #define HIMASK hmskA
  47 #endif
  48
  49 #ifndef LOMASK
  50 #define LOMASK lmskA
  51 #endif
  52
  53 static void
  54 TEST (void)
  55 {
  56   union
  57     {
  58       __m128 x;
  59       float f[4];
  60     } val1[16], val2[16], res[16], chk[16];
  61   int i,j;
  62   float tmp;
  63
  64   for (i = 0; i < 16; i++)
  65     {
  66       val1[i].f[0] = 2.;
  67       val1[i].f[1] = 3.;
  68       val1[i].f[2] = 4.;
  69       val1[i].f[3] = 5.;
  70
  71       val2[i].f[0] = 10.;
  72       val2[i].f[1] = 100.;
  73       val2[i].f[2] = 1000.;
  74       val2[i].f[3] = 10000.;
  75
  76       tmp = 0.;
  77       for (j = 0; j < 4; j++)
  78         if ((HIMASK & (0x10 << j)))
  79           tmp += val1[i].f [j] * val2[i].f [j];
  80
  81       for (j = 0; j < 4; j++)
  82         if ((LOMASK & (1 << j)))
  83           chk[i].f[j] = tmp;
  84     }
  85
  86    for (i = 0; i < 16; i++)
  87      {
  88        res[i].x = _mm_dp_ps (val1[i].x, val2[i].x, HIMASK | LOMASK);
  89        if (memcmp (&res[i], &chk[i], sizeof (chk[i])))
  90          abort ();
  91      }
  92 }