Tue, 28 Jan 2014

14:30 - 15:00
L5

An algorithm for the convolution of Legendre expansions

Nick Hale
(University of Oxford)
Abstract

Convolution is widely-used and fundamental mathematical operation
in signal processing, statistics, and PDE theory.

Unfortunately the CONV() method in Chebfun for convolving two chebfun 
objects has long been one of the most disappointingly slow features of 
the project. In this talk we will present a new algorithm, which shows 
performance gains on the order of a factor 100.

The key components of the new algorithm are:
* a convolution theorem for Legendre polynomials 
* recurrence relations satisfied by spherical Bessel functions
* recent developments in fast Chebyshev-Legendre transforms [1]

Time-permitting, we shall end with an application from statistics,
using the fact that the probability distribution of the sum of two 
independent random variables is the convolution of their individual 
distributions.

[1] N. Hale and A. Townsend, "A fast, simple, and stable Chebyshev-
Legendre transform using an asymptotic formula”, SISC (to appear).

Mon, 24 Feb 2014

14:15 - 15:15
Eagle House

The splitting method for SPDEs: from robustness to applications in financial engineering, nonlinear filtering and optimal control

HARALD OBERHAUSER
(University of Oxford)
Abstract

The splitting-up method is a powerful tool to solve (SP)DEs by dividing the equation into a set of simpler equations that are easier to handle. I will speak about how such splitting schemes can be derived and extended by insights from the theory of rough paths.

Finally, I will discuss numerics for real-world applications that appear in the management of risk and engineering applications like nonlinear filtering.

Mon, 10 Feb 2014

14:15 - 15:15
Eagle House

Discretely sampled signals and the rough Hoff path

GUY FLINT
(University of Oxford)
Abstract

Sampling a $d$-dimensional continuous signal (say a semimartingale) $X:[0,T] \rightarrow \mathbb{R}^d$ at times $D=(t_i)$, we follow the recent papers [Gyurko-Lyons-Kontkowski-Field-2013] and [Lyons-Ni-Levin-2013] in constructing a lead-lag path; to be precise, a piecewise-linear, axis-directed process $X^D: [0,1] \rightarrow
\mathbb{R}^{2d}$ comprised of a past and future component. Lifting $X^D$ to its natural rough path enhancement, we can consider the question of convergence as
the latency of our sampling becomes finer.

Mon, 03 Feb 2014

14:15 - 15:15
Eagle House

TBC

DANYU YANG
(University of Oxford)
Tue, 21 Jan 2014

14:00 - 15:00
L5

Numerical solution of Hamilton—Jacobi—Bellman equations

Iain Smears
(University of Oxford)
Abstract

Hamilton—Jacobi—Bellman (HJB) equations are a class of fully nonlinear second-order partial differential equations (PDE) of elliptic or parabolic type that originate from Stochastic Optimal Control Theory. These PDE are fully nonlinear in the sense that the nonlinear terms include the second partial derivatives of the unknown solution; this strong nonlinearity severely restricts the range of numerical methods that are known to be convergent. These problems have traditionally been solved with low order monotone schemes, often of finite difference type, which feature certain limitations in terms of efficiency and practicability.
In this summary talk of my DPhil studies, we will be interested in the development of hp-version discontinuous Galerkin finite element methods (DGFEM) for the class of HJB equations that satisfy a Cordès condition. First, we will show the novel techniques of analysis used to find a stable and convergent scheme in the elliptic setting, and then we will present recent work on their extension to parabolic problems. The resulting method is very nonstandard, provably of high order, and it even allows for exponential convergence under hp-refinement. We present numerical experiments showing the accuracy, computational efficiency and flexibility of the scheme
Thu, 13 Feb 2014

12:00 - 13:00
L6

Modelling collective motion in biology

Prof. Philip Maini
(University of Oxford)
Abstract

We will present three different recent applications of cell motion in biology: (i) Movement of epithelial sheets and rosette formation, (ii) neural crest cell migrations, (iii) acid-mediated cancer cell invasion. While the talk will focus primarily on the biological application, it will be shown that all of these processes can be represented by reaction-diffusion equations with nonlinear diffusion term.

Tue, 26 Nov 2013

14:30 - 15:00
L5

Small dot, big challenging: on the new benchmark of Top500 and Green500

Shengxin (Jude) Zhu
(University of Oxford)
Abstract

A new benchmark, High Performance Conjugate Gradient (HPCG), finally was introduced recently for the Top500 list and the Green500 list. This will draw more attention to performance of sparse iterative solvers on distributed supercomputers and energy efficiency of hardware and software. At the same time, this will more widely promote the concept that communications are the bottleneck of performance of iterative solvers on distributed supercomputers, here we will go a little deeper, discussing components of communications and discuss which part takes a dominate share. Also discussed are mathematics tricks to detect some metrics of an underlying supercomputer.

Subscribe to University of Oxford