Related papers: High-order geometric integrators for the variation…
Two recent landmark experiments have performed Gaussian boson sampling (GBS) with a non-programmable linear interferometer and threshold detectors on up to 144 output modes (see Refs.~\onlinecite{zhong_quantum_2020,zhong2021phase}). Here we…
The elucidation of many physical problems in science and engineering is subject to the accurate numerical modelling of complex wave propagation phenomena. Over the last decades, high-order numerical approximation for partial differential…
This study proposes an extension of the high-order compact gas-kinetic scheme (CGKS) to compressible flow simulation in an arbitrary Lagrangian-Eulerian (ALE) formulation in unstructured mesh. The ALE method is achieved by subdividing…
We introduce a nonconforming hybrid finite element method for the two-dimensional vector Laplacian, based on a primal variational principle for which conforming methods are known to be inconsistent. Consistency is ensured using penalty…
We study the convergence of higher order schemes for the Cauchy problem associated to the KdV equation. More precisely, we design a Galerkin type implicit scheme which has higher order accuracy in space and first order accuracy in time. The…
A global solution of the Schr\"odinger equation for explicitly time-dependent Hamiltonians is derived by integrating the non-linear differential equation associated with the time-dependent wave operator. A fast iterative solution method is…
High-order implicit shock tracking (fitting) is a class of high-order, optimization-based numerical methods to approximate solutions of conservation laws with non-smooth features by aligning elements of the computational mesh with…
Computationally efficient numerical methods for high-order approximations of convolution integrals involving weakly singular kernels find many practical applications including those in the development of fast quadrature methods for…
This paper presents a Lie-Trotter splitting for inertial Langevin equations (Geometric Langevin Algorithm) and analyzes its long-time statistical properties. The splitting is defined as a composition of a variational integrator with an…
We present an extension to high-order of a first-order Lagrange-projection like method for the approximation of the Euler equations introduced in Coquel {\it et al.} (Math. Comput., 79 (2010), pp.~1493--1533). The method is based on a…
We consider a simple initial-boundary-value problem for the shallow water equations in one space dimension, and also the analogous problem for a symmetric variant of the system. Assuming smoothness of solutions, we discretize these problems…
A variational formulation for accelerated optimization on normed vector spaces was recently introduced in Wibisono et al., and later generalized to the Riemannian manifold setting in Duruisseaux and Leok. This variational framework was…
We introduce a hybrid method to couple continuous Galerkin finite element methods and high-order finite difference methods in a nonconforming multiblock fashion. The aim is to optimize computational efficiency when complex geometries are…
Most numerical methods for time integration use real-valued time steps. Complex time steps, however, can provide an additional degree of freedom, as we can select the magnitude of the time step in both the real and imaginary directions. We…
We design a novel, exactly energy-conserving implicit non-symplectic integration method for an eight-dimensional Hamiltonian system with four degrees of freedom. In our algorithm, each partial derivative of the Hamiltonian with respect to…
We propose a second-order (Hessian or Hessian-free) based optimization method for variational inference inspired by Gaussian backpropagation, and argue that quasi-Newton optimization can be developed as well. This is accomplished by…
We develop an explicit, second-order, variational time integrator for full body dynamics that preserves the momenta of the continuous dynamics, such as linear and angular momenta, and exhibits near-conservation of total energy over…
A numerical integrator is presented that computes a symmetric or skew-symmetric low-rank approximation to large symmetric or skew-symmetric time-dependent matrices that are either given explicitly or are the unknown solution to a matrix…
Inspired by so-called TVD limiter-based second-order schemes for hyperbolic conservation laws, we develop a second-order accurate numerical method for multi-dimensional aggregation equations. The method allows for simulations to be…
Geometrical design is a crucial and challenging strategy for improving the performance of type-II superconductors, because the proper placement of intended defects in the current path contribute to flux pinning, a reduction in dissipation,…