Source author record

David Kelly

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record

math.PR math.NA Distributed, Parallel, and Cluster Computing math.AP math.CA math.DS Quantitative Methods Software Engineering

Catalog footprint

Research overview

8works

8topics

4close collaborators

Actions

Is this your work?

Open graph Browse papers

preprint2016arXiv

Fluctuations in the heterogeneous multiscale methods for fast-slow systems

How heterogeneous multiscale methods (HMM) handle fluctuations acting on the slow variables in fast-slow systems is investigated. In particular, it is shown via analysis of central limit theorems (CLT) and large deviation principles (LDP) that the standard version of HMM artificially amplifies these fluctuations. A simple modification of HMM, termed parallel HMM, is introduced and is shown to remedy this problem, capturing fluctuations correctly both at the level of the CLT and the LDP. Similar type of arguments can also be used to justify that the tau-leaping method used in the context of Gillespie's stochastic simulation algorithm for Markov jump processes also captures the right CLT and LDP for these processes.

preprint2016arXiv

Rough path recursions and diffusion approximations

In this article, we consider diffusion approximations for a general class of stochastic recursions. Such recursions arise as models for population growth, genetics, financial securities, multiplicative time series, numerical schemes and MCMC algorithms. We make no particular probabilistic assumptions on the type of noise appearing in these recursions. Thus, our technique is well suited to recursions where the noise sequence is not a semi-martingale, even though the limiting noise may be. Our main theorem assumes a weak limit theorem on the noise process appearing in the random recursions and lifts it to diffusion approximation for the recursion itself. To achieve this, we approximate the recursion (pathwise) by the solution to a stochastic equation driven by piecewise smooth paths; this can be thought of as a pathwise version of backward error analysis for SDEs. We then identify the limit of this stochastic equation, and hence the original recursion, using tools from rough path theory. We provide several examples of diffusion approximations, both new and old, to illustrate this technique.

preprint2016arXiv

Smooth approximation of stochastic differential equations

Consider an Itô process $X$ satisfying the stochastic differential equation $dX=a(X)\,dt+b(X)\,dW$ where $a,b$ are smooth and $W$ is a multidimensional Brownian motion. Suppose that $W_n$ has smooth sample paths and that $W_n$ converges weakly to $W$ . A central question in stochastic analysis is to understand the limiting behavior of solutions $X_n$ to the ordinary differential equation $dX_n=a(X_n)\,dt+b(X_n)\,dW_n$ . The classical Wong--Zakai theorem gives sufficient conditions under which $X_n$ converges weakly to $X$ provided that the stochastic integral $\int b(X)\,dW$ is given the Stratonovich interpretation. The sufficient conditions are automatic in one dimension, but in higher dimensions the correct interpretation of $\int b(X)\,dW$ depends sensitively on how the smooth approximation $W_n$ is chosen. In applications, a natural class of smooth approximations arise by setting $W_n(t)=n^{-1/2}\int_0^{nt}v\circϕ_s\,ds$ where $ϕ_t$ is a flow (generated, e.g., by an ordinary differential equation) and $v$ is a mean zero observable. Under mild conditions on $ϕ_t$ , we give a definitive answer to the interpretation question for the stochastic integral $\int b(X)\,dW$ . Our theory applies to Anosov or Axiom A flows $ϕ_t$ , as well as to a large class of nonuniformly hyperbolic flows (including the one defined by the well-known Lorenz equations) and our main results do not require any mixing assumptions on $ϕ_t$ . The methods used in this paper are a combination of rough path theory and smooth ergodic theory.

preprint2015arXiv

Nonlinear stability of the ensemble Kalman filter with adaptive covariance inflation

