9.2. Forest Fire Simulation Example

Note

You should use the VNC connection to the head node of your Raspberry Pi cluster, because this example produces graphical output.

This example was ported to Python by Libby Shoop (Macalester College), from the original Shodor foundation C code example.

According to the above Shodor page:

“You can run a single instance of a forest fire simulation in which a forest is modeled as an \(N \times N\) grid of trees. One tree starts to smolder, and each iteration nearby trees have some chance of catching fire. The model follows the following rules:”

The main input parameters for the model are:

The simulation starts with the tree in the center of the grid smoldering, and all other trees alive.

The main outputs for the single fire model are:

One single instance of this model can produce a different result each time it is run because of the randomness of the probability of unburnt trees catching fire.

9.2.1. Overview of Code

You can access the forest fire simulation from the head node of the CSinParallel Raspberry Pi cluster by navigating to the fire directory under Exemplars/mpi4py-examples:

cd CSinParallel/Exemplars/mpi4py-examples/fire

The folder called fire has a simple simulation of the spread of a wildfire in a forest.

  • fire_sequential_once.py runs a single model of a fire starting from one tree and burning until it is out.

  • fire_sequential_simulate.py runs many of the above single model cases for many ‘trials’ and averages the results to get a reasonable overall behavior of a fire at different thresholds of probability that fire will spread from one tree to another.

  • fire_mpi_simulate.py does the same simulation as the above sequential version, but uses parallelism to split the number of trials among processes.

9.2.2. Patterns used

The mpi version of the simulation uses these patterns from the mpi4py patternlets:

  • SPMD, Single program, multiple data

  • Decompose the work using parallel loop of equal chunks

  • Broadcast

  • Message passing

Study the file fire_mpi_simulate.py in an editor to find these.

9.2.3. Running the Code (Sequential Version)

The code in fire_sequential_once.py creates a visualization serially that looks like the following, with \(N\) = 25 and the probability threshold = 0.4.

../_images/single_forest_fire_simulation.png

Output like this can be obtained by running the single model like this:

python fire_sequential_once.py 25 0.4

Note

The resulting plot can be closed by clicking the X icon in the upper left of the window that pops up. The terminal prompt will not return until you do.

Try running this several times at the same 0.4 threshold. Then try varying the threshold from 0.2 to 1.0 by 0.1 increments.

Each time the code is run, the result could be different. In addition, even if we ran several trials, the resulting percent of trees burned and number of iterations before the fire burned out on average would be different, depending on the input probability threshold. Because of this, a more realistic simulation requires that many instances of the above single simulation be run in this way:

  • Keep the size of the grid of trees the same.

  • Start with a low probability threshold.

  • Run a given number of trials at a fixed probability threshold.

  • Repeat for another probability threshold, some increment larger than the previous one, until the threshold is 1.0.

This simulation of multiple trials at a range of different probability thresholds has a known interesting output, which can be graphed as follows:

many simulations over a range of thresholds

In this case, we ran 20 trials on a single Raspberry Pi 3B, with the probability threshold starting at 0.1 and incrementing by 0.1. We did this running the code file fire_sequential_simulate.py on a cluster head node like this:

python fire_sequential_simulate.py 10 0.1 20

As the size of the grid changes and the probability points increase, this curve will look roughly the same, although it should get smoother as the number of trials increases and the increment value is smaller. But these more accurate simulations take a long time to run.

9.2.4. The parallel MPI version

The desired outcome of the parallel version is to also produce a plot of average percent burns as a function of probability of spreading, as quickly and as accurately as possible. This should take into account that the probability of the fire spreading will affect not only how long it takes for the fire to burn out but also the number of iterations required to reached an accurate representation of the average.

If we put more processes to work on the problem, we should be able to complete a more accurate simulation in less time than the sequential version. Even the same problem as above produced the same results running on 4 processes on different nodes of a cluster in almost 1/4 of the time. Its output looks like this:

many simulations over a range of thresholds

This can be run like this:

mpirun -np 4 -hostfile ~/hostfile python fire_mpi_simulate.py 10 0.1 20

The parallelization happens by splitting up the number of trials to be run among the processes. Each process completes the range of probabilities for its portion of the trials, sending the results back to the conductor process.

9.2.4.1. Try some other cases to observe how it scales

Ideally, as you double the number of workers on the same problem, the time should be cut in half. This is called strong scalability. But there is some overhead from the message passing, so we don’t often see perfect strong scalability.

Try running these tests (be patient, since 40 trials is a fair amount for these little Raspberry Pis):

-np

tree row size

probability increment

number of trials

running time

3

20

0.1

40

6

20

0.1

40

12

20

0.1

40

What do you observe about the time as you double the number of workers?

When does the message passing cause the most overhead, which adds to the running time?

Try some other cases of your own design.

You have attempted of activities on this page