Optimal Stability Polynomials for Numerical Integration of Initial Value Problems

By Ketcheson David I. and Ahmadia Aron J.
Communications in Applied Mathematics and Computational Science (2012)

  • David Ketcheson

    King Abdullah University of Science and Technology

    Saudi Arabia

  • Aron Ahmadia

    U.S. Army Corps of Engineers



November 6, 2013

Last update

November 6, 2013










This code reproduces the figures 3a, 3b, 4a, 4b, 6a, 6b and 7b of the article “Optimal stability polynomials for numerical integration of initial value problems” (David I. Ketcheson and Aron J. Ahmadia, 2012). The user can fix the number of stages (s) and the order p, and then retrieve (1) the scaled size of real axis interval inclusion for optimized methods Hopt/s^2 (Table 1, page 12), (2) the scaled size of imaginary axis inclusion for optimized methods Hopt/s (Table 2, page 13), (3) the relative size of largest disk that can be included in the stability region scaled by the number of stages (Figure 5, page 15) and (4) the optimal effective step size (Figure 7a, page 16). Please note that updated code is available at http://numerics.kaust.edu.sa/RK-opt/

On Reduced Input-Output Dynamic Mode Decomposition

Working Paper (2017)

Himpe Christian, Benner Peter, and Mitchell Tim

Fast Low-Rank Empirical Cross Gramians

Working Paper (2017)

Himpe Christian, Leibner Tobias, Rave Stephan, and Saak Jens

emgr - The Empirical Gramian Framework

Working Paper (2016)

Himpe Christian

Cross-Gramian-Based Model Reduction: A Comparison

Working Paper (2016)

Himpe Christian and Ohlberger Mario

Parareal in Time Intermediate Targets Methods for Optimal Control Problems

International Series of Numerical Mathematics (2013)

Mohamed Kamel Riahi, Maday Yvon, and Julien salomon

The Empirical Cross Gramian for Parametrized Nonlinear Systems

Working Paper (2015)

Himpe Christian and Ohlberger Mario

A Note on the Cross Gramian for Non-Symmetric Systems

Working Paper (2015)

Himpe Christian and Ohlberger Mario

Data-Driven Combined State and Parameter Reduction for Inverse Problems

Working Paper (2014)

Himpe Christian and Ohlberger Mario

Model Reduction for Complex Hyperbolic Networks

arXiv (2013)

Himpe Christian and Ohlberger Mario

0 comment

Add comment

You need to log in to post a comment.