Past

Discussion of https://arxiv.org/abs/1911.07895

The task of solving large scale linear algebraic problems such as factorising matrices or solving linear systems is of central importance in many areas of scientific computing, as well as in data analysis and computational statistics. The talk will describe how randomisation can be used to design algorithms that in many environments have both better asymptotic complexities and better practical speed than standard deterministic methods.

The talk will in particular focus on randomised algorithms for solving large systems of linear equations. Both direct solution techniques based on fast factorisations of the coefficient matrix, and techniques based on randomised preconditioners, will be covered.

Note: There is a related talk in the Random Matrix Seminar on Tuesday Feb 25, at 15:30 in L4. That talk describes randomised methods for computing low rank approximations to matrices. The two talks are independent, but the Tuesday one introduces some of the analytical framework that supports the methods described here.

Sustainability is a highly complex topic, containing interwoven economic, ecological, and social aspects.  Simply defining the concept of sustainability is a challenge in itself.  Assessing the impact of sustainability efforts and generating effective policy requires analyzing the interactions and challenges presented by these different aspects. To address this challenge, it is necessary to develop methods that bridge fields and take into account phenomena that range from physical analysis of climate to network analysis of societal phenomena. In this talk, I will give an insight into areas of mathematical research that try to account for these inter-dependencies. The aim of this talk is to provide a critical discussion of the challenges in a joint discussion and emphasize the importance of multi-disciplinary approaches.

I will describe new solutions to the stationary Ginzburg-Landau equation in 3 dimensions with vortex lines given by interacting helices, with degree one around each filament and total degree an arbitrary positive integer. I will also present results on the asymptotic behavior of vortices in the entire plane for a dissipative Ginzburg-Landau equation. This is work in collaboration with Manuel del Pino, Remy Rodiac, Maria Medina, Monica Musso and Juncheng Wei.

We will provide a general overview on some recent work on non-archimedean parametrizations and their applications. We will start by presenting our work with Cluckers and Comte on the existence of good Yomdin-Gromov parametrizations in the non-archimedean context and a $p$-adic Pila-Wilkie theorem.   We will then explain how this is used in our work with Chambert-Loir to prove bialgebraicity results in products of Mumford curves. 
 

There is a class $\mathcal{B}$ of analytic Besov functions on a half-plane, with a very simple description.   This talk will describe a bounded functional calculus $f \in \mathcal{B} \mapsto f(A)$ where $-A$ is the generator of either a bounded $C_0$-semigroup on Hilbert space or a bounded analytic semigroup on a Banach space.    This calculus captures many known results for such operators in a unified way, and sometimes improves them.   A discrete version of the functional calculus was shown by Peller in 1983.

The talk will describe how ideas from random matrix theory can be leveraged to effectively, accurately, and reliably solve important problems that arise in data analytics and large scale matrix computations. We will focus in particular on accelerated techniques for computing low rank approximations to matrices. These techniques rely on randomised embeddings that reduce the effective dimensionality of intermediate steps in the computation. The resulting algorithms are particularly well suited for processing very large data sets.

The algorithms described are supported by rigorous analysis that depends on probabilistic bounds on the singular values of rectangular Gaussian matrices. The talk will briefly review some representative results.

Note: There is a related talk in the Computational Mathematics and Applications seminar on Thursday Feb 27, at 14:00 in L4. There, the ideas introduced in this talk will be extended to the problem of solving large systems of linear equations.

Robust principal component analysis and low-rank matrix completion are extensions of PCA that allow for outliers and missing entries, respectively. Solving these problems requires a low coherence between the low-rank matrix and the canonical basis. However, in both problems the well-posedness issue is even more fundamental; in some cases, both Robust PCA and matrix completion can fail to have any solutions due to the fact that the set of low-rank plus sparse matrices is not closed. Another consequence of this fact is that the lower restricted isometry property (RIP) bound cannot be satisfied for some low-rank plus sparse matrices unless further restrictions are imposed on the constituents. By restricting the energy of one of the components, we close the set and are able to derive the RIP over the set of low rank plus sparse matrices and operators satisfying concentration of measure inequalities. We show that the RIP of an operator implies exact recovery of a low-rank plus sparse matrix is possible with computationally tractable algorithms such as convex relaxations or line-search methods. We propose two efficient iterative methods called Normalized Iterative Hard Thresholding (NIHT) and Normalized Alternative Hard Thresholding (NAHT) that provably recover a low-rank plus sparse matrix from subsampled measurements taken by an operator satisfying the RIP.
 

