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
|
/* { dg-do compile } */
#include <stdio.h>
extern float average (float, float, float);
void
a12 (float *x, float *xold, int n, float tol)
{
int c, i, toobig;
float error, y;
c = 0;
#pragma omp parallel
{
do
{
#pragma omp for private(i)
for (i = 1; i < n - 1; ++i)
{
xold[i] = x[i];
}
#pragma omp single
{
toobig = 0;
}
#pragma omp for private(i,y,error) reduction(+:toobig)
for (i = 1; i < n - 1; ++i)
{
y = x[i];
x[i] = average (xold[i - 1], x[i], xold[i + 1]);
error = y - x[i];
if (error > tol || error < -tol)
++toobig;
}
#pragma omp master
{
++c;
printf ("iteration %d, toobig=%d\n", c, toobig);
}
}
while (toobig > 0);
}
}
|
