1 /* { dg-do run } */
2 /* { dg-additional-options "-msse2" { target sse2_runtime } } */
3 /* { dg-additional-options "-mavx" { target avx_runtime } } */
4
5 #define N 128
6 #define M 16
7 #define EPS 0.0000000000000001
8 #define SAFELEN 16
9
10 #include <stdlib.h>
11
init(double * a,double * b,int n)12 void init(double *a, double *b, int n)
13 {
14 int i, s = -1;
15 for ( i = 0; i < n; i++ )
16 {
17 a[i] = i*i*s;
18 b[i] = a[i];
19 s = -s;
20 }
21 }
22
work(double * b,int n,int m)23 void work( double *b, int n, int m )
24 {
25 int i;
26 #pragma omp simd safelen(SAFELEN)
27 for (i = m; i < n; i++)
28 b[i] = b[i-m] - 1.0f;
29 }
30
work_ref(double * b,int n,int m)31 void work_ref( double *b, int n, int m )
32 {
33 int i;
34 for (i = m; i < n; i++)
35 b[i] = b[i-m] - 1.0f;
36 }
37
check(double * a,double * b)38 void check (double *a, double *b)
39 {
40 int i;
41 for (i = 0; i < N; i++)
42 if (a[i] - b[i] > EPS || b[i] - a[i] > EPS)
43 abort ();
44 }
45
main()46 int main ()
47 {
48 double b[N], b_ref[N];
49
50 init(b, b_ref, N);
51
52 work(b, N, M);
53 work(b_ref, N, M);
54
55 check(b, b_ref);
56
57 return 0;
58 }
59