2,025 research outputs found
Multilevel quadrature for elliptic problems on random domains by the coupling of FEM and BEM
Elliptic boundary value problems which are posed on a random domain can be
mapped to a fixed, nominal domain. The randomness is thus transferred to the
diffusion matrix and the loading. While this domain mapping method is quite
efficient for theory and practice, since only a single domain discretisation is
needed, it also requires the knowledge of the domain mapping.
However, in certain applications, the random domain is only described by its
random boundary, while the quantity of interest is defined on a fixed,
deterministic subdomain. In this setting, it thus becomes necessary to compute
a random domain mapping on the whole domain, such that the domain mapping is
the identity on the fixed subdomain and maps the boundary of the chosen fixed,
nominal domain on to the random boundary.
To overcome the necessity of computing such a mapping, we therefore couple
the finite element method on the fixed subdomain with the boundary element
method on the random boundary. We verify the required regularity of the
solution with respect to the random domain mapping for the use of multilevel
quadrature, derive the coupling formulation, and show by numerical results that
the approach is feasible
Fast derivatives of likelihood functionals for ODE based models using adjoint-state method
We consider time series data modeled by ordinary differential equations
(ODEs), widespread models in physics, chemistry, biology and science in
general. The sensitivity analysis of such dynamical systems usually requires
calculation of various derivatives with respect to the model parameters.
We employ the adjoint state method (ASM) for efficient computation of the
first and the second derivatives of likelihood functionals constrained by ODEs
with respect to the parameters of the underlying ODE model. Essentially, the
gradient can be computed with a cost (measured by model evaluations) that is
independent of the number of the ODE model parameters and the Hessian with a
linear cost in the number of the parameters instead of the quadratic one. The
sensitivity analysis becomes feasible even if the parametric space is
high-dimensional.
The main contributions are derivation and rigorous analysis of the ASM in the
statistical context, when the discrete data are coupled with the continuous ODE
model. Further, we present a highly optimized implementation of the results and
its benchmarks on a number of problems.
The results are directly applicable in (e.g.) maximum-likelihood estimation
or Bayesian sampling of ODE based statistical models, allowing for faster, more
stable estimation of parameters of the underlying ODE model.Comment: 5 figure
Variational Integrators for Nonvariational Partial Differential Equations
Variational integrators for Lagrangian dynamical systems provide a systematic
way to derive geometric numerical methods. These methods preserve a discrete
multisymplectic form as well as momenta associated to symmetries of the
Lagrangian via Noether's theorem. An inevitable prerequisite for the derivation
of variational integrators is the existence of a variational formulation for
the considered problem. Even though for a large class of systems this
requirement is fulfilled, there are many interesting examples which do not
belong to this class, e.g., equations of advection-diffusion type frequently
encountered in fluid dynamics or plasma physics. On the other hand, it is
always possible to embed an arbitrary dynamical system into a larger Lagrangian
system using the method of formal (or adjoint) Lagrangians. We investigate the
application of the variational integrator method to formal Lagrangians, and
thereby extend the application domain of variational integrators to include
potentially all dynamical systems. The theory is supported by physically
relevant examples, such as the advection equation and the vorticity equation,
and numerically verified. Remarkably, the integrator for the vorticity equation
combines Arakawa's discretisation of the Poisson brackets with a symplectic
time stepping scheme in a fully covariant way such that the discrete energy is
exactly preserved. In the presentation of the results, we try to make the
geometric framework of variational integrators accessible to non specialists.Comment: 49 page
The complex step method for approximating the Fréchet derivative of matrix functions in automorphism groups
We show, that the Complex Step approximation to the Fréchet derivative of matrix functions is applicable to the matrix sign, square root and polar mapping using iterative schemes. While this property was already discovered for the matrix sign using Newtons method, we extend the research to the family of Padé iterations, that allows us to introduce iterative schemes for finding function and derivative values while approximately preserving automorphism group structure
- …