Source author record

Heike Faßbender

Heike Faßbender appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math.NA Numerical Analysis

Catalog footprint

What is connected

4works

2topics

4close 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

The Hamiltonian Extended Krylov Subspace Method

An algorithm for constructing a $J$-orthogonal basis of the extended Krylov subspace $\mathcal{K}_{r,s}=\operatorname{range}\{u,Hu, H^2u,$ $ \ldots, $ $H^{2r-1}u, H^{-1}u, H^{-2}u, \ldots, H^{-2s}u\},$ where $H \in \mathbb{R}^{2n \times 2n}$ is a large (and sparse) Hamiltonian matrix is derived (for $r = s+1$ or $r=s$). Surprisingly, this allows for short recurrences involving at most five previously generated basis vectors. Projecting $H$ onto the subspace $\mathcal{K}_{r,s}$ yields a small Hamiltonian matrix. The resulting HEKS algorithm may be used in order to approximate $f(H)u$ where $f$ is a function which maps the Hamiltonian matrix $H$ to, e.g., a (skew-)Hamiltonian or symplectic matrix. Numerical experiments illustrate that approximating $f(H)u$ with the HEKS algorithm is competitive for some functions compared to the use of other (structure-preserving) Krylov subspace methods.

preprint2020arXiv

A link between gramian based model order reduction and moment matching

We analyze a family of Runge-Kutta based quadrature algorithms for the approximation of the gramians of linear time invariant dynamical systems. The approximated gramians are used to obtain an approximate balancing transformation similar to the approach used in balanced POD. It is shown that hereby rational interpolation is performed, as the approximants span certain Krylov subspaces. The expansion points are mainly determined by the time step sizes and the eigenvalues of the matrices given by the Butcher tableaus.

preprint2020arXiv

Riccati ADI: Existence, uniqueness and new iterative methods

The approximate solution of large-scale algebraic Riccati equations is considered. We are interested in approximate solutions which yield a Riccati residual matrix of a particular small rank. It is assumed that such approximate solutions can be written in factored form $ZYZ^*$ with a rectangular matrix $Z$ and a small quadratic matrix $Y$. We propose to choose $Z$ such that its columns span a certain rational Krylov subspace. Conditions under which such an approximate solution exists and is unique are determined. It is shown that the proposed method can be interpreted as an oblique projection method. Two new iterative procedures with efficient updates of the solution and the residual factor are derived. With our approach complex system matrices can be handled, realification is provided and parallelization is introduced.

preprint2016arXiv

Symmetric and skew-symmetric block-Kronecker linearizations

Many applications give rise to structured matrix polynomials. The problem of constructing structure-preserving strong linearizations of structured matrix polynomials is revisited in this work and in the forthcoming ones \cite{PartII,PartIII}. With the purpose of providing a much simpler framework for structure-preserving linearizations for symmetric and skew-symmetric matrix polynomial than the one based on Fiedler pencils with repetition, we introduce in this work the families of (modified) symmetric and skew-symmetric block Kronecker pencils. These families provide a large arena of structure-preserving strong linearizations of symmetric and skew-symmetric matrix polynomials. When the matrix polynomial has degree odd, these linearizations are strong regardless of whether the matrix polynomial is regular or singular, and many of them give rise to structure-preserving companion forms. When some generic nonsingularity conditions are satisfied, they are also strong linearizations for even-degree regular matrix polynomials. Many examples of structure-preserving linearizations obtained from Fiedler pencils with repetitions found in the literature are shown to belong (modulo permutations) to these families of linearizations. In particular, this is shown to be true for the well-known block-tridiagonal symmetric and skew-symmetric companion forms. Since the families of symmetric and skew-symmetric block Kronecker pencils belong to the recently introduced set of minimal bases pencils \cite{Fiedler-like}, they inherit all its desirable properties for numerical applications. In particular, it is shown that eigenvectors, minimal indices, and minimal bases of matrix polynomials are easily recovered from those of any of the linearizations constructed in this work.