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/

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

LSD: A Fast Line Segment Detector with a False Detection Control

IEEE Transactions on Pattern Analysis and Machine Intelligence (2010)

Grompone von Gioi Rafael, Jakubowicz Jérémie, Morel Jean-Michel, and Randall Gregory

0 comment

Add comment

You need to log in to post a comment.