Related papers: Efficient discretization of the Laplacian on compl…
This paper extends previous work on finitedifference schemes over staggered grids for infinite-dimensional port-Hamiltonian systems. In the one-dimensional setting, it generalizes the discretization approach originally developed for the…
In this paper, we discuss the stable discretisation of the double layer boundary integral operator for the wave equation in $1d$. For this, we show that the boundary integral formulation is $L^2$-elliptic and also inf-sup stable in standard…
We consider a vector-Laplace problem posed on a 2D surface embedded in a 3D domain, which results from the modeling of surface fluids based on exterior Cartesian differential operators. The main topic of this paper is the development and…
In this work is considered a spectral problem, involving a second order term on the domain boundary: the Laplace-Beltrami operator. A variational formulation is presented, leading to a finite element discretization. For the Laplace-Beltrami…
We show that discrete schemes developed for lattice hydrodynamics provide an elegant and physically transparent way of deriving Laplacians with isotropic discretisation error. Isotropy is guaranteed whenever the Laplacian weights follow…
Many computational algorithms applied to geometry operate on discrete representations of shape. It is sometimes necessary to first simplify, or coarsen, representations found in modern datasets for practicable or expedited processing. The…
We construct new, efficient, and accurate high-order finite differencing operators which satisfy summation by parts. Since these operators are not uniquely defined, we consider several optimization criteria: minimizing the bandwidth, the…
This work aims to construct an efficient and highly accurate numerical method to address the time singularity at $t=0$ involved in a class of time-fractional parabolic integro-partial differential equations in one and two dimensions. The…
The Stokes system with constant viscosity can be cast into different formulations by exploiting the incompressibility constraint. For instance the strain in the weak formulation can be replaced by the gradient to decouple the velocity…
We propose highly accurate finite-difference schemes for simulating wave propagation problems described by linear second-order hyperbolic equations. The schemes are based on the summation by parts (SBP) approach modified for applications…
This paper is concerned with the construction of discrete counterparts of the Laplace-Beltrami operator on Riemannian manifolds that can be effectively used in the numerical solution of partial differential equations. Since existing…
We build convergent discretizations and semi-implicit solvers for the Infinity Laplacian and the game theoretical $p$-Laplacian. The discretizations simplify and generalize earlier ones. We prove convergence of the solution of the Wide…
In this work we explore the fidelity of numerical approximations to the analytic spectra of hyperbolic partial differential equation systems with variable coefficients. We are particularly interested in the ability of discrete methods to…
The scalar, one-dimensional advection equation and heat equation are considered. These equations are discretized in space, using a finite difference method satisfying summation-by-parts (SBP) properties. To impose the boundary conditions,…
The data input model is a fundamental component of every quantum algorithm, as its efficiency is crucial for achieving potential speed-ups over classical methods. For quantum linear algebra tasks that utilize quantum eigenvalue or singular…
We consider approximations to the solutions of differential Riccati equations in the context of linear quadratic regulator problems, where the state equation is governed by a multiscale operator. Similarly to elliptic and parabolic…
We propose a novel and simple spectral method based on the semi-discrete Fourier transforms to discretize the fractional Laplacian $(-\Delta)^\frac{\alpha}{2}$. Numerical analysis and experiments are provided to study its performance. Our…
We investigate the use of renormalisation group methods to solve partial differential equations (PDEs) numerically. Our approach focuses on coarse-graining the underlying continuum process as opposed to the conventional numerical analysis…
We employ the summation-by-parts (SBP) framework to extend the recent domain-of-dependence (DoD) stabilization for cut cells to linear kinetic models in diffusion scaling. Numerical methods for these models are challenged by increased…
In this work we extend the shifted Laplacian approach to the elastic Helmholtz equation. The shifted Laplacian multigrid method is a common preconditioning approach for the discretized acoustic Helmholtz equation. In some cases, like…