Computational Mathematics and Applications Seminar

Professor Nick Trefethen will speak about: 'Quadrature = rational approximation'

Whenever you see a string of quadrature nodes, you can consider it as a branch cut defined by the poles of a rational approximation to the Cauchy transform of a weight function.  The aim of this talk is to explain this strange statement and show how it opens the way to calculation of targeted quadrature formulas for all kinds of applications.  Gauss quadrature is an example, but it is just the starting point, and many more examples will be shown.  I hope this talk will change your understanding of quadrature formulas. 

This is joint work with Andrew Horning. 
 

TBA

Speaker Estefania Loayza Romero will talk about:  A Riemannian Approach for PDE-Constrained Shape Optimization Using Outer Metrics

In PDE-constrained shape optimisation, shapes are traditionally viewed as elements of a Riemannian manifold, specifically as embeddings of the unit circle into the plane, modulo reparameterizations. The standard approach employs the Steklov-Poincaré metric to compute gradients for Riemannian optimisation methods. A significant limitation of current methods is the absence of explicit expressions for the geodesic equations associated with this metric. Consequently, algorithms have relied on retractions (often equivalent to the perturbation of identity method in shape optimisation) rather than true geodesic paths. Previous research suggests that incorporating geodesic equations, or better approximations thereof, can substantially enhance algorithmic performance. This talk presents numerical evidence demonstrating that using outer metrics, defined on the space of diffeomorphisms with known geodesic expressions, improves Riemannian gradient-based optimisation by significantly reducing the number of required iterations and preserving mesh quality throughout the optimisation process.

This talk is hosted at RAL. 

Professor Ana Djurdjevac will talk about; 'The Dean–Kawasaki Equation: Theory, Numerics, and Applications'

The Dean–Kawasaki equation provides a stochastic partial differential equation description of interacting particle systems at the level of empirical densities and has attracted considerable interest in statistical physics, stochastic analysis, and applied modeling. In this work, we study analytical and numerical aspects of the Dean–Kawasaki equation, with a particular focus on well-posedness, structure preservation, and possible discretization strategies. In addition, we extend the framework to the Dean–Kawasaki equation posed on smooth hypersurfaces. We discuss applications of the Dean–Kawasaki framework to particle-based models arising in biological systems and modeling social dynamics.

Speaker Yongho Park will talk about 'Subspace Correction Methods for Convex Optimization: Algorithms, Theory, and Applications'

This talk considers a framework of subspace correction methods for convex optimization, which provides a unified perspective for the design and analysis of a wide range of iterative methods, including advanced domain decomposition and multigrid methods. We first develop a convergence theory for parallel subspace correction methods based on the observation that these methods can be interpreted as nonlinearly preconditioned gradient descent methods. This viewpoint leads to a simpler and sharper analysis compared with existing approaches. We further show how the theory can be extended to semicoercive and nearly semicoercive problems. In addition, we explore connections between subspace correction methods and other classes of iterative algorithms, such as alternating projection methods, through the lens of convex duality, thereby enabling a unified treatment. Several applications are presented, including nonlinear partial differential equations, variational inequalities, and mathematical imaging problems. The talk concludes with a discussion of relevant and emerging research directions.

TBA

Speaker Katharina Schratz will talk about 'Resonances as a computational tool'

A large toolbox of numerical schemes for dispersive equations has been established, based on different discretization techniques such as discretizing the variation-of-constants formula (e.g., exponential integrators) or splitting the full equation into a series of simpler subproblems (e.g., splitting methods). In many situations these classical schemes allow a precise and efficient approximation. This, however, drastically changes whenever non-smooth phenomena enter the scene such as for problems at low regularity and high oscillations. Classical schemes fail to capture the oscillatory nature of the solution, and this may lead to severe instabilities and loss of convergence. In this talk I present a new class of resonance based schemes. The key idea in the construction of the new schemes is to tackle and deeply embed the underlying nonlinear  structure of resonances into the numerical discretization. As in the continuous case, these terms are central to structure preservation and offer the new schemes strong geometric properties at low regularity.

TBA

TBA; this talk is hosted at RAL. 

Placeholder entry; Fernando is visiting April-July so we'll ask him to speak some time in Trinity. 

TBA

Speaker Prof Dr Maria Lukacova will talk about 'Numerical analysis of oscillatory solutions of compressible flows'

Oscillatory solutions of compressible flows arise in many practical situations.  An iconic example is the Kelvin-Helmholtz problem, where standard numerical methods yield oscillatory solutions. In such a situation,  standard tools of numerical analysis for partial differential equations are not applicable. 

We will show that structure-preserving numerical methods converge in general to generalised solutions, the so-called dissipative solutions. 
The latter describes the limits of oscillatory sequences. We will concentrate on the inviscid flows, the Euler equations of gas dynamics, and mention also the relevant results obtained for the viscous compressible flows, governed by the Navier-Stokes equations.

We discuss a concept of K-convergence that turns a weak convergence of numerical solutions into the strong convergence of
their empirical means to a dissipative solution. The latter satisfies a weak formulation of the Euler equations modulo the Reynolds turbulent stress.  We will also discuss suitable selection criteria to recover well-posedness of the Euler equations of gas dynamics. Theoretical results will be illustrated by a series of numerical simulations.  

TBA

Speaker Professor Luis Nunes Vicente will talk about 'Reducing Sample Complexity in Stochastic Derivative-Free Optimization via Tail Bounds and Hypothesis Testing';

We introduce and analyze new probabilistic strategies for enforcing sufficient decrease conditions in stochastic derivative-free optimization, with the goal of reducing sample complexity and simplifying convergence analysis. First, we develop a new tail bound condition imposed on the estimated reduction in function value, which permits flexible selection of the power used in the sufficient decrease test, q in (1,2]. This approach allows us to reduce the number of samples per iteration from the standard O(delta^{−4}) to O(delta^{-2q}), assuming that the noise moment of order q/(q-1) is bounded. Second, we formulate the sufficient decrease condition as a sequential hypothesis testing problem, in which the algorithm adaptively collects samples until the evidence suffices to accept or reject a candidate step. This test provides statistical guarantees on decision errors and can further reduce the required sample size, particularly in the Gaussian noise setting, where it can approach O(delta^{−2-r}) when the decrease is of the order of delta^r. We incorporate both techniques into stochastic direct-search and trust-region methods for potentially non-smooth, noisy objective functions, and establish their global convergence rates and properties. 

This is joint work with Anjie Ding, Francesco Rinaldi, and Damiano Zeffiro.

TBA 

TBA

TBA