University of Oxford

We propose several optimal preconditioners for systems defined by some functions $g$ of Toeplitz matrices $T_n$. In this paper we are interested in solving $g(T_n)x=b$ by the preconditioned conjugate method or the preconditioned minimal residual method, namely in the cases when $g(T_n)$ are the analytic functions $e^{T_n}$, $\sin{T_n}$ and $\cos{T_n}$. Numerical results are given to show the effectiveness of the proposed preconditioners.

Kerdock matrices are an attractive choice as deterministic measurement matrices for compressive sensing. I'll explain how Kerdock matrices are constructed, and then show how they can be adapted to one particular  strategy for quantizing measurements, in which measurements exceeding the desired dynamic range are rejected.

Efficient factorization or efficient computation of class 
numbers would both suffice to break RSA.  However the talk lies more in 
computational number theory rather than in cryptography proper. We will 
address two questions: (1) How quickly can one construct a factor table 
for the numbers up to x?, and (2) How quickly can one do the same for the 
class numbers (of imaginary quadratic fields)? Somewhat surprisingly, the 
approach we describe for the second problem is motivated by the classical 
Hardy-Littlewood method.

Operators, functions, and functionals are combined in many problems of computational science in a fashion that has the same logical structure as is familiar for block matrices and vectors.  It is proposed that the explicit consideration of such block structures at the continuous as opposed to discrete level can be a useful tool.  In particular, block operator diagrams provide templates for spectral discretization by the rectangular differentiation, integration, and identity matrices introduced by Driscoll and Hale.  The notion of the rectangular shape of a linear operator can be made rigorous by the theory of Fredholm operators and their indices, and the block operator formulations apply to nonlinear problems too, where the convention is proposed of representing nonlinear blocks as shaded.  At each step of a Newton iteration, the structure is linearized and the blocks become unshaded, representing Fréchet derivative operators, square or rectangular.  The use of block operator diagrams makes it possible to precisely specify discretizations of even rather complicated problems with just a few lines of pseudocode.

[Joint work with Nick Trefethen]

The Faraday cage effect is the phenomenon whereby electrostatic and electromagnetic fields are shielded by a wire mesh "cage". Nick Trefethen, Jon Chapman and I recently carried out a mathematical analysis of the two-dimensional electrostatic problem with thin circular wires, demonstrating that the shielding effect is not as strong as one might infer from the physics literature. In this talk I will present new results generalising the previous analysis to the electromagnetic case, and to wires of arbitrary shape. The main analytical tool is the asymptotic method of multiple scales, which is used to derive continuum models for the shielding, involving homogenized boundary conditions on an effective cage boundary. In the electromagnetic case one observes interesting resonance effects, whereby at frequencies close to the natural frequencies of the equivalent solid shell, the presence of the cage actually amplifies the incident field, rather than shielding it. We discuss applications to radiation containment in microwave ovens and acoustic scattering by perforated shells. This is joint work with Ian Hewitt.

Accurate methods for the first-order advection equation, used for example in tracking contaminants in fluids, usually exploit the theory of characteristics. Such methods are described and contrasted with methods that do not make use of characteristics.

Then the second-order wave equation, in the form of a first-order system, is considered. A review of the one-dimensional theory using solutions of various Riemann problems will be provided. In the special case that the medium has the ‘Goupillaud’ property, that waves take the same time to travel through each layer, one can derive exact solutions even when the medium is spatially heterogeneous. The extension of this method to two-dimensional problems will then be discussed. In two-dimensions it is not apparent that exact solutions can be found, however by exploiting a generalised Goupillaud property, it is possible to calculate approximate solutions of high accuracy, perhaps sufficient to be of benchmark quality. Some two-dimensional simulations, using exact one-dimensional solutions and operator splitting, will be described and a numerical evaluation of accuracy will be given.

In this talk I will review mathematical models used to describe the dynamics of ice sheets, and highlight some current areas of active research.  Melting of glaciers and ice sheets causes an increase in global sea level, and provides many other feedbacks on isostatic adjustment, the dynamics of the ocean, and broader climate patterns.  The rate of melting has increased in recent years, but there is still considerable uncertainty over  why this is, and whether the increase will continue.  Central to these questions is understanding the physics of how the ice intereacts with the atmosphere, the ground on which it rests, and with the ocean at its margins.  I will given an overview of the fluid mechanical problems involved and the current state of mathematical/computational modelling.  I will focus particularly on the issue of changing lubrication due to water flowing underneath the ice, and discuss how we can use models to rationalise observations of ice speed-up and slow-down.