1 /******************************************************************************
2 * OpenMP Example - Combined Parallel Loop Work-sharing - C/C++ Version
3 * FILE: omp_workshare4.c
4 * DESCRIPTION:
5 *   This is a corrected version of the omp_workshare3.c example. Corrections
6 *   include removing all statements between the parallel for construct and
7 *   the actual for loop, and introducing logic to preserve the ability to
8 *   query a thread's id and print it from inside the for loop.
9 * SOURCE: Blaise Barney  5/99
10 * LAST REVISED: 03/03/2002
11 ******************************************************************************/
12 
13 #include <omp.h>
14 #include <stdio.h>
15 #define N       50
16 #define CHUNKSIZE   5
17 
18 main ()  {
19 
20 int i, chunk, tid;
21 float a[N], b[N], c[N];
22 char first_time;
23 
24 /* Some initializations */
25 for (i=0; i < N; i++)
26   a[i] = b[i] = i * 1.0;
27 chunk = CHUNKSIZE;
28 first_time = 'y';
29 
30 #pragma omp parallel for     \
31   shared(a,b,c,chunk)            \
32   private(i,tid)             \
33   schedule(static,chunk)     \
34   firstprivate(first_time)
35 
36   for (i=0; i < N; i++)
37     {
38     if (first_time == 'y')
39       {
40       tid = omp_get_thread_num();
41       first_time = 'n';
42       }
43     c[i] = a[i] + b[i];
44     printf("tid= %d i= %d c[i]= %f\n", tid, i, c[i]);
45     }
46 
47   return 0;
48 }
49