Source author record

Bernhard A. Schmitt

Bernhard A. Schmitt appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math.NA math.DS math.OC Numerical Analysis physics.ao-ph

Catalog footprint

What is connected

4works

5topics

3close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

Implicit A-Stable Peer Triplets for ODE Constrained Optimal Control Problems

This paper is concerned with the construction and convergence analysis of novel implicit Peer triplets of two-step nature with four stages for nonlinear ODE constrained optimal control problems. We combine the property of superconvergence of some standard Peer method for inner grid points with carefully designed starting and end methods to achieve order four for the state variables and order three for the adjoint variables in a first-discretize-then-optimize approach together with A-stability. The notion triplets emphasizes that these three different Peer methods have to satisfy additional matching conditions. Four such Peer triplets of practical interest are constructed. Also as a benchmark method, the well-known backward differentiation formula BDF4, which is only $A(73.35^o)$-stable, is extended to a special Peer triplet to supply an adjoint consistent method of higher order and BDF type with equidistant nodes. Within the class of Peer triplets, we found a diagonally implicit $A(84^o)$-stable method with nodes symmetric in $[0,1]$ to a common center that performs equally well. Numerical tests with three well established optimal control problems confirm the theoretical findings also concerning A-stability.

preprint2022arXiv

Krylov Methods for Adjoint-Free Singular Vector Based Perturbations in Dynamical Systems

The estimation of weather forecast uncertainty with ensemble systems requires a careful selection of perturbations to establish a reliable sampling of the error growth potential in the phase space of the model. Usually, the singular vectors of the tangent linear model propagator are used to identify the fastest growing modes (classical singular vector perturbation (SV) method). In this paper we present an efficient matrix-free block Krylov method for generating fast growing perturbations in high dimensional dynamical systems. A specific matrix containing the non-linear evolution of perturbations is introduced, which we call Evolved Increment Matrix (EIM). Instead of solving an equivalent eigenvalue problem, we use the Arnoldi method for a direct approximation of the leading singular vectors of this matrix, which however is never computed explicitly. This avoids linear and adjoint models but requires forecasts with the full non-linear system. The performance of the approximated perturbations is compared with singular vectors of a full EIM (not with the classical SV method). We show promising results for the Lorenz96 differential equations and a shallow water model, where we obtain good approximations of the fastest growing perturbations by using only a small number of Arnoldi iterations.

preprint2020arXiv

Discrete Adjoint Implicit Peer Methods in Optimal Control

It is well known that in the first-discretize-then-optimize approach in the control of ordinary differential equations the adjoint method may converge under additional order conditions only. For Peer two-step methods we derive such adjoint order conditions and pay special attention to the boundary steps. For $s$-stage methods, we prove convergence of order $s$ for the state variables if the adjoint method satisfies the conditions for order $s\!-\!1$, at least. We remove some bottlenecks at the boundaries encountered in an earlier paper of the first author et al. [J. Comput. Appl. Math., 262:73-86, 2014] and discuss the construction of 3-stage methods for the order pair (3,2) in detail including some matrix background for the combined forward and adjoint order conditions. The impact of nodes having equal differences is highlighted. It turns out that the most attractive methods are related to BDF. Three 3-stage methods are constructed which show the expected orders in numerical tests.

preprint2015arXiv

Algebraic criteria for A-stability of peer two-step methods

A new criterion for A-stability of peer two-step methods is presented which is verifiable exactly in exact arithmetic by checking semi-definiteness of a certain test matrix. It depends on the existence of two positive definite weight matrices for a given method. Although the initial approach is different using properties of the numerical radius the criterion itself resembles the one from algebraic stability of General Linear Methods. Known numerical algorithms for the computation of the unknown weight matrices suffer from rank deficiencies of the test matrix. For $s$-stage peer methods of order $s-1$ this rank defect is identified with an explicit block diagonal decomposition of the test matrix in trivial and definite blocks. In the design of methods its coefficients are unknown and an explicit parametrization of A-stable peer methods of order $s-1$ is presented with a weight matrix as parameter. This leads to a general existence result for any number of stages. The restrictions for efficient L-stable peer methods like diagonally-implicit and parallel ones are also discussed and such methods with $3$ and $4$ stages are constructed.

Institution

Affiliation not imported yet

This author record came from a source that does not expose affiliation metadata. Once the author claims the profile or we enrich the record from another provider, this section will link to the concrete institution.

Topic footprint