File: e.55.1.c

package info (click to toggle)
gcc-avr 1%3A5.4.0%2BAtmel3.6.2-1
  • links: PTS, VCS
  • area: main
  • in suites: bullseye
  • size: 589,872 kB
  • sloc: ansic: 2,775,581; ada: 756,757; cpp: 723,977; f90: 117,673; asm: 66,898; makefile: 62,755; xml: 44,466; sh: 29,549; exp: 23,315; objc: 15,216; fortran: 10,901; pascal: 4,185; python: 4,093; perl: 2,969; awk: 2,811; ml: 2,385; cs: 879; yacc: 316; lex: 198; haskell: 112; lisp: 8
file content (69 lines) | stat: -rw-r--r-- 1,045 bytes parent folder | download | duplicates (6) 
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
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
 
/* { dg-do run } */
/* { dg-additional-options "-DCHUNKSZ=5000" { target { ! run_expensive_tests } } } */
/* { dg-additional-options "-DCHUNKSZ=1000" { target run_expensive_tests } } */

#include <stdlib.h>

#define EPS 0.00001
#define N 100000

float Y[N];
float Z[N];

#pragma omp declare target
float F (float a)
{
  return -a;
}
#pragma omp end declare target

void pipedF_ref ()
{
  int i;
  for (i = 0; i < N; i++)
    Y[i] = F (Y[i]);
}

void pipedF ()
{
  int i, C;
  for (C = 0; C < N; C += CHUNKSZ)
    {
      #pragma omp task
	#pragma omp target map(Z[C:CHUNKSZ])
	  #pragma omp parallel for
	    for (i = C; i < C + CHUNKSZ; i++)
	      Z[i] = F (Z[i]);
    }
  #pragma omp taskwait
}

void init ()
{
  int i;
  for (i = 0; i < N; i++)
    Y[i] = Z[i] = 0.1 * i;
}

void check ()
{
  int i;
  for (i = 0; i < N; i++)
    {
      float err = (Z[i] == 0.0) ? Y[i] : (Y[i] - Z[i]) / Z[i];
      if (((err > 0) ? err : -err) > EPS)
	abort ();
    }
}

int main ()
{
  init ();

  pipedF_ref ();
  pipedF ();

  check ();

  return 0;
}