103 lines
2.3 KiB
C
103 lines
2.3 KiB
C
// RUN: %libomp-compile-and-run
|
|
/*
|
|
* Threadprivate is tested in 2 ways:
|
|
* 1. The global variable declared as threadprivate should have
|
|
* local copy for each thread. Otherwise race condition and
|
|
* wrong result.
|
|
* 2. If the value of local copy is retained for the two adjacent
|
|
* parallel regions
|
|
*/
|
|
#include "omp_testsuite.h"
|
|
#include <stdlib.h>
|
|
#include <stdio.h>
|
|
|
|
static int sum0=0;
|
|
static int myvalue = 0;
|
|
|
|
#pragma omp threadprivate(sum0)
|
|
#pragma omp threadprivate(myvalue)
|
|
|
|
int test_omp_threadprivate()
|
|
{
|
|
int sum = 0;
|
|
int known_sum;
|
|
int i;
|
|
int iter;
|
|
int *data;
|
|
int size;
|
|
int num_failed = 0;
|
|
int my_random;
|
|
omp_set_dynamic(0);
|
|
|
|
#pragma omp parallel private(i)
|
|
{
|
|
sum0 = 0;
|
|
#pragma omp for
|
|
for (i = 1; i <= LOOPCOUNT; i++) {
|
|
sum0 = sum0 + i;
|
|
} /*end of for*/
|
|
#pragma omp critical
|
|
{
|
|
sum = sum + sum0;
|
|
} /*end of critical */
|
|
} /* end of parallel */
|
|
known_sum = (LOOPCOUNT * (LOOPCOUNT + 1)) / 2;
|
|
if (known_sum != sum ) {
|
|
fprintf (stderr, " known_sum = %d, sum = %d\n", known_sum, sum);
|
|
}
|
|
|
|
/* the next parallel region is just used to get the number of threads*/
|
|
omp_set_dynamic(0);
|
|
#pragma omp parallel
|
|
{
|
|
#pragma omp master
|
|
{
|
|
size=omp_get_num_threads();
|
|
data=(int*) malloc(size*sizeof(int));
|
|
}
|
|
}/* end parallel*/
|
|
|
|
srand(45);
|
|
for (iter = 0; iter < 100; iter++) {
|
|
my_random = rand(); /* random number generator is
|
|
called inside serial region*/
|
|
|
|
/* the first parallel region is used to initialiye myvalue
|
|
and the array with my_random+rank */
|
|
#pragma omp parallel
|
|
{
|
|
int rank;
|
|
rank = omp_get_thread_num ();
|
|
myvalue = data[rank] = my_random + rank;
|
|
}
|
|
|
|
/* the second parallel region verifies that the
|
|
value of "myvalue" is retained */
|
|
#pragma omp parallel reduction(+:num_failed)
|
|
{
|
|
int rank;
|
|
rank = omp_get_thread_num ();
|
|
num_failed = num_failed + (myvalue != data[rank]);
|
|
if(myvalue != data[rank]) {
|
|
fprintf (stderr, " myvalue = %d, data[rank]= %d\n",
|
|
myvalue, data[rank]);
|
|
}
|
|
}
|
|
}
|
|
free (data);
|
|
return (known_sum == sum) && !num_failed;
|
|
} /* end of check_threadprivate*/
|
|
|
|
int main()
|
|
{
|
|
int i;
|
|
int num_failed=0;
|
|
|
|
for(i = 0; i < REPETITIONS; i++) {
|
|
if(!test_omp_threadprivate()) {
|
|
num_failed++;
|
|
}
|
|
}
|
|
return num_failed;
|
|
}
|