PETSc/TS: A Modern Scalable ODE/DAE Solver Library

Abstract

High-quality ordinary differential equation (ODE) solver libraries have a long history, going back to the 1970s. Over the past several years we have implemented, on top of the PETSc linear and nonlinear solver package, a new general-purpose, extensive, extensible library for solving ODEs and differential algebraic equations (DAEs). Package includes support for both forward and adjoint sensitivities that can be easily utilized by the TAO optimization package, which is also part of PETSc. The ODE/DAE integrator library strives to be highly scalable but also to deliver high efficiency for modest-sized problems. The library includes explicit solvers, implicit solvers, and a collection of implicit-explicit solvers, all with a common user interface and runtime selection of solver types, adaptive error control, and monitoring of solution progress. The library also offers enormous flexibility in selection of nonlinear and linear solvers, including the entire suite of PETSc iterative solvers, as well as several parallel direct solvers.

Cited by in corpus (28)

• MFEM: a modular finite element methods library
• A phase field model for cohesive fracture in micropolar continua
• Multi-phase-field microporomechanics model for simulating ice lens growth and thaw in frozen soil
• Fully-Discrete Explicit Locally Entropy-Stable Schemes for the Compressible Euler and Navier-Stokes Equations
• Towards an open-source landscape for 3D CSEM modelling
• Conservative and entropy stable solid wall boundary conditions for the compressible Navier-Stokes equations: Adiabatic wall and heat entropy transfer
• On the robustness and performance of entropy stable discontinuous collocation methods for the compressible Navier-Stokes equations
• A Graph-Based Modeling Abstraction for Optimization: Concepts and Implementation in Plasmo.jl
• Elastic wave propagation in anisotropic solids using energy-stable finite differences with weakly enforced boundary and interface conditions
• Phase-field modeling of paramagnetic austenite-ferromagnetic martensite transformation coupled with mechanics and micromagnetics
• The Linear Algebra Mapping Problem. Current state of linear algebra languages and libraries
• Entropy Stable p-Nonconforming Discretizations with the Summation-by-Parts Property for the Compressible Euler equations
• Randomized maximum likelihood based posterior sampling
• A 55-line code for large-scale parallel topology optimization in 2D and 3D
• GWsim: a code to simulate gravitational waves propagating in a potential well
• A Stable FE Method For the Space-Time Solution of the Cahn-Hilliard Equation
• A Keller-Segel model for C elegans L1 aggregation
• An Unconditionally Stable Space-Time FE Method for the Korteweg-de Vries Equation
• An Adaptive Stable Space-Time FE Method for the Shallow Water Equations
• Investigating the potential for a limited quantum speedup on protein lattice problems
• Multi-Objective Control Co-design Using Graph-Based Optimization for Offshore Wind Farm Grid Integration
• A Discrete Probabilistic Approach to Dense Flow Visualization
• Computational investigations of a multi-class traffic flow model: mean-field and microscopic dynamics
• Computing Derivatives for PETSc Adjoint Solvers using Algorithmic Differentiation
• Solving stochastic inverse problems for property-structure linkages using data-consistent inversion and machine learning
• A p-adaptive, implicit-explicit mixed finite element method for reaction-diffusion problems
• Base flow decomposition for complex moving objects in linear hydrodynamics: application to helix-shaped flagellated microswimmers
• Graph-Based Modeling and Simulation of Complex Systems