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.
Contents |
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.
