gcc/testsuite/gcc.dg/vect/pr79920.c

   1 /* { dg-additional-options "-O3 -fno-fast-math" } */
   2
   3 #include "tree-vect.h"
   4
   5 double __attribute__((noinline,noclone))
   6 compute_integral (double w_1[18])
   7 {
   8   double A = 0;
   9   double t33[2][6] = {{0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
  10         {0.0, 0.0, 0.0, 0.0, 0.0, 0.0}};
  11   double t43[2] = {0.0, 0.0};
  12   double t31[2][2] = {{1.0, 1.0}, {1.0, 1.0}};
  13   double t32[2][3] = {{0.0, 0.0, 1.0}, {0.0, 0.0, 1.0}};
  14
  15   for (int ip_1 = 0; ip_1 < 2; ++ip_1)
  16     {
  17 #pragma GCC unroll 0
  18       for (int i_0 = 0; i_0 < 6; ++i_0)
  19         t33[ip_1][i_0] = ((w_1[i_0*3] * t32[ip_1][0])
  20                           + (w_1[i_0*3+2] * t32[ip_1][2]));
  21       t43[ip_1] = 2.0;
  22     }
  23   for (int i_0 = 0; i_0 < 6; ++i_0)
  24     A += t43[1]*t33[1][i_0];
  25   return A;
  26 }
  27
  28 int main()
  29 {
  30   check_vect ();
  31
  32   double w_1[18] = {0., 1.0, 1.0,
  33       0., 1.0, 1.0,
  34       0., 1.0, 1.0,
  35       0., 1.0, 1.0,
  36       0., 1.0, 1.0,
  37       0., 1.0, 1.0};
  38   double A = compute_integral(w_1);
  39   if (A != 12.0)
  40     __builtin_abort ();
  41   return 0;
  42 }
  43
  44 /* { dg-final { scan-tree-dump-times {using an in-order \(fold-left\) reduction} 1 "vect" { target vect_double } } } */
  45 /* { dg-final { scan-tree-dump-times "vectorized 2 loops" 1 "vect" { target { vect_double && { vect_perm && vect_hw_misalign } } } } } */