Positively folded galleries arise as images of retractions of buildings onto a fixed apartment and play a role in many areas of maths (such as in the study of affine Hecke algebras, Macdonald polynomials, MV-polytopes, and affine Deligne-Lusztig varieties). In this talk, we will define positively folded galleries, and then look at how these can be used to study affine flag varieties. We will also look at a new recursive description of the set of end alcoves of folded galleries with respect to alcove-induced orientations, which gives us a combinatorial description of certain double coset intersections in these affine flag varieties. This talk is based on joint work with Elizabeth Milićević, Petra Schwer and Anne Thomas.

For a family $A$ in $\{0,...,k\}^n$, its deletion shadow is the set obtained from $A$ by deleting from any of its vectors one coordinate. Given the size of $A$, how should we choose $A$ to minimise its deletion shadow? And what happens if instead we may delete only a coordinate that is zero? We discuss these problems, and give an exact solution to the second problem.

Randomized SVD has become an extremely successful approach for efficiently computing a low-rank approximation of matrices. In particular the paper by Halko, Martinsson (who is speaking twice this week), and Tropp (SIREV 2011) contains extensive analysis, and made it a very popular method. 
The complexity for $m\times n$ matrices is $O(Nr+(m+n)r^2)$ where $N$ is the cost of a (fast) matrix-vector multiplication; which becomes $O(mn\log n+(m+n)r^2)$ for dense matrices. This work uses classical results in numerical linear algebra to reduce the computational cost to $O(Nr)$ without sacrificing numerical stability. The cost is essentially optimal for many classes of matrices, including $O(mn\log n)$ for dense matrices. The method can also be adapted for updating, downdating and perturbing the matrix, and is especially efficient relative to previous algorithms for such purposes.  

Stress perfusion cardiac magnetic resonance (CMR) imaging has been shown to be highly accurate for the detection of coronary artery disease. However, a major limitation is that the accuracy of the visual assessment of the images is challenging and thus the accuracy of the diagnosis is highly dependent on the training and experience of the reader. Quantitative perfusion CMR, where myocardial blood flow values are inferred directly from the MR images, is an automated and user-independent alternative to the visual assessment.

This talk will focus on addressing the main technical challenges which have hampered the adoption of quantitative myocardial perfusion MRI in clinical practice. The talk will cover the problem of respiratory motion in the images and the use of dimension reduction techniques, such as robust principal component analysis, to mitigate this problem. I will then discuss our deep learning-based image processing pipeline that solves the necessary series of computer vision tasks required for the blood flow modelling and introduce the Bayesian inference framework in which the kinetic parameter values are inferred from the imaging data.

A well-known theorem of Choquet-Bruhat and Geroch states that for given smooth initial data for the Einstein equations there exists a unique maximal globally hyperbolic development. In particular, time evolution of globally hyperbolic solutions is unique. This talk investigates whether the same result holds for quasilinear wave equations defined on a fixed background. After recalling the notion of global hyperbolicity, we first present an example of a quasilinear wave equation for which unique time evolution in fact fails and contrast this with the Einstein equations. We then proceed by presenting conditions on quasilinear wave equations which ensure uniqueness. This talk is based on joint work with Harvey Reall and Felicity Eperon.
 

Multilayer networks are a way to represent dependent connectivity patterns — e.g., time-dependence, multiple types of interactions, or both — that arise in many applications and which are difficult to incorporate into standard network representations. In the study of multilayer networks, it is important to investigate mesoscale (i.e., intermediate-scale) structures, such as communities, to discover features that lie between the microscale and the macroscale. We introduce a framework for the construction of generative models for mesoscale structure in multilayer networks.  We model dependency at the level of partitions rather than with respect to edges, and treat the process of generating a multilayer partition separately from the process of generating edges for a given multilayer partition. Our framework can admit many features of empirical multilayer networks and explicitly incorporates a user-specified interlayer dependency structure. We discuss the parameters and some properties of our framework, and illustrate an example of its use with benchmark models for multilayer community-detection tools. 

Spatial navigation in preclinical and clinical Alzheimer’s disease - Relevance for topological data analysis?

