Source author record

Gabriel Lord

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

ResearcherUnclaimed source record

math.NA math.SP Numerical Analysis math.CA math.DS math.FA q-fin.CP

Catalog footprint

What is connected

7works

7topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2022arXiv

The role of adaptivity in a numerical method for the Cox-Ingersoll-Ross model

We demonstrate the effectiveness of an adaptive explicit Euler method for the approximate solution of the Cox-Ingersoll-Ross model. This relies on a class of path-bounded timestepping strategies which work by reducing the stepsize as solutions approach a neighbourhood of zero. The method is hybrid in the sense that a convergent backstop method is invoked if the timestep becomes too small, or to prevent solutions from overshooting zero and becoming negative. Under parameter constraints that imply Feller's condition, we prove that such a scheme is strongly convergent, of order at least 1/2. Control of the strong error is important for multi-level Monte Carlo techniques. Under Feller's condition we also prove that the probability of ever needing the backstop method to prevent a negative value can be made arbitrarily small. Numerically, we compare this adaptive method to fixed step implicit and explicit schemes, and a novel semi-implicit adaptive variant. We observe that the adaptive approach leads to methods that are competitive in a domain that extends beyond Feller's condition, indicating suitability for the modelling of stochastic volatility in Heston-type asset models.

preprint2020arXiv

Selection of quasi-stationary states in the stochastically forced Navier-Stokes equation on the torus

The stochastically forced vorticity equation associated with the two dimensional incompressible Navier-Stokes equation on $D_δ:=[0,2πδ]\times [0,2π]$ is considered for $δ\approx 1$, periodic boundary conditions, and viscocity $0<ν\ll 1$. An explicit family of quasi-stationary states of the deterministic vorticity equation is known to play an important role in the long-time evolution of solutions both in the presence of and without noise. Recent results show the parameter $δ$ plays a central role in selecting which of the quasi-stationary states is most important. In this paper, we aim to develop a finite dimensional model that captures this selection mechanism for the stochastic vorticity equation. This is done by projecting the vorticity equation in Fourier space onto a center manifold corresponding to the lowest eight Fourier modes. Through Monte Carlo simulation, the vorticity equation and the model are shown to be in agreement regarding key aspects of the long-time dynamics. Following this comparison, perturbation analysis is performed on the model via averaging and homogenization techniques to determine the leading order dynamics for statistics of interest for $δ\approx1$.

preprint2020arXiv

Strong convergence of a Verlet integrator for the semi-linear stochastic wave equation

The full discretization of the semi-linear stochastic wave equation is considered. The discontinuous Galerkin finite element method is used in space and analyzed in a semigroup framework, and an explicit stochastic position Verlet scheme is used for the temporal approximation. We study the stability under a CFL condition and prove optimal strong convergence rates of the fully discrete scheme. Numerical experiments illustrate our theoretical results. Further, we analyze and bound the expected energy and numerically show excellent agreement with the energy of the exact solution.

preprint2016arXiv

A positivity preserving convergent event based asynchronous PDE solver

A new numerical scheme for conservation equations based on evolution by asynchronous discrete events is presented. During each event of the scheme only two cells of the underlying Cartesian grid are active, and an event is processed as the exact evolution of this subsystem. This naturally leads to and adaptive scheme in space and time. Numerical results are presented which show that the error of the asynchronous scheme decreases to zero as a control parameter is reduced. The construction of the scheme allows it to be expressed as repeated multiplications of matrix exponentials on an initial state vector; thus techniques such as the Goldberg series and the Baker Campbell Hausdorff (BCH) formula can be used to explore the theoretical properties of the scheme. We present the framework of a convergence proof in this manner.

preprint2016arXiv

New efficient substepping methods for exponential timestepping

Exponential integrators are time stepping schemes which exactly solve the linear part of a semilinear ODE system. This class of schemes requires the approxima- tion of a matrix exponential in every step, and one successful modern method is the Krylov subspace projection method. We investigate the effect of breaking down a single timestep into arbitrary multiple substeps, recycling the Krylov subspace to minimise costs. For these recyling based schemes we analyse the lo- cal error, investigate them numerically and show they can be applied to a large system with 106 unknowns. We also propose a new second order integrator that is found using the extra information from the substeps to form a corrector to increase the overall order of the scheme. This scheme is seen to compare favourably with other order two integrators.

preprint2014arXiv

Basis properties of the p,q-sine functions

We improve the currently known thresholds for basisness of the family of periodically dilated p,q-sine functions. Our findings rely on a Beurling decomposition of the corresponding change of coordinates in terms of shift operators of infinite multiplicity. We also determine refined bounds on the Riesz constant associated to this family. These results seal mathematical gaps in the existing literature on the subject.

preprint2011arXiv

Approximation properties of the $q$-sine bases

For $q>12/11$ the eigenfunctions of the non-linear eigenvalue problem associated to the one-dimensional $q$-Laplacian are known to form a Riesz basis of $L^2(0,1)$. We examine in this paper the approximation properties of this family of functions and its dual, in order to establish non-orthogonal spectral methods for the $p$-Poisson boundary value problem and its corresponding parabolic time evolution initial value problem. The principal objective of our analysis is the determination of optimal values of $q$ for which the best approximation is achieved for a given $p$ problem.