{"id":3483,"date":"2024-05-23T10:16:55","date_gmt":"2024-05-23T05:16:55","guid":{"rendered":"https:\/\/afzalbadshah.com\/?p=3483"},"modified":"2024-05-23T10:19:14","modified_gmt":"2024-05-23T05:19:14","slug":"monte-carlo-simulation-mpi4py","status":"publish","type":"post","link":"https:\/\/afzalbadshah.com\/index.php\/2024\/05\/23\/monte-carlo-simulation-mpi4py\/","title":{"rendered":"Monte Carlo Simulation: MPI4Py"},"content":{"rendered":"\n

Monte Carlo simulations are a statistical technique that allows for solving problems through random sampling. They are widely used in various fields such as physics, finance, and engineering to understand the impact of risk and uncertainty in prediction and forecasting models. The core idea is to use randomness to solve problems that might be deterministic in nature. You can visit the detailed tutorial here. <\/a><\/p>\n\n\n\n

For example, to estimate the value of Pi ((\\pi)), we can use the Monte Carlo method. We generate random points in a unit square and count how many fall within a quarter-circle inscribed within the square. The ratio of the points inside the quarter-circle to the total number of points approximates (\\pi\/4).<\/p>\n\n\n\n

Detailed Program<\/h3>\n\n\n\n

Below is a Python program used mpi4py<\/code> for parallel processing in a Monte Carlo simulation to estimate the value of Pi.<\/p>\n\n\n\n

from mpi4py import MPI\nimport numpy as np\nimport random\n\ndef monte_carlo_pi(num_samples):\n    count = 0\n    for _ in range(num_samples):\n        x, y = random.random(), random.random()\n        if x**2 + y**2 <= 1:\n            count += 1\n    return count\n\ndef main():\n    comm = MPI.COMM_WORLD\n    rank = comm.Get_rank()\n    size = comm.Get_size()\n\n    num_samples_per_proc = 1000000  # Adjust as needed\n    total_samples = num_samples_per_proc * size\n\n    # Each process runs its own Monte Carlo simulation\n    local_count = monte_carlo_pi(num_samples_per_proc)\n\n    # Gather all local counts to the root process\n    total_count = comm.reduce(local_count, op=MPI.SUM, root=0)\n\n    if rank == 0:\n        # Root process calculates the final estimate of Pi\n        pi_estimate = (4.0 * total_count) \/ total_samples\n        print(f\"Estimated Pi: {pi_estimate}\")\n\nif __name__ == \"__main__\":\n    main()<\/code><\/pre>\n\n\n\n
Line-by-Line Explanation<\/h3>\n\n\n\n
from mpi4py import MPI<\/code><\/pre>\n\n\n\n
This line imports the MPI<\/code> module from the mpi4py<\/code> package, which provides support for parallel processing.<\/p>\n\n\n\n
import numpy as np\nimport random<\/code><\/pre>\n\n\n\n
These lines import the numpy<\/code> and random<\/code> modules. numpy<\/code> is a powerful library for numerical computations, although it’s not explicitly used in this example, and random<\/code> provides functions for generating random numbers.<\/p>\n\n\n\n
def monte_carlo_pi(num_samples):<\/code><\/pre>\n\n\n\n
This line defines a function named monte_carlo_pi<\/code> that takes one parameter num_samples<\/code>, representing the number of random points to generate.<\/p>\n\n\n\n
    count = 0<\/code><\/pre>\n\n\n\n
Initializes a counter count<\/code> to zero. This will be used to count the number of points that fall inside the quarter-circle.<\/p>\n\n\n\n
    for _ in range(num_samples):<\/code><\/pre>\n\n\n\n
Starts a loop that iterates num_samples<\/code> times to generate random points.<\/p>\n\n\n\n
        x, y = random.random(), random.random()<\/code><\/pre>\n\n\n\n
Generates two random floating-point numbers x<\/code> and y<\/code> between 0 and 1, representing the coordinates of a point.<\/p>\n\n\n\n
        if x**2 + y**2 <= 1:<\/code><\/pre>\n\n\n\n
Checks if the point (x, y) lies inside the quarter-circle by verifying if the sum of the squares of x<\/code> and y<\/code> is less than or equal to 1.<\/p>\n\n\n\n
            count += 1<\/code><\/pre>\n\n\n\n
If the point lies inside the quarter-circle, increments the counter count<\/code> by 1.<\/p>\n\n\n\n
    return count<\/code><\/pre>\n\n\n\n
Returns the total count of points that lie inside the quarter-circle.<\/p>\n\n\n\n
def main():<\/code><\/pre>\n\n\n\n
Defines the main function which will be executed when the script is run.<\/p>\n\n\n\n
    comm = MPI.COMM_WORLD<\/code><\/pre>\n\n\n\n
Initializes the MPI communication world, allowing communication between different processes.<\/p>\n\n\n\n
    rank = comm.Get_rank()<\/code><\/pre>\n\n\n\n
Gets the rank (ID) of the current process. The rank is used to identify each process uniquely.<\/p>\n\n\n\n
    size = comm.Get_size()<\/code><\/pre>\n\n\n\n
Gets the total number of processes running in the MPI world.<\/p>\n\n\n\n
    num_samples_per_proc = 1000000  # Adjust as needed<\/code><\/pre>\n\n\n\n
Defines the number of samples each process will generate. This can be adjusted based on the desired accuracy and computational resources.<\/p>\n\n\n\n
    total_samples = num_samples_per_proc * size<\/code><\/pre>\n\n\n\n
Calculates the total number of samples by multiplying the number of samples per process by the total number of processes.<\/p>\n\n\n\n
    local_count = monte_carlo_pi(num_samples_per_proc)<\/code><\/pre>\n\n\n\n
Each process runs the monte_carlo_pi<\/code> function independently to simulate its portion of the samples and stores the result in local_count<\/code>.<\/p>\n\n\n\n
    total_count = comm.reduce(local_count, op=MPI.SUM, root=0)<\/code><\/pre>\n\n\n\n
Gathers all local counts to the root process (rank 0) and sums them up using the reduce operation.<\/p>\n\n\n\n
    if rank == 0:<\/code><\/pre>\n\n\n\n
Checks if the current process is the root process.<\/p>\n\n\n\n
        pi_estimate = (4.0 * total_count) \/ total_samples<\/code><\/pre>\n\n\n\n
The root process calculates the final estimate of Pi using the total count of points inside the quarter-circle and the total number of samples.<\/p>\n\n\n\n
        print(f\"Estimated Pi: {pi_estimate}\")<\/code><\/pre>\n\n\n\n
Prints the estimated value of Pi.<\/p>\n\n\n\n
if __name__ == \"__main__\":\n    main()<\/code><\/pre>\n\n\n\n
Checks if the script is being run directly and, if so, calls the main<\/code> function to start the program.<\/p>\n\n\n\n
Running the Program<\/h3>\n\n\n\n
To run this program with MPI, use the following command, specifying the number of processes you want to use (e.g., 4):<\/p>\n\n\n\n
mpiexec -n 4 python your_script_name.py<\/code><\/pre>\n\n\n\n
This command will start the script with 4 parallel processes, each performing part of the Monte Carlo simulation and collaborating to estimate the value of Pi.<\/p>\n\n\n\n
\n