Assignment 4: Name: Harshil Shah Reg No: 201070033 Subject: Parallel Computing

Name: Harshil Shah
Reg No: 201070033

Subject: Parallel Computing
ASSIGNMENT 4
AIM: Implement the following programs using openmp:
1. Fork join model
2. Producer consumer problem
3. Matrix multiplication
4. Find prime numbers
5. Largest element in an array
6. Calculate the value of pi
THEORY: OpenMP is a set of compiler directives as well as an API for programs

written in C, C++, or FORTRAN that provides support for parallel programming in
shared-memory environments. OpenMP identifies parallel regions as blocks of code that
may run in parallel. Application developers insert compiler directives into their code at
parallel regions, and these directives instruct the OpenMP run-time library to execute
the region in parallel.
When OpenMP encounters the directive
#pragma omp parallel
It creates as many threads which are processing cores in the system. Thus, for a dual-
core system, two threads are created, for a quad-core system, four are created; and so
forth. Then all the threads simultaneously execute the parallel region. When each
thread exits the parallel region, it is terminated. OpenMP provides several additional
directives for running code regions in parallel, including parallelizing loops.
In addition to providing directives for parallelization, OpenMP allows developers to
choose among several levels of parallelism. E.g., they can set the number of threads
manually. It also allows developers to identify whether data are shared between
threads or are private to a thread. OpenMP is available on several open-source and
commercial compilers for Linux, Windows, and Mac OS X systems
We must include the omp header file in the beginning of the program.
The syntax to compile an openmp code (parallel) is:
gcc -o obj -fopenmp file.c
OUTPUT :
1) Fork Join
// OpenMP program to print Hello World
#include <omp.h>
#include <stdio.h>
#include <stdlib.h>
int main()
int tid, nthreads;
// Beginning of parallel region
#pragma omp parallel private(nthreads, tid)
tid = omp_get_thread_num();
printf("Hello World... from thread = %d\n",
tid);
// Only master thread action
if(tid == 0){
nthreads = omp_get_num_threads();
printf("Number of threads = %d\n", nthreads);
}
}
// Ending of parallel region
2) Producer Consumer Problem
#include<stdio.h>
#include<omp.h>
int main()
{
int i=0; int x=0;

#pragma omp parallel shared(i)
{
if(omp_get_thread_num()==0)
{
printf("Master thread with Thread ID:%d\n", omp_get_thread_num());
printf("Since it is the producer thread It is adding some data to be consumed by other consumer threads\n");
i+=10; x=1;
}
else
{
while(x==0)
printf("Waiting for buffer to be filled. Thread ID:
%d\n",omp_get_thread_num());
#pragma critical
{
if(i>0){
printf("Data is consumed by Consumer with Thread ID: %d\n",omp_get_thread_num());
}
else {
i-=5;
printf("Could not find any data for thread ID: %d\n",omp_get_thread_num());

}}
}
}
}
3) Matrix multiplication
#include<stdio.h>
#include<omp.h>
#include<stdlib.h>
int main(){
int i,j,k,m,n,p;
printf("Enter the number of rows in Matrix 1:");
scanf("%d",&m);
int *matrixA[m];
printf("Enter the number of columns in Matrix 1:");
scanf("%d",&n);
for(i=0;i<m;i++){
matrixA[i] = (int *)malloc(n*sizeof(int));
}
printf("<--Now Input the values for matrix 1 row-wise-->\n");
for(i=0;i<m;i++){
for(j=0;j<n;j++){
scanf("%d",&matrixA[i][j]);
}
}
printf("Enter the number of columns in Matrix 2:");

scanf("%d",&p);
int *matrixB[n];
for(i=0;i<n;i++){
matrixB[i] = (int *)malloc(p*sizeof(int));
}
printf("<--Now Input the values for matrix 2 row-wise-->\n");
for(i=0;i<n;i++){
for(j=0;j<p;j++){
scanf("%d",&matrixB[i][j]);
}
}
int matrixC[m][p];
#pragma omp parallel private(i,j,k) shared(matrixA,matrixB,matrixC)
{
#pragma omp for schedule(static)

for (i=0; i<m; i=i+1){
for (j=0; j<p; j=j+1){
matrixC[i][j] = 0;
for (k=0; k<n; k=k+1){
matrixC[i][j]=(matrixC[i][j])+((matrixA[i][k])*(matrixB[k][j]));
}
}
}
}
printf("The output after Matrix Multiplication is: \n"); for(i=0;i<m;i++){
for(j=0;j<p;j++)
printf("%d \t",matrixC[i][j]);
printf("\n");
}
return 0;
4) Find the prime number

