summaryrefslogtreecommitdiff
path: root/gcc/testsuite/gcc.dg/gomp/appendix-a/a.12.1.c
blob: fab948dd73686dfa4b37439c56c3195e85f98270 (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
/* { 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);
  }
}