Variational Ensemble Filters for Sequential Inverse Problems
Abstract
Given a model dynamical system, a model of any measuring apparatus relating states to observations, and a prior assessment of uncertainty, the probability density of subsequent system states, conditioned upon the history of the observations, is of some practical interest.
When observations are made at discrete times, it is known that the evolving probability density is a solution of the Bayesian filtering equations. This talk will describe the difficulties in approximating the evolving probability density using a Gaussian mixture (i.e. a sum of Gaussian densities). In general this leads to a sequence of optimisation problems and related high-dimensional integrals. There are other problems too, related to the necessity of using a small number of densities in the mixture, the requirement to maintain sparsity of any matrices and the need to compute first and, somewhat disturbingly, second derivatives of the misfit between predictions and observations. Adjoint methods, Taylor expansions, Gaussian random fields and Newton’s method can be combined to, possibly, provide a solution. The approach is essentially a combination of filtering methods and '4-D Var’ methods and some recent progress will be described.
What is the mathematics of the Faraday cage?
Abstract
Everybody has heard of the Faraday cage effect, in which a wire mesh does a good job of blocking electric fields and electromagnetic waves. For example, the screen on the front of your microwave oven keeps the microwaves from getting out, while light with its smaller wavelength escapes so you can see your burrito. Surely the mathematics of such a famous and useful phenomenon has been long ago worked out and written up in the physics books, right?
Well, maybe. Dave Hewett and I have communicated with dozens of mathematicians, physicists, and engineers on this subject so far, and we've turned up amazingly little. Everybody has a view of why the Faraday cage mathematics is obvious, and most of their views are different. Feynman discusses the matter in his Lectures on Physics, but so far as we can tell, he gets it wrong.
For the static case at least (the Laplace equation), Hewett and I have made good progress with numerical explorations based on Mikhlin's method backed up by a theorem. The effect seems to much weaker than we had imagined -- are we missing something? For time-harmonic waves (the Helmholtz equation), our simulations lead to further puzzles. We need advice! Where in the world is the literature on this problem?
Integral representation of martingales motivated by the problem of market completion with derivative securities.
Abstract
A model of a financial market is complete if any payoff can be obtained as the terminal value of a self-financing trading strategy. It is well known that numerous models, for example stochastic volatility models, are however incomplete. We present conditions, which, in a general diffusion framework, guarantee that in such cases the market of primitive assets enlarged with an appropriate number of traded derivative contracts is complete. From a purely mathematical point of view we prove an integral representation theorem which guarantees that every local Q-martingale can be represented as a stochastic integral with respect to the vector of primitive assets and derivative contracts.
We consider the short time asymptotics of the heat content $E(s)$ of a domain $D$ of $\mathbb{R}^d$, where $D$ has a random boundary.
Abstract
When $\partial D$ is spatially homogeneous, we show that we can recover the lower and upper Minkowski dimensions of $\partial D$ from the sort time behaviour of $E(s)$. Furthermore, when the Minkowski dimension exists, finer geometric fluctuations can be recovered and $E(s)$ is controlled by $s^\alpha e^{f(\log(1/s))}$ for small $s$, for some $\alpha \in (0, \infty)$ and some regularly varying function $f$. The function $f$ is not constant is general and carries some geometric information.
When $\partial D$ is statistically self-similar, the Minkowski dimension and content of $\partial D$ typically exist and can be recovered from $E(s)$. Furthermore, $E(s)$ has an almost sure expansion $E(s) = c s^{\alpha} N_\infty + o(s^\alpha)$ for small $s$, for some $c$ and $\alpha \in (0, \infty)$ and some positive random variable $N_\infty$ with unit expectation arising as the limit of some martingale. In some cases, we can show that the fluctuations around this almost sure behaviour are governed by a central limit theorem, and conjecture that this is true more generally.
This is based on joint work with David Croydon and Ben Hambly.
The decay rate of the expected signature of a stopped Brownian motion
Abstract
In this presentation, we focus on the decay rate of the expected signature of a stopped Brownian motion; more specifically we consider two types of the stopping time: the first one is the Brownian motion up to the first exit time from a bounded domain $\Gamma$, denoted by $\tau_{\Gamma}$, and the other one is the Brownian motion up to $min(t, \tau_{\Gamma\})$. For the first case, we use the Sobolev theorem to show that its expected signature is geometrically bounded while for the second case we use the result in paper (Integrability and tail estimates for Gaussian rough differential equation by Thomas Cass, Christian Litterer and Terry Lyons) to show that each term of the expected signature has the decay rate like 1/ \sqrt((n/p)!) where p>2. The result for the second case can imply that its expected signature determines the law of the signature according to the paper (Unitary representations of geometric rough paths by Ilya Chevyrev)
The algebraicity of sieved sets and rational points on curves
Abstract
Using multiple frequencies to satisfy local constraints in PDE and applications to hybrid inverse problems
Abstract
In this talk I will describe a multiple frequency approach to the boundary control of Helmholtz and Maxwell equations. We give boundary conditions and a finite number of frequencies such that the corresponding solutions satisfy certain non-zero constraints inside the domain. The suitable boundary conditions and frequencies are explicitly constructed and do not depend on the coefficients, in contrast to the illuminations given as traces of complex geometric optics solutions. This theory finds applications in several hybrid imaging modalities. Some examples will be discussed.
Particle Methods for Inference in Non-linear Non-Gaussian State-Space Models
Abstract
State-space models are a very popular class of time series models which have found thousands of applications in engineering, robotics, tracking, vision, econometrics etc. Except for linear and Gaussian models where the Kalman filter can be used, inference in non-linear non-Gaussian models is analytically intractable. Particle methods are a class of flexible and easily parallelizable simulation-based algorithms which provide consistent approximations to these inference problems. The aim of this talk is to introduce particle methods and to present the most recent developments in this area.
Towards realistic performance for iterative methods on shared memory machines
Abstract
This talk introduces a random linear model to investigate the memory bandwidth barrier effect on current shared memory computers. Based on the fact that floating-point operations can be hidden by implicit compiling techniques, the runtime for memory intensive applications can be modelled by memory reference time plus a random term. The random term due to cache conflicts, data reuse and other environmental factors is proportional to memory reference volume. Statistical techniques are used to quantify the random term and the runtime performance parameters. Numerical results based on thousands representative matrices from various applications are presented, compared, analysed and validated to confirm the proposed model. The model shows that a realistic and fair metric for performance of iterative methods and other memory intensive applications should consider the memory bandwidth capability and memory efficiency.