#include<stdio.h>
#include<omp.h>
#include<stdlib.h>
int IsPrime(int number)

{
int i;
for (i = 2; i < number; i++) {
if (number % i == 0 && i != number) return 0;
}
return 1;
}
int main()
{
int noOfThreads,valueN,indexCount=0,arrayVal[10000],tempValue;
printf("Enter the Number of threads: ");
scanf("%d",&noOfThreads);
printf("Enter the value of N: ");

scanf("%d",&valueN);
omp_set_num_threads(noOfThreads);
#pragma omp parallel for reduction(+:indexCount)
for(tempValue=2;tempValue<=valueN;tempValue++)
{
if(IsPrime(tempValue)){ arrayVal[indexCount] = tempValue;
indexCount++;
}
}
printf("Number of prime numbers between 2 and %d: %d\n",valueN,indexCount);
return 0;
}
5 ) Largest element in array

#include<stdio.h>
#include<omp.h>
#include<stdlib.h>
int main()
int
numberOfElements,currentMax=-1,iIterator,arrayInput[10000]; printf("Enter the

Number of Elements: "); scanf("%d",&numberOfElements);
for(iIterator=0;iIterator<numberOfElements;iIterator++){
scanf("%d",&arrayInput[iIterator]);
#pragma omp parallel for shared(currentMax)
for(iIterator=0;iIterator<numberOfElements;iIterator++){
#pragma omp critical if(arrayInput[iIterator] >
currentMax){
currentMax = arrayInput[iIterator];
printf("The Maximum Element is:
%d\n",currentMax); return 0;
}
6) Pi calculation
#include<stdio.h>
#include<omp.h>
#include<stdlib.h>
int main(){
int num_steps=10000,i;
double aux,pi,step = 1.0/(double) num_steps,x=0.0,sum = 0.0;
#pragma omp parallel private(i,x,aux) shared(sum) {
#pragma omp for schedule(static)
for (i=0; i<num_steps; i=i+1){
x=(i+0.5)*step;
aux=4.0/(1.0+x*x);
#pragma omp critical sum = sum + aux;
}
}
pi=step*sum;
printf("The Value of PI is %lf\n",pi);
return 0;
CONCLUSION: In this experiment, we became more familiar with programming in

openmp andbecame more accustomed to parallel programming in C. Hence the
objectives of this experiment are satisfied.

Assignment 4: Name: Harshil Shah Reg No: 201070033 Subject: Parallel Computing

Uploaded by

Copyright:

Available Formats

Assignment 4: Name: Harshil Shah Reg No: 201070033 Subject: Parallel Computing

Uploaded by

Document Information

Original Description:

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Assignment 4: Name: Harshil Shah Reg No: 201070033 Subject: Parallel Computing

Uploaded by

Copyright:

Available Formats

Name: Harshil Shah

Reg No: 201070033

AIM: Implement the following programs using openmp:

1. Fork join model

2. Producer consumer problem

4. Find prime numbers

5. Largest element in an array

6. Calculate the value of pi

THEORY: OpenMP is a set of compiler directives as well as an API for programs

When OpenMP encounters the directive

#pragma omp parallel

The syntax to compile an openmp code (parallel) is:

gcc -o obj -fopenmp file.c

int tid, nthreads;

// Beginning of parallel region

#pragma omp parallel private(nthreads, tid)

printf("Hello World... from thread = %d\n",

// Only master thread action

printf("Number of threads = %d\n", nthreads);

int i=0; int x=0;

printf("Could not find any data for thread ID: %d\n",omp_get_thread_num());

printf("Enter the number of columns in Matrix 2:");

#pragma omp for schedule(static)

4) Find the prime number

int IsPrime(int number)

printf("Enter the value of N: ");

5 ) Largest element in array

numberOfElements,currentMax=-1,iIterator,arrayInput[10000]; printf("Enter the

#pragma omp parallel for shared(currentMax)

#pragma omp critical if(arrayInput[iIterator] >

printf("The Maximum Element is:

CONCLUSION: In this experiment, we became more familiar with programming in

You might also like