blob: fab948dd73686dfa4b37439c56c3195e85f98270 [file] [log] [blame]
/* { 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);
}
}