blob: b2e1f1feef81a298a343297f6ab86e057735f8fc (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
|
/* { dg-do run } */
/* { dg-additional-options "-DCHUNKSZ=5000" { target { ! run_expensive_tests } } } */
/* { dg-additional-options "-DCHUNKSZ=1000" { target run_expensive_tests } } */
#include <stdlib.h>
#define EPS 0.00001
#define N 100000
float Y[N];
float Z[N];
#pragma omp declare target
float F (float a)
{
return -a;
}
#pragma omp end declare target
void pipedF_ref ()
{
int i;
for (i = 0; i < N; i++)
Y[i] = F (Y[i]);
}
void pipedF ()
{
int i, C;
for (C = 0; C < N; C += CHUNKSZ)
{
#pragma omp task
#pragma omp target map(Z[C:CHUNKSZ])
#pragma omp parallel for
for (i = C; i < C + CHUNKSZ; i++)
Z[i] = F (Z[i]);
}
#pragma omp taskwait
}
void init ()
{
int i;
for (i = 0; i < N; i++)
Y[i] = Z[i] = 0.1 * i;
}
void check ()
{
int i;
for (i = 0; i < N; i++)
{
float err = (Z[i] == 0.0) ? Y[i] : (Y[i] - Z[i]) / Z[i];
if (((err > 0) ? err : -err) > EPS)
abort ();
}
}
int main ()
{
init ();
pipedF_ref ();
pipedF ();
check ();
return 0;
}
|
