Subversion Repositories openrisc

/* { 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);

  }

}