// implementation of threading to estimate the value of pi
// using the instructions given in Q 4.24
#include <stdio.h>
#include <pthread.h>
#include <math.h>
#include <stdlib.h>
#include <limits.h>
#define SEED 122
#define data_size 8000000
unsigned int seed = 122;
// defining the total number of points from which random points are selected
//int N = 8000000;
//global variable for number of points inside unit circle
int totalPts = 0;
//Random number generator with linear congruential generator
double RandUint (long i)
{
seed = seed * 2896512345 + 427738;
return seed / (double)UINT_MAX;
}
//function to find the number of random points inside the unit circle
void * Count (void *X)
{
/* initialize random numbers */
for (int i = 0; i < N; i++)
{
double x = (double)RandUint (i);
double y = (double)RandUint (i);
// checking if it is a unit circle
if (((x * x) + (y * y)) <= 1)
{
totalPts++;
}
}
return NULL;
}
//function to estimate pi using pthreads
void EstimateParallel ()
{
srand (SEED);
// defining the thread variable as pthread_t datatype (4 bytes)
pthread_t thread;
// creating a pthread which computes the number of points inside the unit circle
pthread_create (&thread, NULL, &Count, NULL);
//waiting for the termination of the pthread
pthread_join (thread, NULL);
// estimating the value of pi using the formula given in the question
double estPiValue = (4.0 * totalPts) / N;
printf ("No of trials: %d \nNo of points inside unit circle: %d \nEst value of
pi:%lf\n",N,totalPts,estPiValue);
int main ()
{
/ defining the total number of points from which random points are selected
int N = 8000000;
// calculating the time taken by the function estimate()
clock_t t;
t = clock ();
EstimateParallel ();
t = clock () - t;
double timeTaken = ((double)t) / CLOCKS_PER_SEC; // in seconds
printf ("EstimateParallel() took %f seconds to execute \n", timeTaken);
char data[data_size];
FILE * fPtr;
fPtr = fopen("inp.txt", "w");
if(fPtr == NULL)
{
/* File not created hence exit */
printf("Unable to create file.\n");
exit(EXIT_FAILURE);
}
fgets(N, DATA_SIZE, stdin);
fputs(N, fPtr);
fgets(totalPts, DATA_SIZE, stdin);
fputs(totalPts, fPtr);
fgets(estPivalue, DATA_SIZE, stdin);
fputs(estPivalue, fPtr);
fclose(fPtr);
return 0;
}
