131
|
1 /* { dg-additional-options "-O3 -fno-fast-math" } */
|
111
|
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 for (int i_0 = 0; i_0 < 6; ++i_0)
|
|
18 t33[ip_1][i_0] = ((w_1[i_0*3] * t32[ip_1][0])
|
|
19 + (w_1[i_0*3+2] * t32[ip_1][2]));
|
|
20 t43[ip_1] = 2.0;
|
|
21 }
|
|
22 for (int i_0 = 0; i_0 < 6; ++i_0)
|
|
23 A += t43[1]*t33[1][i_0];
|
|
24 return A;
|
|
25 }
|
|
26
|
|
27 int main()
|
|
28 {
|
|
29 check_vect ();
|
|
30
|
|
31 double w_1[18] = {0., 1.0, 1.0,
|
|
32 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 double A = compute_integral(w_1);
|
|
38 if (A != 12.0)
|
|
39 __builtin_abort ();
|
|
40 return 0;
|
|
41 }
|
|
42
|
131
|
43 /* { dg-final { scan-tree-dump-times {using an in-order \(fold-left\) reduction} 1 "vect" { target vect_double } } } */
|
|
44 /* { dg-final { scan-tree-dump-times "vectorized 2 loops" 1 "vect" { target { vect_double && { vect_perm && vect_hw_misalign } } } } } */
|