Refine
Year of publication
Language
- English (97)
Keywords
- optimal control (11)
- differential-algebraic equation (4)
- finite elements (4)
- metastability (4)
- trajectory storage (4)
- Lyapunov exponent (3)
- Sacker-Sell spectrum (3)
- exponential dichotomy (3)
- geodesic finite elements (3)
- rare events (3)
- spectral interval (3)
- state constraints (3)
- Bohl exponent (2)
- Cosserat rods (2)
- Currents (2)
- Dune (2)
- Hamiltonian matrix (2)
- Krylov subspace methods (2)
- MINRES (2)
- Matching (2)
- Newton-CG (2)
- Schrödinger equation (2)
- Steklov function (2)
- cluster analysis (2)
- compression (2)
- deflation (2)
- domain decomposition (2)
- dynamical systems (2)
- flux coupling analysis (2)
- generic programming (2)
- grid-based methods (2)
- hierarchical grids (2)
- index (2)
- metabolic network (2)
- molecular dynamics (2)
- nonlinear elasticity (2)
- path-following (2)
- perturbation theory (2)
- regularity (2)
- strangeness (2)
- strangeness index (2)
- $CUR$ decomposition (1)
- $H_\infty$ control (1)
- Acetabulum (1)
- Adaptive Finite Elements (1)
- Adaptive finite elements (1)
- Atomistic-Continuum Coupling (1)
- BVD pencil (1)
- Boolean network (1)
- Bridging Domain (1)
- C-stationarity (1)
- CG (1)
- Cahn-Larche system (1)
- Cayley transformation (1)
- Conformation Dynamics; Markov State Models; Molecular Kinetics; Transition Rates (1)
- Constraint-based modeling, metabolic network, thermodynamic constraints (1)
- Cosserat rod (1)
- Deformable Surfaces (1)
- Diffeomorphic Registration (1)
- Diffeomorphic registration (1)
- Dirichlet-Neumann (1)
- Discrete model (1)
- Discrete model; Gene regulatory network; Piecewise affine model (1)
- Dynamical contact problems (1)
- Elementary flux mode (1)
- Elliptic div-grad operators (1)
- Flux balance analysis (1)
- GMRES (1)
- Gene regulatory network (1)
- Ginzburg-Landau equations (1)
- Hamiltonian eigenvalue problem (1)
- Hamiltonian pseudospectra (1)
- Hamiltonian real Schur form (1)
- Klein-Gordon equation (1)
- Large Deformation (1)
- Large deformation (1)
- Lithography (1)
- Mandible (1)
- Markov Random Field (1)
- Maxwell's equations (1)
- Meshfree (1)
- Metabolic network (1)
- Monte Carlo simulation (1)
- Multiscale Simulation (1)
- Newmark method (1)
- Newton methods in function space (1)
- Nonlinear eigenvalue problems (1)
- PDE constrained optimization (1)
- Perron-Frobenius theorem (1)
- Piecewise affine model (1)
- Polynomial eigenvalue problem (1)
- Positive descriptor systems (1)
- QR algorithm (1)
- Richards equation (1)
- Riemannian center of mass (1)
- Riemannian manifolds (1)
- Riemannian trust region method (1)
- Riemannian trust-region methods (1)
- Ritz vector recycling (1)
- Rothe method (1)
- SQP (1)
- Schroedinger equation (1)
- Schur-Newton multigrid (1)
- Segmentation (1)
- Signorini condition (1)
- TR1 Update (1)
- Tucker decomposition (1)
- Variational inequalities (1)
- abstract interface (1)
- adaptive finite element method (1)
- adaptive finite elements (1)
- adaptive mesh refinement (1)
- adaptive methods (1)
- adaptivity (1)
- adjoint based updating (1)
- adjoint gradient computation (1)
- anisotropic ellipticity in three dimensions (1)
- approximation (1)
- augmentation (1)
- backward error (1)
- biological network (1)
- biomechanics (1)
- bounded-realness (1)
- c++ (1)
- chattering (1)
- chemical master equation (1)
- complex skew-symmetric matrix polynomial (1)
- complex symmetric matrix polynomial (1)
- complexity reduction (1)
- compressed sensing (1)
- conformation dynamics (1)
- conformational space (1)
- conical diffraction (1)
- consistency (1)
- constraint-based analysis (1)
- constraint-based analysis; metabolic networks; elementary flux modes (1)
- constraint-based modeling (1)
- contact problem (1)
- continuity of Nemyckii Operators (1)
- continuous QR algorithm (1)
- control constraints (1)
- convergence order (1)
- coupling (1)
- coupling conditions (1)
- defibrillation (1)
- descriptor system (1)
- dictionary (1)
- differential-algebraic equations (1)
- differentiation index (1)
- directors (1)
- discrete QR algorithm (1)
- discrete dynamical system (1)
- disk function method (1)
- distance to bounded-realness (1)
- drift-diffusion equation (1)
- dune (1)
- eigenfunctions (1)
- eigenvalue cluster (1)
- eigenvalue problem (1)
- electro-magnetic scattering (1)
- elementary flux modes (1)
- energy minimization (1)
- equivalence (1)
- even matrix polynomial (1)
- even matrix polynomials (1)
- even pencil (1)
- exact residual (1)
- existence and uniqueness of solutions (1)
- explicit formula (1)
- finite dynamical systems (1)
- finite state projections (1)
- flux balance analysis (1)
- function space (1)
- generalised projected Lyapunov equations (1)
- goal oriented error estimation (1)
- hierarchical basis (1)
- homotopy (1)
- hybrid system (1)
- hypercubes (1)
- implant design (1)
- index reduction (1)
- inequality constraints (1)
- interior point methods (1)
- interior point methods in function space (1)
- interpolation error (1)
- kinematic (1)
- kinematic equivalence (1)
- knockout simulation (1)
- lattice theory (1)
- leading direction (1)
- least squares (1)
- linear programming (1)
- linear-quadratic control (1)
- liquid crystals (1)
- local convergence (1)
- logical modeling (1)
- lossy compression (1)
- manifold (1)
- matrix pencils (1)
- mesh compression (1)
- meta grid (1)
- metabolic networks (1)
- metabolic networks, (1)
- metastable conformation (1)
- micro-fe (1)
- mixed Dirichlet-Neumann boundary conditions (1)
- mixed derivatives (1)
- modules (1)
- molecular dynamics simulations (1)
- monodomain model (1)
- monotone multigrid (1)
- monotone multigrid method (1)
- morphing (1)
- mutual incoherence (1)
- nonlinear Schrödinger equations (1)
- nonlinear configuration space (1)
- nonlinear equations (1)
- odd matrix polynomial (1)
- odd matrix polynomials (1)
- optimal Sobolev regularity (1)
- palindromic matrix polynomials (1)
- parametrizations (1)
- partial differential equation (1)
- partial differential equations (1)
- passive system (1)
- pencil arithmetic (1)
- peptides (1)
- perfectly matched layer method (1)
- phase separation (1)
- piecewise linear 3D flattening (1)
- piecewise linear maps (1)
- pole condition (1)
- polyconvex elasticity (1)
- prime implicand (1)
- projector chain (1)
- purely imaginary eigenvalues (1)
- quai-Newton (1)
- quasi-Newton methods (1)
- rank $k$ (1)
- rate matrix (1)
- regularization (1)
- restricted isometry property (1)
- rough potential energy function (1)
- semilinear parabolic PDEs (1)
- sign characteristic (1)
- simplicial surfaces (1)
- singular value decomposition (1)
- sliding mode simulation (1)
- smooth $QR$ factorization (1)
- smooth QR factorization (1)
- smooth singular value decomposition (1)
- sparse solution (1)
- spectral approximation (1)
- spectral clustering (1)
- spectral deferred correction (1)
- stability (1)
- stationary distribution (1)
- structured backward error (1)
- structured condition number (1)
- structured doubling algorithm (1)
- structured mapping problem (1)
- structured staircase form, structured Kronecker canonical form, symmetric pencil, Hermitian pencil, alternating pencil, palindromic pencil, linear quadratic control, $H_\infty$ control (1)
- subgrid (1)
- subspace recycling (1)
- switched positive systems (1)
- switched system (1)
- symplectic matrix (1)
- symplectic pencil (1)
- systems biology (1)
- thermodynamic constraints (1)
- time integration (1)
- tractability index (1)
- transition state theory (1)
- transmission at material interfaces (1)
- transparent boundary conditions (1)
- trust region (1)
- unit sphere (1)
- variational analysis (1)
- viscoelasticity (1)
- wave equation (1)
- zero curvature (1)
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.
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.
Deflated and augmented Krylov subspace methods: Basic Facts and a Breakdown-free deflated MINRES
(2011)
In this paper we consider deflation and augmentation techniques for accelerating
the convergence of Krylov subspace methods for the solution of nonsingular linear
algebraic systems. The two techniques are conceptually different from
preconditioning. Deflation "removes" certain parts from the operator, while
augmentation adds a subspace to the Krylov subspace. Both approaches have been
used in a variety of methods and settings. For Krylov subspace methods that
satisfy a (Petrov-) Galerkin condition we show that augmentation can in general
be achieved implicitly by projecting the residuals appropriately and correcting
the approximate solutions in a final step. In this context, we analyze known
methods to deflate CG, GMRes and MinRes. Our analysis reveals that the recently
proposed RMinRes method can break down. We show how such breakdowns can be
avoided by choosing a special initial guess, and we derive a breakdown-free
deflated MinRes method. In numerical experiments we study the properties of
different variants of MinRes analyzed in this paper.
One of the most challenging problems in dynamic concurrent multiscale
simulations is the reflectionless transfer of physical quantities between the
different scales. In particular, when coupling molecular dynamics and finite
element discretizations in solid body mechanics, often spurious wave reflections
are introduced by the applied coupling technique. The reflected waves are
typically of high frequency and are arguably of little importance in the domain
where the finite element discretization drives the simulation.
In this work, we provide an analysis of this phenomenon.
Based on the gained
insight, we derive a new coupling approach, which neatly separates high and low
frequency waves. Whereas low frequency waves are permitted to
bridge the scales, high frequency waves can be removed by applying damping techniques without affecting the coupled share of the solution. As a consequence, our new method almost completely eliminates unphysical wave reflections and deals in a consistent way with waves of arbitrary frequencies. The separation of
wavelengths is achieved by employing a discrete $L^2$-projection, which acts as a
low pass filter. Our coupling constraints enforce matching in the range of this projection. With respect to the numerical realization this approach
has the advantage of a small number of constraints, which is computationally
efficient. Numerical results in one and two dimensions confirm our theoretical
findings and illustrate the performance of our new weak coupling approach.
We investigate optimal elliptic
regularity (within the scale of Sobolev spaces) of anisotropic
div--grad operators in three dimensions at a multi-material vertex on
the Neumann boundary part of a polyhedral spatial domain. The
gradient of a solution to the corresponding elliptic PDE (in a
neighbourhood of the vertex) is integrable to an index greater than
three.
We consider Large Deformation Diffeomorphic Metric Mapping of general $m$-currents. After stating an optimization algorithm in the function space of admissable morph generating velocity fields, two innovative aspects in this framework are presented and numerically investigated: First, we spatially discretize the velocity field with conforming adaptive finite elements and discuss advantages of this new approach. Second, we directly compute the temporal evolution of discrete $m$-current attributes.
We consider the mechanical coupling of a geometrically exact Cosserat rod to a linear elastic continuum. The coupling conditions are formulated in the nonlinear rod configuration space. We describe a Dirichlet--Neumann algorithm for the coupled system, and use it to simulate the static stresses in a human knee joint, where the Cosserat rods are models for the ligaments.
In optimal control problems with nonlinear time-dependent 3D PDEs, full 4D discretizations are usually prohibitive due to the storage requirement. For this reason gradient and Newton type methods working on the reduced functional are often employed. The computation of the reduced gradient requires one solve of the state equation forward in time, and one backward solve of the adjoint equation. The state enters into the adjoint equation, again requiring the storage of a full 4D data set. We propose a lossy compression algorithm using an inexact but cheap predictor for the state data, with additional entropy coding of prediction errors. As the data is used inside a discretized, iterative algorithm, lossy compression
maintaining a certain error bound turns out to be sufficient.
We introduce geodesic finite elements as a conforming way to discretize partial differential equations for functions $v : \Omega \to M$, where $\Omega$ is an open subset of $\R^d$ and $M$ is a Riemannian manifold. These geodesic finite elements naturally generalize standard first-order
finite elements for Euclidean spaces. They also generalize the geodesic finite elements proposed for $d=1$ by the author. Our formulation is equivariant under isometries of $M$, and hence preserves objectivity of continuous problem formulations. We concentrate on partial differential equations that can be formulated as minimization problems. Discretization leads to algebraic minimization problems on product manifolds $M^n$. These can be solved efficiently
using a Riemannian trust-region method. We propose a monotone multigrid method to solve the constrained inner problems with linear multigrid speed. As an example we numerically compute harmonic maps from a domain in $\R^3$
to $S^2$.
Logical modeling of biological regulatory networks gives rise to a representation of the system's dynamics as a so-called state transition graph. Analysis of such a graph in its entirety allows for a comprehensive understanding of the functionalities and behavior of the modeled system. However, the size of the vertex set of the graph is exponential in the number of the network components making analysis costly, motivating development of reduction methods. In this paper, we present results allowing for a complete description of an asynchronous state transition graph of a Thomas network solely based on the analysis of the subgraph induced by certain extremal states. Utilizing this notion, we compare the behavior of a simple multi-valued network and a corresponding Boolean network and analyze the conservation of dynamical properties between them. Understanding the relation between such coarser and finer models is a necessary step towards meaningful network reduction as well as model refinement methods.