Related papers: Numerical computation of the half Laplacian by mea…
Compared to the the classical first-order Gr\"unwald-Letnikov formula at time $t_{k+1} (\textmd{or}\, t_{k})$, we firstly propose a second-order numerical approximate scheme for discretizing the Riemann-Liouvile derivative at time…
Randomized approximation algorithms for many #P-complete problems (such as the partition function of a Gibbs distribution, the volume of a convex body, the permanent of a $\{0,1\}$-matrix, and many others) reduce to creating random…
Integral transforms are invaluable mathematical tools to map functions into spaces where they are easier to characterize. We introduce the hyperdimensional transform as a new kind of integral transform. It converts square-integrable…
We develop a novel a posteriori error estimator for the $L^2$ error committed by the finite element discretization of the solution of the fractional Laplacian. Our a posteriori error estimator takes advantage of the semi-discretization…
We present three equivalent definitions of the fractional $p$-Laplacian $(-\Delta_{\mathbb{H}^{n}})^{s}_{p}$, $0<s<1$, $p>1$, with normalizing constants, on hyperbolic spaces. The explicit values of the constants enable us to study the…
We present an exponentially convergent numerical method to approximate the solution of the Cauchy problem for the inhomogeneous fractional differential equation with an unbounded operator coefficient and Caputo fractional derivative in…
In this paper, we are concerned with the following equation involving higher-order fractional Lapalacian \begin{equation*} \left\{\begin{aligned} &(-\Delta)^{p+{\frac{\alpha}{2}}}u(x)=u_+^\gamma~~ \mbox{ in }\mathbb{R}^n,\\…
We present a fast and numerically accurate method for expanding digitized $L \times L$ images representing functions on $[-1,1]^2$ supported on the disk $\{x \in \mathbb{R}^2 : |x|<1\}$ in the harmonics (Dirichlet Laplacian eigenfunctions)…
In this paper we treat the numerical approximation of the two-phase parabolic obstacle-like problem: \[\Delta u -u_t=\lambda^+\cdot\chi_{\{u>0\}}-\lambda^-\cdot\chi_{\{u<0\}},\quad (t,x)\in (0,T)\times\Omega,\] where $T < \infty, \lambda^+…
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…
In this paper we are interested in the numerical approximation of the marginal distributions of the Hilbert space valued solution of a stochastic Volterra equation driven by an additive Gaussian noise. This equation can be written in the…
The Mittag-Leffler function is computed via a quadrature approximation of a contour integral representation. We compare results for parabolic and hyperbolic contours, and give special attention to evaluation on the real line. The main point…
Product-to-sum identities for trigonometric functions play a fundamental role in function theory and numerous applications. In this spirit, we present convolution-to-sum identities for Mittag-Leffler type functions. Using a Laplace domain…
We study the symmetry properties for solutions of elliptic systems of the type (-\Delta)^{s_1} u = F_1(u, v), (-\Delta)^{s_2} v= F_2(u, v), where $F\in C^{1,1}_{loc}(\R^2)$, $s_1,s_2\in (0,1)$ and the operator $(-\Delta)^s$ is the so-called…
In this paper we present a variant of the proximal forward-backward splitting iteration for solving nonsmooth optimization problems in Hilbert spaces, when the objective function is the sum of two nondifferentiable convex functions. The…
We develop an $L^p(\mathbb{R}^n)$-functional calculus appropriated for interpreting "non-classical symbols" of the form $a(-\Delta)$, and for proving existence in $L^q(\mathbb{R}^n)$, some $q > p$, of solutions to nonlinear…
Let $f$ be a fixed (holomorphic or Maass) modular cusp form. Let $\cq$ be a Dirichlet character mod $q$. We describe a fast algorithm that computes the value $L(1/2,f\times\chi_q)$ up to any specified precision. In the case when $q$ is…
We introduce a simple method for nearly simultaneous computation of all moments needed for quasi maximum likelihood estimation of parameters in discretely observed stochastic differential equations commonly seen in finance. The method…
While exploiting the generalized Parseval equality for the Mellin transform, we derive the reciprocal inverse operator in the weighted L_2-space related to the Hilbert transform on the nonnegative half-axis. Moreover, employing the…
We study the extremal solution for the problem $(-\Delta)^s u=\lambda f(u)$ in $\Omega$, $u\equiv0$ in $\R^n\setminus\Omega$, where $\lambda>0$ is a parameter and $s\in(0,1)$. We extend some well known results for the extremal solution when…