Scribd
Upload
14 views

Editedcode

Uploaded by

cs20b1023
Copyright
© © All Rights Reserved
Available Formats
Download as TXT, PDF, TXT or read online on Scribd
14 views

Editedcode

Uploaded by

cs20b1023
Copyright
© © All Rights Reserved
Available Formats
Download as TXT, PDF, TXT or read online on Scribd
You are on page 1/ 3

// 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;
}

You might also like