Dirk Blömker

2papers

2 Papers

NAFeb 26, 2019
Well Posedness and Convergence Analysis of the Ensemble Kalman Inversion

Dirk Blömker, Claudia Schillings, Philipp Wacker et al.

The ensemble Kalman inversion is widely used in practice to estimate unknown parameters from noisy measurement data. Its low computational costs, straightforward implementation, and non-intrusive nature makes the method appealing in various areas of application. We present a complete analysis of the ensemble Kalman inversion with perturbed observations for a fixed ensemble size when applied to linear inverse problems. The well-posedness and convergence results are based on the continuous time scaling limits of the method. The resulting coupled system of stochastic differential equations allows to derive estimates on the long-time behaviour and provides insights into the convergence properties of the ensemble Kalman inversion. We view the method as a derivative free optimization method for the least-squares misfit functional, which opens up the perspective to use the method in various areas of applications such as imaging, groundwater flow problems, biological problems as well as in the context of the training of neural networks.

NANov 13, 2017
Numerically Computable A Posteriori-Bounds for stochastic Allen-Cahn equation

Dirk Blömker, Minoo Kamrani

The aim of this paper is the derivation of an a-posteriori error estimate for a numerical method based on an exponential scheme in time and spectral Galerkin methods in space. We obtain analytically a rigorous bound on the mean square error conditioned to the calculated data, which is numerically computable and uses the given numerical approximation. Thus one can check a-posteriori the error for a given numerical computation without relying on an asymptotic result. All estimates are only based on the numerical data and the structure of the equation, but they do not use any a-priori information of the solution, which makes the approach applicable to equations where global existence of solutions is not known. For simplicity of presentation, we develop the method here in a relatively simple situation of a stable one-dimensional Allen-Cahn equation with additive forcing.