CSC/ECE 506 Spring 2015/37 mr

1. Example 1
2. Example 2

Parallelize the code below using loop distribution

---

for (i=0; i<n; i++) { 
S1:  a[i] = b[i-1] * a[i+1];
S2:  b[i] = b[i] * 3;
S3:  c[i] = .2 * (c[i-1] * c[i]);
} 

After loop distribution

---

for(i=0; i < n; i++) { S1: a[i] = b[i-1] * a[i+1]; } for(i=0; i < n; i++) S3: c[i] = .2* c[i-1] * c[i]; for(i=0; i < n; i++) S2: b[i] = b[i] * 3;

---

Now write the parallelized code in openMP using C++

---

1. pragma omp parallel shared(a,b,c)private(i) {

#pragma omp sections nowait { #pragma omp section for(i=0; i < n; i++) a[i] = b[i-1] * a[i+1]; #pragma omp section for(i=0; i<n; i++) c[i] = .2* c[i-1]*c[i]; #pragma omp section for(i=0; i < n; i++) b[i] = b[i] * 3; }//end omp sections } //end omp parallel