相关论文: Optimal convergence of a second order low-regulari…
We prove strong convergence of a semi-discrete finite difference method for the KdV and modified KdV equations. We extend existing results to non-smooth data (namely, in $L^2$), without size restrictions. Our approach uses a fourth order…
We analyze a divergence based first order system least squares method applied to a second order elliptic model problem with homogeneous boundary conditions. We prove optimal convergence in the $L^2(\Omega)$ norm for the scalar variable.…
A hybrid LDG-HWENO scheme is proposed for the numerical solution of KdV-type partial differential equations. It evolves the cell averages of the physical solution and its moments (a feature of Hermite WENO) while discretizes high order…
We present a novel class of oscillatory integrators for the Klein-Gordon-Zakharov system which are uniformly accurate with respect to the plasma frequency $c$. Convergence holds from the slowly-varying low-plasma up to the highly…
A linear implicit finite difference method is proposed for the approximation of the solution to a periodic, initial value problem for a Schrodinger-Hirota equation. Optimal, second order convergence in the discrete $H^1-$norm is proved,…
We establish a novel numerical and analytical framework for solving the Korteweg--de Vries (KdV) equation in the negative Sobolev spaces, where classical numerical methods fail due to their reliance on high regularity and inability to…
We analyze a fully discrete finite element numerical scheme for the Cahn-Hilliard-Stokes-Darcy system that models two-phase flows in coupled free flow and porous media. To avoid a well-known difficulty associated with the coupling between…
In this paper, we consider the numerical solution of the continuous disordered nonlinear Schr\"odinger equation, which contains a spatial random potential. We address the finite time accuracy order reduction issue of the usual numerical…
An implicit finite difference scheme based on the $L2$-$1_{\sigma}$ formula is presented for a class of one-dimensional time fractional reaction-diffusion equations with variable coefficients and time drift term. The unconditional stability…
In this paper, we study a second-order accurate and linear numerical scheme for the nonlocal Cahn-Hilliard equation. The scheme is established by combining a modified Crank-Nicolson approximation and the Adams-Bashforth extrapolation for…
We develop efficient numerical integration methods for computing an integral whose integrand is a product of a smooth function and the Gaussian function with a small standard deviation. Traditional numerical integration methods applied to…
An "exact discretization" of the Schroedinger operator is considered and its direct and inverse scattering problems are solved. It is shown that a differential-difference nonlinear evolution equation depending on two arbitrary constants can…
The filtered Lie splitting scheme is an established method for the numerical integration of the periodic nonlinear Schr\"{o}dinger equation at low regularity. Its temporal convergence was recently analyzed in a framework of discrete…
Among the family of fourth-order time integration schemes, the two-stage Gauss--Legendre method, which is an implicit Runge--Kutta method based on collocation, is the only superconvergent. The computational cost of this implicit scheme for…
We study the kernel instrumental variable (KIV) algorithm, a kernel-based two-stage least-squares method for nonparametric instrumental variable regression. We provide a convergence analysis covering both identified and non-identified…
In this paper, we propose and analyse a novel third-order low-regularity trigonometric integrator for the semilinear Klein-Gordon equation with non-smooth solution in the $d$-dimensional space, where $d=1,2,3$. The integrator is constructed…
In this paper we devise and analyze an unconditionally stable, second-order-in-time numerical scheme for the Cahn-Hilliard equation in two and three space dimensions. We prove that our two-step scheme is unconditionally energy stable and…
High order exponential integrators require computing linear combination of exponential like $\varphi$-functions of large matrices $A$ times a vector $v$. Krylov projection methods are the most general and remain an efficient choice for…
Estimation of the precision matrix (or inverse covariance matrix) is of great importance in statistical data analysis and machine learning. However, as the number of parameters scales quadratically with the dimension $p$, computation…
In 2013, Cuturi [Cut13] introduced the Sinkhorn algorithm for matrix scaling as a method to compute solutions to regularized optimal transport problems. In this paper, aiming at a better convergence rate for a high accuracy solution, we…