aboutsummaryrefslogtreecommitdiff
path: root/rc3/runtime/test/worksharing/for/omp_for_private.c
blob: 1f537b90c8f212e9c1249a4a9ab7ccd343fe1fa9 (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
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
// RUN: %libomp-compile-and-run
#include <stdio.h>
#include <math.h>
#include "omp_testsuite.h"

/* Utility function do spend some time in a loop */
static void do_some_work()
{
  int i;
  double sum = 0;
  for(i = 0; i < 1000; i++){
  sum += sqrt ((double) i);
  }
}

int sum1;
#pragma omp threadprivate(sum1)

int test_omp_for_private()
{
  int sum = 0;
  int sum0;
  int known_sum;

  sum0 = 0;  /* setting (global) sum0 = 0 */

  #pragma omp parallel
  {
    sum1 = 0;  /* setting sum1 in each thread to 0 */
    {  /* begin of orphaned block */
      int i;
      #pragma omp for private(sum0) schedule(static,1)
      for (i = 1; i <= LOOPCOUNT; i++) {
        sum0 = sum1;
        #pragma omp flush
        sum0 = sum0 + i;
        do_some_work ();
        #pragma omp flush
        sum1 = sum0;
      }
    }  /* end of orphaned block */

    #pragma omp critical
    {
      sum = sum + sum1;
    }  /*end of critical*/
  }  /* end of parallel*/
  known_sum = (LOOPCOUNT * (LOOPCOUNT + 1)) / 2;
  return (known_sum == sum);
}

int main()
{
  int i;
  int num_failed=0;

  for(i = 0; i < REPETITIONS; i++) {
    if(!test_omp_for_private()) {
      num_failed++;
    }
  }
  return num_failed;
}