Related papers: High-order geometric integrators for the variation…
Hamiltonian Monte Carlo (HMC) is a popular method in sampling. While there are quite a few works of studying this method on various aspects, an interesting question is how to choose its integration time to achieve acceleration. In this…
A new conservative symmetry-preserving second-order time-accurate PISO-based pressure-velocity coupling for solving the incompressible Navier-Stokes equations on unstructured collocated grids is presented in this paper. This new method for…
We propose an efficient approach for time integration of Klein-Gordon equations with highly oscillatory in time input terms. The new methods are highly accurate in the entire range, from slowly varying up to highly oscillatory regimes. Our…
We discuss the efficient implementation of a high-performance second-order collocation-type finite-element scheme for solving the compressible Euler equations of gas dynamics on unstructured meshes. The solver is based on the convex…
In this work, we propose an adaptive spectral element algorithm for solving nonlinear optimal control problems. The method employs orthogonal collocation at the shifted Gegenbauer-Gauss points combined with very accurate and stable…
A novel class of high-order linearly implicit energy-preserving integrating factor Runge-Kutta methods are proposed for the nonlinear Schr\"odinger equation. Based on the idea of the scalar auxiliary variable approach, the original equation…
A high-order well-balanced scheme for the Euler equations with gravitation is presented. The scheme is able to preserve a spatially high-order accurate discrete representation of a large class of hydrostatic equilibria. It is based on a…
Vecchia's approximate likelihood for Gaussian process parameters depends on how the observations are ordered, which can be viewed as a deficiency because the exact likelihood is permutation-invariant. This article takes the alternative…
In this work, a second-order approximation of the fractional substantial derivative is presented by considering a modified shifted substantial Gr\"{u}nwald formula and its asymptotic expansion. Moreover, the proposed approximation is…
We develop sampling methods, which consist of Gaussian invariant versions of random walk Metropolis (RWM), Metropolis adjusted Langevin algorithm (MALA) and second order Hessian or Manifold MALA. Unlike standard RWM and MALA we show that…
A new second-order method for approximating the compressible Euler equations is introduced. The method preserves all the known invariant domains of the Euler system: positivity of the density, positivity of the internal energy and the local…
Elegant integration schemes of second and fourth order for simulations of rigid body systems are presented which treat translational and rotational motion on the same footing. This is made possible by a recent implementation of the exact…
Symplectic integration of autonomous Hamiltonian systems is a well-known field of study in geometric numerical integration, but for non-autonomous systems the situation is less clear, since symplectic structure requires an even number of…
In this paper we construct higher-order variational integrators for a class of degenerate systems described by Lagrangians that are linear in velocities. We analyze the geometry underlying such systems and develop the appropriate theory for…
With increasing engineering demands, there need high order accurate schemes embedded with precise physical information in order to capture delicate small scale structures and strong waves with correct "physics". There are two families of…
We provide an exact analysis of a class of randomized algorithms for solving overdetermined least-squares problems. We consider first-order methods, where the gradients are pre-conditioned by an approximation of the Hessian, based on a…
Following on our previous work [S. Delong and B. E. Griffith and E. Vanden-Eijnden and A. Donev, Phys. Rev. E, 87(3):033302, 2013], we develop temporal integrators for solving Langevin stochastic differential equations that arise in…
We study the order of convergence of Galerkin variational integrators for ordinary differential equations. Galerkin variational integrators approximate a variational (Lagrangian) problem by restricting the space of curves to the set of…
A high order discontinuous Galerkin method for the material transport of thermodynamic tracers is coupled to a low order mixed finite element solver in the context of the thermal shallow water equations. The coupling preserves the energy…
We propose a high-order numerical methodology for computing the ground state and time evolution of the two-dimensional Gross-Pitaevskii equation with harmonic trapping potential. The ground state is obtained by combining normalized gradient…