Related papers: Correction of high-order $L_k$ approximation for s…
In this work, we present numerical analysis for a distributed optimal control problem, with box constraint on the control, governed by a subdiffusion equation which involves a fractional derivative of order $\alpha\in(0,1)$ in time. The…
In this article, two kinds of numerical algorithms are derived for the ultra-slow (or superslow) diffusion equation in one and two space dimensions, where the ultra-slow diffusion is characterized by the Caputo-Hadamard fractional…
We prove that the superconvergence of $C^0$-$Q^k$ finite element method at the Gauss Lobatto quadrature points still holds if variable coefficients in an elliptic problem are replaced by their piecewise $Q^k$ Lagrange interpolant at the…
We provide the first generic exact simulation algorithm for multivariate diffusions. Current exact sampling algorithms for diffusions require the existence of a transformation which can be used to reduce the sampling problem to the case of…
We present a novel numerical scheme to approximate the solution map $s\mapsto u(s) := \mathcal{L}^{-s}f$ to partial differential equations involving fractional elliptic operators. Reinterpreting $\mathcal{L}^{-s}$ as interpolation operator…
Due to the intrinsically initial singularity of solution and the discrete convolution form in numerical Caputo derivatives, the traditional $H^1$-norm analysis (corresponding to the case for a classical diffusion equation) to the time…
We present a method to rapidly approximate convolution quadrature (CQ) approximations, based on a piecewise polynomial interpolation of the Laplace domain operator, which we call the \emph{parsimonious} convolution quadrature method. For…
The infinitesimal generator (fractional Laplacian) of a process obtained by subordinating a killed Brownian motion catches the power-law attenuation of wave propagation. This paper studies the numerical schemes for the stochastic wave…
Discrete diffusion models have emerged as a powerful generative modeling framework for discrete data with successful applications spanning from text generation to image synthesis. However, their deployment faces challenges due to the high…
In this work, we theoretically and numerically discuss the time fractional subdiffusion-normal transport equation, which depicts a crossover from sub-diffusion (as $t\rightarrow 0$) to normal diffusion (as $t\rightarrow \infty$). Firstly,…
For quasi-linear interface problems with discontinuous diffusion coefficients, the nonconvex objective functional often leads to optimization stagnation in randomized neural network approximations. This paper Proposes a…
We couple the L1 discretization for Caputo derivative in time with spectral Galerkin method in space to devise a scheme that solves quasilinear subdiffusion equations. Both the diffusivity and the source are allowed to be nonlinear…
Time harmonic inverse scattering using accurate forward models is often computationally expensive. On the other hand, the use of computationally efficient solvers, such as the Born approximation, may fail if the targets do not satisfy the…
We consider a class of numerical approximations to the Caputo fractional derivative. Our assumptions permit the use of nonuniform time steps, such as is appropriate for accurately resolving the behavior of a solution whose derivatives are…
We develop a gauge-invariant frozen Gaussian approximation (GIFGA) method for the linear Schr\"odinger equation (LSE) with periodic potentials in the semiclassical regime. The method generalizes the Herman-Kluk propagator for LSE to the…
We develop a general framework for estimating the $L_\infty(\mathbb{T}^d)$ error for the approximation of multivariate periodic functions belonging to specific reproducing kernel Hilbert spaces (RHKS) using approximants that are…
We give a new algorithm for learning mixtures of $k$ Gaussians (with identity covariance in $\mathbb{R}^n$) to TV error $\varepsilon$, with quasi-polynomial ($O(n^{\text{poly\,log}\left(\frac{n+k}{\varepsilon}\right)})$) time and sample…
In this paper, two kinds of high-order compact finite difference schemes for second-order derivative are developed. Then a second-order numerical scheme for Riemann-Liouvile derivative is established based on fractional center difference…
In this paper, we consider the adaptive Eulerian--Lagrangian method (ELM) for linear convection-diffusion problems. Unlike the classical a posteriori error estimations, we estimate the temporal error along the characteristics and derive a…
The main purpose of this paper is to give a solution to a long-standing unsolved problem concerning the pathwise strong approximation of stochastic differential equations with respect to the global error in the $L_{\infty}$-norm. Typically,…