Refine
Year of publication
- 2007 (7) (remove)
Language
- English (7)
Keywords
- Dune (2)
- grid-based methods (2)
- hierarchical grids (2)
- Bohl exponent (1)
- Lyapunov exponent (1)
- Nonlinear eigenvalue problems (1)
- Perron-Frobenius theorem (1)
- Positive descriptor systems (1)
- Sacker-Sell spectrum (1)
- Steklov function (1)
We present a new extension of the well-known
Perron-Frobenius theorem to regular matrix pairs $(E,A)$.
The new extension is based on projector chains and is motivated from
the solution of positive differential-algebraic systems or descriptor
systems. We present several examples where the new condition holds, whereas conditions
in previous literature are not satisfied.
Lyapunov and exponential dichotomy spectral theory is extended
from ordinary differential equations (ODEs) to nonautonomous
differential-algebraic equations (DAEs). By using orthogonal
changes of variables, the original DAE system is transformed into
appropriate condensed forms, for which concepts such as Lyapunov
exponents, Bohl exponents, exponential dichotomy and spectral
intervals of various kinds can be analyzed via the resulting
underlying ODE. Some essential differences between the spectral
theory for ODEs and that for DAEs are pointed out. Numerical
methods for computing the spectral intervals associated with
Lyapunov and Sacker-Sell (exponential dichotomy) spectra are
derived by modifying and extending those methods proposed for ODEs. Perturbation theory and error analysis are discussed, as
well. Finally, some numerical examples are presented to illustrate
the theoretical results and the properties of the numerical
methods.
A Generic Grid Interface for Parallel and Adaptive Scientific Computing. Part I: Abstract Framework
(2007)
We give a mathematically rigorous definition of a grid for algorithms solving
partial differential equations. Unlike previous approaches, our grids have a
hierarchical structure. This makes them suitable for geometric multigrid
algorithms and hierarchical local grid refinement. The description is also
general enough to include geometrically nonconforming grids. The definitions
in this article serve as the basis for an implementation of an abstract grid
interface as C++ classes in the DUNE.
In a companion paper [Matheon-Preprint 403] we introduced an abstract definition of a parallel and adaptive hierarchical grid for scientific computing. Based on this
definition we derive an efficient interface specification as a set of C++ classes.
This interface separates the applications from the grid data structures.
Thus, user implementations become independent of the underlying grid
implementation. Modern C++ template techniques are used to provide an
interface implementation without big performance losses.
The implementation is realized as part of the
software environment DUNE.
Numerical tests demonstrate the flexibility and the efficiency of our approach.
In this paper, we discuss stability properties of positive descriptor systems in the continuous-time as well as in the discrete-time case. We present different characterisations of positivity and establish generalised stability criteria for the case of positive descriptor systems. We show that if the spectral projector onto the right finite deflating subspace of the matrix pair $(E,A)$ is non-negative, then all stability criteria for standard positive systems take a comparably simple form in the positive descriptor case. Furthermore, we provide sufficient conditions that guarantee entry-wise non-negativity along with positive semi-definiteness of solutions of generalised projected Lyapunov equations. As an application of the framework established throughout this paper, we exemplarily generalise two criteria for the stability of two switched standard positive systems under arbitrary switching to the descriptor case.
Adjoint Broyden a la GMRES
(2007)
It is shown here that a compact storage implementation of a quasi-Newton
method based on the adjoint Broyden update reduces in the affine
case exactly to the well established GMRES procedure. Generally,
storage and linear algebra effort per step are small multiples of $n\cdot k$,
where $n$ is the number of variables and $k$ the number of steps taken
in the current cycle. In the affine case the storage is exactly $(n+k)\cdot k$
and in the nonlinear case the same bound can be achieved if adjoints,
i.e. transposed Jacobian-vector products are available. A
transposed-free variant that relies exclusively on Jacobian-vector
products (or possibly their approximation by divided differences)
requires roughly twice the storage and turns out to be somewhat slower
in our numerical experiments reported at the end.
Structured eigenvalue conditioning and backward error of a class of polynomial eigenvalue problems
(2007)
Characterisations of simple eigenvalues of complex matrix polynomials with *-even/odd and *-palindromic/antipalindromic structures that have the same normwise condition number with respect to structure preserving and arbitrary perturbations are obtained. Here * denotes either the transpose T or the conjugate tranpose *. In the process we obtain formulae for the normwise structured condition number of simple eigenvalues of T-palindromic/antipalindromic and *-even/odd polynomials. Moreover, conditions under which the normwise structured backward error of approximate eigenvalues of such polynomials is equal to the unstructured error are also derived. These lead to complete characterisations of approximate eigenvalues that have the same structured and unstructured backward errors for the *-even/odd and T-even/odd polynomials.