Refine
Year of publication
- 2010 (70) (remove)
Language
- English (70)
Keywords
- dissipativity (3)
- optimal control (3)
- backward error (2)
- behavior approach (2)
- even pencil (2)
- even polynomial (2)
- finite elements (2)
- linear quadratic optimal control (2)
- problem (2)
- structured eigenvalue problem (2)
This paper presents concepts and implementation of the finite element toolbox Kaskade 7, a flexible C++ code for solving elliptic and parabolic PDE systems. Issues such as problem formulation, assembly and adaptivity are discussed at the example of optimal control problems. Trajectory compression for parabolic optimization problems is considered as a case study.
We propose a generalization of the Structured Doubling Algorithm (SDA) to compute invariant subspaces
of structured matrix pencils
that arise in the context of solving linear quadratic optimal control problems.
The new algorithm is
designed to attain better accuracy when the classical Riccati equation approach for the solution of the optimal control problem is not well suited because
the stable and unstable invariant subspaces are not well separated (due to eigenvalues near or on the imaginary
axis) or in the case when the Riccati solution does not exist at all. We analyze the convergence
of the method and compare the new method with the classical SDA algorithm as well as some recent structured QR-methods.
Recently, the format of TT tensors
\cite{hackbuschHT,osele1,tyrtosele2,tyrtosele3} has turned out to be
a promising new format for the approximation of solutions of high
dimensional problems. In this paper, we prove some new results for
the TT representation of a tensor $U \in \R^{n_1\times \ldots\times
n_d}$ and for the manifold of tensors of TT-rank $\underline{r}$.\As a first result, we prove that the TT (or compression) ranks $r_i$
of a tensor $U$ are unique and equal to the respective separation
ranks of $U$ if the components of the TT decomposition are required to
fulfil a certain maximal rank condition. We then show that the set
$\mathcal{T}$ of TT tensors of fixed rank $\underline{r}$ forms an embedded
manifold in $\R^{n^d}$, therefore preserving the essential theoretical
properties of the Tucker format, but often showing an improved scaling
behaviour. Extending a similar approach for matrices \cite{conte_lub},
we introduce certain gauge conditions to obtain a unique
representation of the tangent space $\cT_U\mathcal{T}$ of $\mathcal{T}$
and deduce a
local parametrization of the TT manifold. The parametrisation of
$\cT_{U}\mathcal{T}$ is often crucial for an algorithmic treatment of
high-dimensional time-dependent PDEs and minimisation problems
\cite{lubuch_blau}. We conclude with remarks on those applications and
present some numerical examples.
We consider the numerical solution of projected Lyapunov equations using Krylov subspace iterative methods. Such equations
play a fundamental role in balanced truncation model reduction of descriptor systems. We present generalizations of
the extended block and global Arnoldi methods to projected Lyapunov equations and compare these methods with the alternating direction implicit method with respect to performance on different examples.
A deflation strategy is also proposed to overcome possible breakdown in the
recurrence.
The PSurface Library
(2010)
We describe psurface, a C++ library that allows to store and access piecewise linear mappings between simplicial surfaces in $\R^2$ and $\R^3$. These mappings are stored in a graph data structure and can be constructed explicitly, by projection, or by surface simplification. Piecewise linear maps can be used, e.g., to construct boundary
approximations for finite element grids, and grid intersections for domain decomposition methods. In computer graphics the mappings allow to build level-of-detail representations as well as texture- and bump maps. We document the data structures and algorithms used and show how \psurface is used in the numerical analysis framework Dune
and the visualization software Amira.
We describe the appearance and stability of spatio-temporal periodic
patterns (rotating waves) in unidirectional rings of coupled oscillators
with delayed couplings. We show how delays in the coupling lead
to a splitting of each rotating wave into several new ones. The appearance
of rotating waves is mediated by Hopf bifurcations of the symmetric
equilibrium.
We also conclude that the coupling delays can be effectively
replaced by increasing the number of oscillators in the chain.
The phenomena are shown for Stuart-Landau
oscillators as well as for coupled FitzHugh-Nagumo systems interacting
via excitatory chemical synapses.
Classical surface parameterization algorithms often place singularities
in order to enhance the quality of the resulting parameter map. Unfortunately, singularities of positive integral index (as the north pole of a sphere) were not handled since they cannot be described with piecewise linear parameter functions on a triangle mesh. Preprocessing is needed to adapt the mesh connectivity. We present an extension to the QuadCover parameterization algorithm [KNP07], which allows to handle those singularities. A singularity of positive integral index can be resolved using bilinear parameter functions on quadrilateral elements. This generalization
of piecewise linear functions for quadrilaterals enriches the space of parameterizations. The resulting parameter map can be visualized by textures using a rendering system which supports quadrilateral elements, or it can be used for remeshing into a pure quad mesh.
We study the class of linear differential-algebraic m-input m-output systems
which have a transfer function with proper inverse.
A sufficient condition for the transfer function to have proper inverse
it that the system has 'strict and non-positive relative degree'.
We present two main results:
First, a so called 'zero dynamics form' is derived: this form is - within the class of system equivalence - a simple "almost normal" form of the DAE; it is a counterpart to the well-known Byrnes-Isidori form for
ODE systems with strictly proper transfer function.
The 'zero dynamics form' is exploited to characterize structural properties such as
asymptotically stable zero dynamics,
minimum phase, and high-gain stabilizability.
The zero dynamics are characterized by (A,E,B)-invariant subspaces.
Secondly, it is shown that the 'funnel controller' (that is a static nonlinear output error feedback) achieves, for all DAE systems with asymptotically stable zero dynamics and transfer function with proper inverse, tracking of a reference signal by the output signal within a pre-specified funnel. This funnel determines the transient behaviour.
We study discrete curvatures computed from nets of curvature lines on a given smooth surface and prove their uniform convergence to smooth principal curvatures. We provide explicit error bounds, with constants depending only on properties of the smooth limit surface and the shape regularity of the discrete net.
We present a novel algorithm for automatic parameterization of tube-like surfaces of arbitrary genus such as the surfaces of knots, trees, blood vessels, neurons, or any tubular graph with a globally consistent stripe texture. We use the principal curvature frame field of the underlying tube-like surface to guide the creation of a global, topologically consistent stripe parameterization of the surface. Our algorithm extends the QuadCover algorithm and is based, first, on the use of so-called projective vector fields instead of frame fields, and second, on different types of branch points. That does not only simplify the mathematical theory, but also reduces computation time by the decomposition of the underlying stiffness matrices.
Riemann surfaces naturally appear in the analysis of complex functions that are branched over the complex plane. However, they usually possess a complicated topology and are thus hard to understand. We present an algorithm for constructing Riemann surfaces as meshes in R3 from explicitly given branch points with corresponding branch indices. The constructed surfaces cover the complex plane by the canonical
projection onto R2 and can therefore be considered as multivalued graphs
over the plane – hence they provide a comprehensible visualization of the
topological structure. Complex functions are elegantly visualized using domain coloring on
a subset of C. By applying domain coloring to the automatically constructed Riemann surface models, we generalize this approach to deal with functions which cannot be entirely visualized in the complex plane.
Diffusion weighted imaging is a magnetic resonance based method to investigate
tissue micro-structure especially in the human brain via water diffusion.
Since the standard diffusion tensor model for the acquired data failes in
large portion of the brain voxel more sophisticated models have bee developed.
Here, we report on the package dti and how some of these models
can be used with the package.
The package fmri is provided for analysis of single run functional
Magnetic Resonance Imaging data. It implements structural adaptive smoothing
methods with signal detection for adaptive noise reduction which avoids blurring
of edges of activation areas. fmri provides fmri analysis from time series
modeling to signal detection and publication-ready images.
Modeling the orientation distribution function by mixtures of angular central Gaussian distributions
(2010)
In this paper we develop a tensor mixture model for diffusion weighted imaging
data using an automatic model selection criterion for the order of tensor
components in a voxel. We show that the weighted orientation distribution
function for this model can be expanded into a mixture of angular central
Gaussian distributions. We show properties of this model in extensive
simulations and in a high angular resolution experimental data set. The results
suggest that the model may improve imaging of cerebral fiber tracts. We
demonstrate how inference on canonical model parameters may give rise to new
clinical applications.
In this work, we study the spectra and eigenmodes of the Hessian of various discrete surface energies and discuss applications to shape analysis. In particular, we consider a physical model that describes the vibration modes and frequencies of a surface through the eigenfunctions and eigenvalues of the Hessian of a deformation energy, and we
derive a closed form representation for the Hessian (at the rest state of the energy) for a general class of deformation energies. Furthermore, we design a quadratic energy, such that the eigenmodes of the Hessian of
this energy are sensitive to the extrinsic curvature of the surface. Based on these spectra and eigenmodes, we derive two shape signatures. One that measures the similarity of points on a surface, and another that
can be used to identify features of the surface. In addition, we discuss a
spectral quadrangulation scheme for surfaces.
We discuss the perturbation analysis for
eigenvalues and eigenvectors of structured homogeneous matrix polynomials with
Hermitian, skew-Hermitian, H-even and H-odd structure.
We construct minimal structured perturbations (structured backward errors) such that an
approximate eigenpair is an exact eigenpair of an appropriate perturbed structured matrix
polynomial. We present various comparisons with unstructured backward
errors and previous error bounds derived for the non-homogeneous case
and show that our bounds present a significant improvement.
Many applications give rise to matrix polynomials whose coefficients have
a kind of reversal symmetry, a structure we call palindromic.
Several properties of scalar palindromic polynomials are derived,
and together with properties of compound matrices, used to
establish the Smith form of regular and singular T-palindromic matrix polynomials,
over arbitrary fields.
The invariant polynomials are shown to
inherit palindromicity,
and their structure is described in detail.
Jordan structures of palindromic matrix polynomials are characterized,
and necessary conditions for the
existence of structured linearizations established.
In the odd degree case, a constructive procedure for building
palindromic linearizations shows that the necessary conditions are sufficient as well.
The Smith form for *-palindromic polynomials is also analyzed. Finally, results for palindromic matrix polynomials over fields of
characteristic two are presented.
We define a risk averse nonanticipative feasible policy for multistage stochastic programs and propose a methodology to implement it. The approach is based on dynamic programming equations written for a risk averse formulation of the problem.
This formulation relies on a new class of multiperiod risk functionals called extended polyhedral risk measures. Dual representations of such risk functionals are given and used to derive conditions of coherence. In the one-period case, conditions for convexity and consistency with second order stochastic dominance are also provided. The risk averse dynamic programming equations are specialized considering convex combinations of one-period extended polyhedral risk measures such as spectral risk measures.
To implement the proposed policy, the approximation of the risk averse recourse functions for stochastic linear programs is discussed. In this context, we detail a stochastic dual dynamic programming algorithm which converges to the optimal value of the risk averse problem.
We derive a Crooks-Jarzynski-type identity for computing free energy differences between metastable states that is based on nonequilibrium diffusion processes. Furthermore we outline a brief derivation of an infinite-dimensional stochastic partial differential equation that can be used to efficiently generate the ensemble of trajectories connecting the metastable states.
To address the plurality of interpretations of the subjective notion of risk, we describe it by means of a risk order and concentrate on the context invariant features of diversification and monotonicity. Our main results are uniquely characterized robust representations of lower semicontinuous risk orders on vector spaces and convex sets. This representation covers most instruments related to risk and allow for a differentiated interpretation depending on the underlying context which is illustrated in different settings: For random variables, risk perception can be interpreted as model risk, and we compute among others the robust representation of the economic index of riskiness. For lotteries, risk perception can be viewed as distributional risk and we study the "Value at Risk". For consumption patterns, which excerpt an intertemporality dimension in risk perception, we provide an interpretation in terms of discounting risk and discuss some examples.