{"title":"Variable-stepsize implicit Peer triplets in ODE constrained optimal control","authors":"Jens Lang , Bernhard A. Schmitt","doi":"10.1016/j.cam.2024.116417","DOIUrl":null,"url":null,"abstract":"<div><div>This paper is concerned with the theory, construction and application of implicit Peer two-step methods that are super-convergent for variable stepsizes, i.e., preserve their classical order achieved for uniform stepsizes when applied to ODE constrained optimal control problems in a first-discretize-then-optimize setting. We upgrade our former implicit two-step Peer triplets constructed in Lang and Schmitt (Algorithms 2022) to get ready for dynamical systems with varying time scales without loosing efficiency. Peer triplets consist of a standard Peer method for interior time steps supplemented by matching methods for the starting and end steps. A decisive advantage of Peer methods is their absence of order reduction since they use stages of the same high stage order. The consistency analysis of variable-stepsize implicit Peer methods results in additional order conditions and severe new difficulties for uniform zero-stability, which intensifies the demands on the Peer triplet. Further, we discuss the construction of 4-stage methods with order pairs (4,3) and (3,3) for state and adjoint variables in detail and provide four Peer triplets of practical interest. We rigorously prove convergence of order <span><math><mrow><mi>s</mi><mo>−</mo><mn>1</mn></mrow></math></span> for <span><math><mi>s</mi></math></span>-stage Peer methods applied on grids with bounded or smoothly changing stepsize ratios. Numerical tests show the expected order of convergence for the new variable-stepsize Peer triplets.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"460 ","pages":"Article 116417"},"PeriodicalIF":2.1000,"publicationDate":"2024-12-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Journal of Computational and Applied Mathematics","FirstCategoryId":"100","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S0377042724006654","RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q1","JCRName":"MATHEMATICS, APPLIED","Score":null,"Total":0}
引用次数: 0
Abstract
This paper is concerned with the theory, construction and application of implicit Peer two-step methods that are super-convergent for variable stepsizes, i.e., preserve their classical order achieved for uniform stepsizes when applied to ODE constrained optimal control problems in a first-discretize-then-optimize setting. We upgrade our former implicit two-step Peer triplets constructed in Lang and Schmitt (Algorithms 2022) to get ready for dynamical systems with varying time scales without loosing efficiency. Peer triplets consist of a standard Peer method for interior time steps supplemented by matching methods for the starting and end steps. A decisive advantage of Peer methods is their absence of order reduction since they use stages of the same high stage order. The consistency analysis of variable-stepsize implicit Peer methods results in additional order conditions and severe new difficulties for uniform zero-stability, which intensifies the demands on the Peer triplet. Further, we discuss the construction of 4-stage methods with order pairs (4,3) and (3,3) for state and adjoint variables in detail and provide four Peer triplets of practical interest. We rigorously prove convergence of order for -stage Peer methods applied on grids with bounded or smoothly changing stepsize ratios. Numerical tests show the expected order of convergence for the new variable-stepsize Peer triplets.
期刊介绍:
The Journal of Computational and Applied Mathematics publishes original papers of high scientific value in all areas of computational and applied mathematics. The main interest of the Journal is in papers that describe and analyze new computational techniques for solving scientific or engineering problems. Also the improved analysis, including the effectiveness and applicability, of existing methods and algorithms is of importance. The computational efficiency (e.g. the convergence, stability, accuracy, ...) should be proved and illustrated by nontrivial numerical examples. Papers describing only variants of existing methods, without adding significant new computational properties are not of interest.
The audience consists of: applied mathematicians, numerical analysts, computational scientists and engineers.