Related papers: Efficient Numerical Method for Models Driven by L\…
Boundary integral equations lead to dense system matrices when discretized, yet they are data-sparse. Using the $\mathcal{H}$-matrix format, this sparsity is exploited to achieve $\mathcal{O}(N\log N)$ complexity for storage and…
In this article, we are concerned with the analysis on the numerical reconstruction of the spatial component in the source term of a time-fractional diffusion equation. This ill-posed problem is solved through a stabilized nonlinear…
The iterated Crank-Nicolson is a predictor-corrector algorithm commonly used in numerical relativity for the solution of both hyperbolic and parabolic partial differential equations. We here extend the recent work on the stability of this…
This article presents a new finite element method for convection-diffusion equations by enhancing the continuous finite element space with a flux space for flux approximations that preserve the important mass conservation locally on each…
A simple yet effective numerical method using orthogonal hybrid functions consisting of piecewise constant orthogonal sample-and-hold functions and piecewise linear orthogonal triangular functions is proposed to solve numerically fractional…
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…
In this work, two novel classes of structure-preserving spectral Galerkin methods are proposed which based on the Crank-Nicolson scheme and the exponential scalar auxiliary variable method respectively, for solving the coupled fractional…
The Fokker-Planck equations for stochastic dynamical systems, with non-Gaussian $\alpha-$stable symmetric L\'evy motions, have a nonlocal or fractional Laplacian term. This nonlocality is the manifestation of the effect of non-Gaussian…
Standard discretization techniques for boundary integral equations, e.g., the Galerkin boundary element method, lead to large densely populated matrices that require fast and efficient compression techniques like the fast multipole method…
In this work, we develop an $\mathcal{O}(N)$ implicit real space method in 1D and 2D for the Cahn Hilliard (CH) and vector Cahn Hilliard (VCH) equations, based on the Method Of Lines Transpose (MOL$^\text{T}$) formulation. This formulation…
This paper is about operator-theoretic methods for solving nonlinear stochastic optimal control problems to global optimality. These methods leverage on the convex duality between optimally controlled diffusion processes and…
In this paper, we introduce and analyse numerical schemes for the homogeneous and the kinetic L\'evy-Fokker-Planck equation. The discretizations are designed to preserve the main features of the continuous model such as conservation of…
The breakthrough ideas in the modern proximal splitting methodologies allow us to express the set of all minimizers of a superposition of multiple nonsmooth convex functions as the fixed point set of computable nonexpansive operators. In…
It is well understood that, when numerically simulating SDEs with general noise, achieving a strong convergence rate better than $O(\sqrt{h})$ (where h is the step size) requires the use of certain iterated integrals of Brownian motion,…
In this paper we extend the hierarchical model reduction framework based on reduced basis techniques for the application to nonlinear partial differential equations. The major new ingredient to accomplish this goal is the introduction of…
A class of Hamiltonian stochastic differential equations with multiplicative L\'{e}vy noise in the sense of Marcus, and the construction and numerical implementation methods of symplectic Euler scheme, are considered. A general symplectic…
We study nonlinear elliptic equations for operators corresponding to non-stable L\'evy diffusions. We include a sum of fractional Laplacians of different orders. Such operators are infinitesimal generators of non-stable (i.e., non…
We give a probabilistic numerical method for solving a partial differential equation with fractional diffusion and nonlinear drift. The probabilistic interpretation of this equation uses a system of particles driven by L\'evy alpha-stable…
An efficient Krylov subspace algorithm for computing actions of the $\varphi$ matrix function for large matrices is proposed. This matrix function is widely used in exponential time integration, Markov chains and network analysis and many…
Stochastic optimal control problems have a long tradition in applied probability, with the questions addressed being of high relevance in a multitude of fields. Even though theoretical solutions are well understood in many scenarios, their…