2011/01/03 by Jan L. Cieśliński, Cieśliński, Jan L. · 1 citation
Computer Science · Engineering · Mathematics · #34K28 #65L12 #65P10 #Advanced Numerical Methods in Computational Mathematics #Computational Fluid Dynamics and Aerodynamics #FOS: Mathematics #Numerical Analysis (math.NA) #Numerical methods for differential equations #cs.NA #math.NA #msc:34K28 #msc:65L12 #msc:65P10
paper · pdf · doi:10.48550/arxiv.1101.0578
33 pages, 0 figures
arxiv created 2011/01/03 · openalex publication_date 2011/01/03 · arxiv updated 2011/04/08 · openalex created_date 2025/10/10 · openalex updated_date 2026/07/28
We present a new class of exponential integrators for ordinary differential equations. They are locally exact, i.e., they preserve the linearization of the original system at every point. Their construction consists in modifying existing numerical schemes in order to make them locally exact. The resulting schemes preserve all fixed points and are A-stable. The most promising results concern the discrete gradient method (modified implicit midpoint rule) where we succeeded to preserve essential geometric properties and the final results have a relatively simple form. In the case of one-dimensional Hamiltonian systems numerical experiments show that our modifications can increase the accuracy by several orders of magnitude. The main result of this paper is the construction of energy-preserving locally exact discrete gradient schemes for arbitrary multidimensional Hamiltonian systems in canonical coordinates.