Python Tools for Reproducible Research on Hyperbolic Problems
by R. J. LeVeque. Submitted to the Computing in Science and Engineering (CiSE) special issue on reproducible research.

Abstract. Reproducible research in computational science is only possible if the computer codes used to generate published results are archived in a form that can later be used to regenerate the results and can be examined to determine details of the method used. Some difficulties in achieving this goal are discussed. A set of Python tools for facilitating reproducible research on finite volume methods for hyperbolic conservation laws using the Clawpack software are briefly surveyed and illustrated on a sample application.

Preprint: cise08draft.pdf

Python and Fortran code to reproduce the results in this paper:

• README file with links to browse

• cise08leveque.tar.gz Self-contained tar file. Unpack using:
  

  tar -zxvf cise08leveque.tar.gz

• Hints on getting python / matplotlib to work

The two-dimensional advection equation is solved on three different types of grids both with and without limiters and the accuracy compared.

• Grid 1: Cartesian
• Grid 2: Polar
• Grid 3: Quadrilateral grid from paper with D. Calhoun and C. Helzel

• Log-log plot of errors

bibtex entry:

@misc{rjl:cise08,
  author="R. J. LeVeque",
  title="Python Tools for Reproducible Research on Hyperbolic Problems",
  howpublished="Submitted to CiSE Magazine",
  year="2008"
}