Spatial navigation changes are one of the first symptoms of Alzheimer’s disease and also lead to significant safeguarding issues in patients after diagnosis. Despite their significant implications, spatial navigation changes in preclinical and clinical Alzheimer’s disease are still poorly understood. In the current talk, I will explain the spatial navigation processes in the brain and their relevance to Alzheimer’s disease. I will then introduce our Sea Hero Quest project, which created the first global benchmark data for spatial navigation in ~4.5 million people worldwide via a VR-based game. I will present data from the game, which has allowed to create personalised benchmark data for at-risk-of-Alzheimer’s people. The final part of my talk will explore how real-world environment & entropy impacts on dementia patients getting lost and how this has relevance for GPS technology based safeguarding and car driving in Alzheimer’s disease.

Understanding the size of the rational points on a curve of higher genus is one of the major open problems in the theory of Diophantine equations. In this talk I will discuss the related problem of understanding how close together rational points can get. I will also discuss the relation to the subject of (generalised) Wieferich primes.

I will present a $\Gamma$-convergence for degenerate integral functionals related to homogenisation problems  in the Heisenberg group. In our  case, both the rescaling and the notion of invariance or periodicity are chosen in a way motivated by the geometry of the Heisenberg group. Without using special geometric features, these functionals would be neither coercive nor periodic, so classic results do not apply.  All the results apply to the more general case of Carnot groups. Joint with Nicolas Dirr, Paola Mannucci and Claudio Marchi.

The square peg problem asks whether every Jordan curve in the
plane contains the vertices of a square. Inspired by Hugelmeyer's approach
for smooth curves, we give a topological proof for "locally 1-Lipschitz"
curves using 4-dimensional topology.

We will talk about some stochastic parabolic and hyperbolic partial differential equations (SPDEs), which arise naturally in the context of Liouville quantum gravity. These dynamics are proposed to preserve the Liouville measure, which has been constructed recently in the series of works by David-Kupiainen-Rhodes-Vargas. We construct global solutions to these equations under some conditions and then show the invariance of the Liouville measure under the resulting dynamics. As a by-product, we also answer an open problem proposed by Sun-Tzvetkov recently.
 

Stochastic processes subject to weak noise often show a metastable
behaviour, meaning that they converge to equilibrium extremely slowly;
typically, the convergence time is exponentially large in the inverse
of the variance of the noise (Arrhenius law).
  
In the case of finite-dimensional Ito stochastic differential
equations, the large-deviation theory developed in the 1970s by
Freidlin and Wentzell allows to prove such Arrhenius laws and compute
their exponent. Sharper asymptotics for relaxation times, including the
prefactor of the exponential term (Eyring–Kramers laws) are known, for
instance, if the stochastic differential equation involves a gradient
drift term and homogeneous noise. One approach that has been very
successful in proving Eyring–Kramers laws, developed by Bovier,
Eckhoff, Gayrard and Klein around 2005, relies on potential theory.
  
I will describe Eyring–Kramers laws for some parabolic stochastic PDEs
such as the Allen–Cahn equation on the torus. In dimension 1, an
Arrhenius law was obtained in the 1980s by Faris and Jona-Lasinio,
using a large-deviation principle. The potential-theoretic approach
allows us to compute the prefactor, which turns out to involve a
Fredholm determinant. In dimensions 2 and 3, the equation needs to be
renormalized, which turns the Fredholm determinant into a
Carleman–Fredholm determinant.
  
Based on joint work with Barbara Gentz (Bielefeld), and with Ajay
Chandra (Imperial College), Giacomo Di Gesù (Vienna) and Hendrik Weber
(Warwick). 

References: 
https://dx.doi.org/10.1214/EJP.v18-1802
https://dx.doi.org/10.1214/17-EJP60

It is well-known that the Teichmüller space of a compact surface can be identified with a connected component of the moduli space of representations of the fundamental group of the surface in PSL(2,R). Higher Teichmüller components are generalizations of this that exist for the moduli space of representations of the fundamental group into certain real simple Lie groups of higher rank. As for the usual Teichmüller space, these components consist entirely  of discrete and faithful representations. Several cases have been identified over the years. First, the Hitchin components for split groups, then the maximal Toledo invariant components for Hermitian groups, and more recently certain components for SO(p,q). In this talk, I will describe a general construction of (still somewhat conjecturally) all possible Teichmüller components, and a parametrization of them using Higgs bundles.

String sigma-models relevant in the AdS/CFT correspondence are highly non-trivial two-dimensional field theories for which predictions at finite coupling exist, assuming integrability and/or the duality itself.  I will discuss general features of the perturbative approach to these models, and present progress on how to go extract finite coupling information in the most possibly general way, namely via the use of lattice field theory techniques. I will also present new results on certain ``defect-CFT’' correlators  at strong coupling.