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 } } } } } */