Optimal stability polynomials for numerical integration of initial value problems

Abstract

We consider the problem of finding optimally stable polynomial approximations to the exponential for application to one-step integration of initial value ordinary and partial differential equations. The objective is to find the largest stable step size and corresponding method for a given problem when the spectrum of the initial value problem is known. The problem is expressed in terms of a general least deviation feasibility problem. Its solution is obtained by a new fast, accurate, and robust algorithm based on convex optimization techniques. Global convergence of the algorithm is proven in the case that the order of approximation is one and in the case that the spectrum encloses a starlike region. Examples demonstrate the effectiveness of the proposed algorithm even when these conditions are not satisfied.

References in corpus (1)

• Optimal stability polynomials for numerical integration of initial value problems

Cited by in corpus (15)

• A purely hyperbolic discontinuous Galerkin approach for self-gravitating gas dynamics
• Optimal stability polynomials for numerical integration of initial value problems
• Spatially partitioned embedded Runge-Kutta methods
• Stability-Optimized High Order Methods and Stiffness Detection for Pathwise Stiff Stochastic Differential Equations
• Many-Stage Optimal Stabilized Runge-Kutta Methods for Hyperbolic Partial Differential Equations
• Exponential Time Differencing for the Tracer Equations Appearing in Primitive Equation Ocean Models
• Multirate Time-Integration based on Dynamic ODE Partitioning through Adaptively Refined Meshes for Compressible Fluid Dynamics
• Learning Runge-Kutta Integration Schemes for ODE Simulation and Identification
• Optimized Runge-Kutta Methods with Automatic Step Size Control for Compressible Computational Fluid Dynamics
• A parallel-in-time approach for wave-type PDEs
• Fourth-Order Paired-Explicit Runge-Kutta Methods
• On Fourier analysis of polynomial multigrid for arbitrary multi-stage cycles
• Walking into the complex plane to "order" better time integrators
• An a posteriori strategy for adaptive schemes in time and space
• A fully adaptive explicit stabilized integrator for advection-diffusion-reaction problems