The Ensemble Kalman filter and Ensemble square root filters are data assimilation methods used to combine high dimensional nonlinear models with observed data. These methods have proved to be indispensable tools in science and engineering as they allow computationally cheap, low dimensional ensemble state approximation for extremely high dimensional turbulent forecast models. From a theoretical perspective, these methods are poorly understood, with the exception of a recently established but still incomplete nonlinear stability theory. Moreover, recent numerical and theoretical studies of catastrophic filter divergence have indicated that stability is a genuine mathematical concern and can not be taken for granted in implementation. In this article we propose a simple modification of ensemble based methods which resolves these stability issues entirely. The method involves a new type of adaptive covariance inflation, which comes with minimal additional cost. We develop a complete nonlinear stability theory for the adaptive method, yielding Lyapunov functions and geometric ergodicity under weak assumptions. We present numerical evidence which suggests the adaptive methods have improved accuracy over standard methods and completely eliminate catastrophic filter divergence. This enhanced stability allows for the use of extremely cheap, unstable forecast integrators, which would otherwise lead to widespread filter malfunction.

preprint2014arXiv

Geometric versus non-geometric rough paths

In this article we consider rough differential equations (RDEs) driven by non-geometric rough paths, using the concept of branched rough paths introduced in Gubinelli (2004). We first show that branched rough paths can equivalently be defined as $γ$ -Hölder continuous paths in some Lie group, akin to geometric rough paths. We then show that every branched rough path can be encoded in a geometric rough path. More precisely, for every branched rough path $\mathbf{X}$ lying above a path $X$ , there exists a geometric rough path $\bar{\mathbf{X}}$ lying above an extended path $\bar X$ , such that $\bar{\mathbf{X}}$ contains all the information of $\mathbf{X}$ . As a corollary of this result, we show that every RDE driven by a non-geometric rough path $\mathbf{X}$ can be rewritten as an extended RDE driven by a geometric rough path $\bar{\mathbf{X}}$ . One could think of this as a generalisation of the Itô-Stratonovich correction formula.

preprint2013arXiv

Reusability in Science: From Initial User Engagement to Dissemination of Results

Effective use of parallel and distributed computing in science depends upon multiple interdependent entities and activities that form an ecosystem. Active engagement between application users and technology catalysts is a crucial activity that forms an integral part of this ecosystem. Technology catalysts play a crucial role benefiting communities beyond a single user group. An effective user-engagement, use and reuse of tools and techniques has a broad impact on software sustainability. From our experience, we sketch a life-cycle for user-engagement activity in scientific computational environment and posit that application level reusability promotes software sustainability. We describe our experience in engaging two user groups from different scientific domains reusing a common software and configuration on different computational infrastructures.

preprint2012arXiv

Stochastic PDEs with multiscale structure

We study the spatial homogenisation of parabolic linear stochastic PDEs exhibiting a two-scale structure both at the level of the linear operator and at the level of the Gaussian driving noise. We show that in some cases, in particular when the forcing is given by space-time white noise, it may happen that the homogenised SPDE is not what one would expect from existing results for PDEs with more regular forcing terms.

preprint2010arXiv

Inferring hidden Markov models from noisy time sequences: a method to alleviate degeneracy in molecular dynamics

We present a new method for inferring hidden Markov models from noisy time sequences without the necessity of assuming a model architecture, thus allowing for the detection of degenerate states. This is based on the statistical prediction techniques developed by Crutchfield et al., and generates so called causal state models, equivalent to hidden Markov models. This method is applicable to any continuous data which clusters around discrete values and exhibits multiple transitions between these values such as tethered particle motion data or Fluorescence Resonance Energy Transfer (FRET) spectra. The algorithms developed have been shown to perform well on simulated data, demonstrating the ability to recover the model used to generate the data under high noise, sparse data conditions and the ability to infer the existence of degenerate states. They have also been applied to new experimental FRET data of Holliday Junction dynamics, extracting the expected two state model and providing values for the transition rates in good agreement with previous results and with results obtained using existing maximum likelihood based methods.

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

David Kelly

Research overview

Is this your work?

8 publications

Fluctuations in the heterogeneous multiscale methods for fast-slow systems

Rough path recursions and diffusion approximations

Smooth approximation of stochastic differential equations

Nonlinear stability of the ensemble Kalman filter with adaptive covariance inflation

Geometric versus non-geometric rough paths

Reusability in Science: From Initial User Engagement to Dissemination of Results

Stochastic PDEs with multiscale structure

Inferring hidden Markov models from noisy time sequences: a method to alleviate degeneracy in molecular dynamics

See the researcher in context

Building this map preview