BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Robert I McLachlan

Robert I McLachlan appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NANumerical AnalysisMachine Learningmath-phmath.MPnlin.SIphysics.comp-ph

Catalog footprint

What is connected

11works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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.

Published work

11 published item(s)

preprint2022arXiv

Backward error analysis for variational discretisations of partial differential equations

In backward error analysis, an approximate solution to an equation is compared to the exact solution to a nearby modified equation. In numerical ordinary differential equations, the two agree up to any power of the step size. If the differential equation has a geometric property then the modified equation may share it. In this way, known properties of differential equations can be applied to the approximation. But for partial differential equations, the known modified equations are of higher order, limiting applicability of the theory. Therefore, we study symmetric solutions of discretized partial differential equations that arise from a discrete variational principle. These symmetric solutions obey infinite-dimensional functional equations. We show that these equations admit second-order modified equations which are Hamiltonian and also possess first-order Lagrangians in modified coordinates. The modified equation and its associated structures are computed explicitly for the case of rotating travelling waves in the nonlinear wave equation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Structure preserving deep learning

Over the past few years, deep learning has risen to the foreground as a topic of massive interest, mainly as a result of successes obtained in solving large-scale image processing tasks. There are multiple challenging mathematical problems involved in applying deep learning: most deep learning methods require the solution of hard optimisation problems, and a good understanding of the tradeoff between computational effort, amount of data and model complexity is required to successfully design a deep learning approach for a given problem. A large amount of progress made in deep learning has been based on heuristic explorations, but there is a growing effort to mathematically understand the structure in existing deep learning methods and to systematically design new deep learning methods to preserve certain types of structure in deep learning. In this article, we review a number of these directions: some deep neural networks can be understood as discretisations of dynamical systems, neural networks can be designed to have desirable properties such as invertibility or group equivariance, and new algorithmic frameworks based on conformal Hamiltonian systems and Riemannian manifolds to solve the optimisation problems have been proposed. We conclude our review of each of these topics by discussing some open problems that we consider to be interesting directions for future research.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Detection of high codimensional bifurcations in variational PDEs

We derive bifurcation test equations for A-series singularities of nonlinear functionals and, based on these equations, we propose a numerical method for detecting high codimensional bifurcations in parameter-dependent PDEs such as parameter-dependent semilinear Poisson equations. As an example, we consider a Bratu-type problem and show how high codimensional bifurcations such as the swallowtail bifurcation can be found numerically. In particular, our original contributions are (1) the use of the Infinite-dimensional Splitting Lemma, (2) the unified and simplified treatment of all A-series bifurcations, (3) the presentation in Banach spaces, i.e. our results apply both to the PDE and its (variational) discretization, (4) further simplifications for parameter-dependent semilinear Poisson equations (both continuous and discrete), and (5) the unified treatment of the continuous problem and its discretisation.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Travelling wave solutions of multisymplectic discretizations of semi-linear wave equations

How well do multisymplectic discretisations preserve travelling wave solutions? To answer this question, the 5-point central difference scheme is applied to the semi-linear wave equation. A travelling wave ansatz leads to an ordinary difference equation, whose solutions correspond to the numerical scheme and can be compared to travelling wave solutions of the corresponding PDE. For a discontinuous nonlinearity the difference equation is solved exactly. For continuous nonlinearities the difference equation is solved using a Fourier series, and resonances that depend on the grid-size are revealed for a smooth nonlinearity. In general, the infinite dimensional functional equation, which must be solved to get the travelling wave solutions, is intractable, but backward error analysis proves to be a powerful tool, as it provides a way to study the solutions of the equation through a simple ODE that describes the behavior to arbitrarily high order. A general framework for using backward error analysis to analyze preservation of travelling waves for other equations and discretisations is presented. Then, the advantages that multisymplectic methods have over other methods are briefly highlighted.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Collective symplectic integrators

We construct symplectic integrators for Lie-Poisson systems. The integrators are standard symplectic (partitioned) Runge--Kutta methods. Their phase space is a symplectic vector space with a Hamiltonian action with momentum map $J$ whose range is the target Lie--Poisson manifold, and their Hamiltonian is collective, that is, it is the target Hamiltonian pulled back by $J$. The method yields, for example, a symplectic midpoint rule expressed in 4 variables for arbitrary Hamiltonians on $\mathfrak{so}(3)^*$. The method specializes in the case that a sufficiently large symmetry group acts on the fibres of $J$, and generalizes to the case that the vector space carries a bifoliation. Examples involving many classical groups are presented.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Discrete gradient methods have an energy conservation law

We show for a variety of classes of conservative PDEs that discrete gradient methods designed to have a conserved quantity (here called energy) also have a time-discrete conservation law. The discrete conservation law has the same conserved density as the continuous conservation law, while its flux is found by replacing all derivatives of the conserved density appearing in the continuous flux by discrete gradients.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Geometric properties of Kahan's method

We show that Kahan's discretization of quadratic vector fields is equivalent to a Runge--Kutta method. When the vector field is Hamiltonian on either a symplectic vector space or a Poisson vector space with constant Poisson structure, the map determined by this discretization has a conserved modified Hamiltonian and an invariant measure, a combination previously unknown amongst Runge--Kutta methods applied to nonlinear vector fields. This produces large classes of integrable rational mappings in two and three dimensions, explaining some of the integrable cases that were previously known.

0 citations0 reviewsNo code linkedOpen access