File: a.12.1.c

package info (click to toggle)
gcc-riscv64-unknown-elf 8.3.0.2019.08%2Bdfsg-1
  • links: PTS, VCS
  • area: main
  • in suites: bullseye
  • size: 680,956 kB
  • sloc: ansic: 3,237,715; cpp: 896,882; ada: 772,854; f90: 144,254; asm: 68,788; makefile: 67,456; sh: 29,743; exp: 28,045; objc: 15,273; fortran: 11,885; python: 7,369; pascal: 5,375; awk: 3,725; perl: 2,872; yacc: 316; xml: 311; ml: 285; lex: 198; haskell: 122
file content (41 lines) | stat: -rw-r--r-- 719 bytes parent folder | download | duplicates (10) 
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);
  }
}