Version 10 (modified by steve, 12 years ago) (diff)


INSTALLING anuga_parallel

First you should install the most uptodate version of the code. Follow the instructions to install Anuga on Ubuntu.


Well first you need to get the anuga_parallel code. You can get this from our svn repository with userid anonymous (blank password)

The location is

(By the way, the most recent version of the development code of anuga is available at )

Setup your PYTHONPATH to point to location of the source directory

For instance I have the following line in my .bashrc file

export PYTHONPATH=/home/steve/anuga/anuga_core/source


Now you need to install MPI on your system. OPENMPI and MPICH2 are supported by pypar (see below) so both should be ok.

Make sure mpi works. You should be able to run a program in parallel. Try something as simple as

mpirun -np 4 pwd

should produce the output of pwd 4 times.


We use pypar as the interface between mpi and python. The most recent version of PYPAR is available from

(There is an old version on sourceforge don't use that)

Install pypar following the instructions in the download. You should be able use the standard command

python install

or maybe

sudo python install

Make sure the pypar examples work


In the anuga_parallel directory there is a subdirectory pymetis.

Follow the instructions in README to install. Essentially just run make.

If you have a 64 bit machine run


From the pymetis directory, test using, ie



Should now be ready to run some parallel anuga code. Go back to the anuga_parallel directory and run

Hopefully that all works.

Example program


First just run it as a sequential program, via


Then try a parallel run using a command like

mpirun -np 4 python

That should run on 4 processors

You should look at the code in

Essentially a fairly standard example, with the extra command

domain = distribute(domain)

which sets up all the parallel stuff.

Also for efficiency reasons we only setup the original full sequential mesh on processor 0, hence the statement

if myid == 0:
     domain = create_domain_from_file(mesh_filename)
     domain.set_quantity('stage', Set_Stage(x0, x1, 2.0))
     domain = None

The output will be an sww file associated to each processor.

There is a script anuga/utilities/ which provides a function to merge sww files into one sww file for viewing with the anuga viewer.

Suppose your parallel code produced 3 sww files, domain_P3_0.sww domain_P3_1.sww and domain_P3_2.sww

The base name would be "domain" and the number of processors would be 3. To stitch these 3 files together either run the as a script with the command

python /dir/to/anuga/utilities/ -f domain -np 3

or add the following command at the end of your simulation script