omp_set_num_threads(p) sets the number of threads to be created to p.

The following statement creates p (set by omp_set_num_threads) parallel processes and assigns n/p iterations to each process in serial order. So, first n/p interations will be done by first process.

#pragma omp parallel for

for(i=0; i<n; i++)

{

     root[i] = sqrt(numbers[i]);

}

The complete code is given below.

#include “stdafx.h”

#include <omp.h>

#include <math.h>

#include <iostream>

#include <stdio.h>

#include <conio.h>

#define n 10

using namespace std; /*used for cout*/

int main()

{

int i;

float numbers[n], root[n];

for(i=0; i<n; i++)

{

numbers[i] = (i+1);

}

omp_set_num_threads(n); /*creates n parallel threads*/

#pragma omp parallel for /*runs the statement following in parallel*/

for(i=0; i<n; i++)

{

root[i] = sqrt(numbers[i]);

}

for(i=0; i<n; i++)

{

cout<<numbers[i]<<” “<<root[i]<<endl;

}

getch(); /*just to pause the output */

}