See figures.py.

To parallelize, start ipcluster engines and then run

import parallel
parallel.go_parallel()