Viewing file:  task-reduction-7.C (3.1 KB) -rw-r--r--
Select action/file-type:
 (+) |  (+) |  (+) | Code (+) | Session (+) |  (+) | SDB (+) |  (+) |  (+) |  (+) |  (+) |  (+) |
typedef __SIZE_TYPE__ size_t;
extern "C" void abort ();
void
bar (int *a, int *b, int *c, int (*d)[2], int (*e)[4], int *f, int *g, size_t n)
{
#pragma omp task in_reduction (*: a[:n], b[3:n], c[n:n], d[0][:n], e[0][1:n], f[:n], g[1:n])
{
a[0] *= 12;
a[1] *= 13;
b[3] *= 14;
b[4] *= 15;
c[n] *= 16;
c[n + 1] *= 17;
d[0][0] *= 18;
d[0][1] *= 19;
e[0][1] *= 20;
e[0][2] *= 21;
f[0] *= 22;
f[1] *= 23;
g[1] *= 24;
g[2] *= 25;
}
}
void
foo (size_t n, void *x, void *y)
{
int a[n], b[n + 3], c[2 * n];
int (*d)[n] = (int (*)[n]) x;
int (*e)[n * 2] = (int (*)[n * 2]) y;
int fb[n], gb[n * 2];
int (&f)[n] = fb;
int (&g)[n * 2] = gb;
int i;
for (i = 0; i < n; i++)
{
a[i] = 1;
b[i + 3] = 1;
c[i + n] = 1;
d[0][i] = 1;
e[0][i + 1] = 1;
f[i] = 1;
g[i + 1] = 1;
}
#pragma omp taskgroup task_reduction (*: a, b[3:n], c[n:n], d[0][:n], e[0][1:n], f, g[1:n])
{
bar (a, b, c, (int (*)[2]) d, (int (*)[4]) e, &f[0], &g[0], n);
#pragma omp task in_reduction (*: a, b[3:n], c[n:n], d[0][:n], e[0][1:n], f, g[1:n])
{
a[0] *= 2;
a[1] *= 3;
b[3] *= 4;
b[4] *= 5;
c[n] *= 6;
c[n + 1] *= 7;
d[0][0] *= 8;
d[0][1] *= 9;
e[0][1] *= 10;
e[0][2] *= 11;
f[0] *= 12;
f[1] *= 13;
g[1] *= 14;
g[2] *= 15;
}
n = 0;
}
if (a[0] != 24 || a[1] != 39 || b[3] != 56 || b[4] != 75)
abort ();
if (c[2] != 96 || c[3] != 119 || d[0][0] != 144 || d[0][1] != 171)
abort ();
if (e[0][1] != 200 || e[0][2] != 231 || f[0] != 264 || f[1] != 299)
abort ();
if (g[1] != 336 || g[2] != 375)
abort ();
}
void
baz (size_t n, void *x, void *y)
{
int a[n], b[n + 3], c[2 * n];
int (*d)[n] = (int (*)[n]) x;
int (*e)[n * 2] = (int (*)[n * 2]) y;
int fb[n], gb[n * 2];
int i;
for (i = 0; i < n; i++)
{
a[i] = 1;
b[i + 3] = 1;
c[i + n] = 1;
d[0][i] = 1;
e[0][i + 1] = 1;
fb[i] = 1;
gb[i + 1] = 1;
}
#pragma omp parallel num_threads(2)
#pragma omp master
{
int (&f)[n] = fb;
int (&g)[n * 2] = gb;
#pragma omp taskgroup task_reduction (*: a, b[3:n], c[n:n], d[0][:n], e[0][1:n], f, g[1:n])
{
bar (a, b, c, (int (*)[2]) d, (int (*)[4]) e, &f[0], &g[0], n);
#pragma omp task in_reduction (*: a, b[3:n], c[n:n], d[0][:n], e[0][1:n], f, g[1:n])
{
a[0] *= 2;
a[1] *= 3;
b[3] *= 4;
b[4] *= 5;
c[n] *= 6;
c[n + 1] *= 7;
d[0][0] *= 8;
d[0][1] *= 9;
e[0][1] *= 10;
e[0][2] *= 11;
f[0] *= 12;
f[1] *= 13;
g[1] *= 14;
g[2] *= 15;
}
n = 0;
}
}
if (a[0] != 24 || a[1] != 39 || b[3] != 56 || b[4] != 75)
abort ();
if (c[2] != 96 || c[3] != 119 || d[0][0] != 144 || d[0][1] != 171)
abort ();
if (e[0][1] != 200 || e[0][2] != 231 || fb[0] != 264 || fb[1] != 299)
abort ();
if (gb[1] != 336 || gb[2] != 375)
abort ();
}
int
main ()
{
int d[2], e[4];
volatile int two;
two = 2;
#pragma omp parallel num_threads (2)
#pragma omp master
foo (two, (void *) d, (void *) e);
baz (two, (void *) d, (void *) e);
return 0;
}
|
