pi_sec_3.c

#include <stdio.h>
#include <stdlib.h>
#include <stdint.h>
#include <string.h>
#include <math.h>
#include <time.h>
#include <mpi.h>

#define SEED     921
#define NUM_ITER 1000000000

int main(int argc, char* argv[]) {
    int provided, rank, size;
    MPI_Init_thread(&argc, &argv, MPI_THREAD_SINGLE, &provided);

    double time = MPI_Wtime();

    int count = 0;
    double x, y, z, pi;
    MPI_Comm_rank(MPI_COMM_WORLD, &rank);
    MPI_Comm_size(MPI_COMM_WORLD, &size);
    srand(SEED * (rank + 1)); // Important: Multiply SEED by "rank" when you introduce MPI!
    
    // Calculate PI following a Monte Carlo method
    for (int iter = rank; iter < NUM_ITER; iter += size)
    {
        // Generate random (X,Y) points
        x = (double)random() / (double)RAND_MAX;
        y = (double)random() / (double)RAND_MAX;
        z = sqrt((x*x) + (y*y));
        
        // Check if point is in unit circle
        if (z <= 1.0)
        {
            count++;
        }
    }

    if (rank == 0) {
        int thread_count[size - 1];
        MPI_Request requests[size - 1];
        for (int i = 1; i < size; i++)
            MPI_Irecv(&thread_count[i - 1], 1, MPI_INT, i, i, MPI_COMM_WORLD, &requests[i - 1]);

        MPI_Waitall(size - 1, requests, MPI_STATUS_IGNORE);

        for (int i = 1; i < size; i++)
            count += thread_count[i - 1];

        // Estimate Pi and display the result
        pi = ((double)count / (double)NUM_ITER) * 4.0;
        time = MPI_Wtime() - time;
        printf("The result is %f, time: %lf sec\n", pi, time);
    }
    else {
        MPI_Request request;
        MPI_Isend(&count, 1, MPI_INT, 0, rank, MPI_COMM_WORLD, &request);
    }

    MPI_Finalize();

    return 0;
}