Portable, Extensible Toolkit for Scientific Computation

infobox software
name = PETSc
latest_release_version = 2.3.3
latest_release_date = 23 May 2007
operating_system = Linux, Unix, Mac OS X, Windows
license = own compatible with [GNU General Public Licenc [e|GPL (version 2)]
language = C (main language),C++, FORTRAN
genre = Scientific simulation software
website = http://www-unix.mcs.anl.gov/petsc/

The Portable, Extensible Toolkit for Scientific Computation (PETSc, pronounced PET-see; the S is silent), is a suite of data structures and routines for the scalable (parallel) solution of scientific applications modeled by partial differential equations. It employs the Message Passing Interface (MPI) standard for all message-passing communication. The current version of PETSc is 2.3.3; released May, 23, 2007.

PETSc is intended for use in large-scale application projects, many ongoing computational science projects are built around the PETSc libraries. Its careful design allows advanced users to have detailed control over the solution process. PETSc includes a large suite of parallel linear and nonlinear equation solvers that are easily used in application codes written in C, C++, Fortran and now Python. PETSc provides many of the mechanisms needed within parallel application code, such as simple parallel matrix and vector assembly routines that allow the overlap of communication and computation. In addition, PETSc includes support for parallel distributed arrays useful for finite difference methods.

Components

PETSc consists of a variety of libraries (similar to classes in C++). Each library manipulates a particular family of objects (for instance, vectors) and the operations one would like to perform on the objects.

Features

PETSc provides the following features:

* Parallel vectors and matrices
* Scatters (handles communicating ghost point information)
* Gathers (the opposite of scatters)
* Several sparse storage formats
* Scalable parallel preconditioners
* Krylov subspace methods
* Parallel Newton-based nonlinear solvers
* Parallel time-stepping (ODE) solvers
* Automatic profiling of floating point and memory usage
* Consistent interface
* Intensive error checking
* Portable to UNIX and Windows

Modules

Some of the PETSc modules deal with

* Index sets, including permutations, for indexing into vectors, renumbering, etc;
* Vectors;
* Matrices (generally sparse);
* Distributed arrays (useful for parallelizing regular grid-based problems);
* Krylov subspace methods;
* Preconditioners, including multigrid and sparse direct solvers;
* Nonlinear solvers; and
* Time-steppers for solving time-dependent (nonlinear) PDEs.

Bibliography

* [ftp://info.mcs.anl.gov/pub/tech_reports/reports/ANL9511.ps.Z "PETSc Users Manual", Satish Balay, Kris Buschelman, Victor Eijkhout, William D. Gropp, Dinesh Kaushik, Matthew G. Knepley, Lois Curfman McInnes, Barry F. Smith, and Hong Zhang, ANL-95/11 Revision 2.3.2, Argonne National Laboratory, September, 2006.]
* "Efficient Management of Parallelism in Object Oriented Numerical Software Libraries", Satish Balay, William D. Gropp, Lois Curfman McInnes, Barry F. Smith, Modern Software Tools in Scientific Computing, ed. Bruaset et.al., pp. 163-202, 1997.
* [http://dx.doi.org/10.1016/j.pepi.2007.04.016 "Numerical simulation of geodynamic processes with the Portable Extensible Toolkit for Scientific Computation", R.F. Katz, M.G. Knepley, B. Smith, M. Spiegelman, and E.T. Coon, Physics of The Earth and Planetary Interiors, 163, pp. 52-68, 2007.]

See also

* list of numerical analysis software

External links

* [http://www.mcs.anl.gov/petsc The Official PETSc web site]
* [http://www.netlib.org/utk/people/JackDongarra/la-sw.html Jack Dongarra's Linear Algebra Software Page]