Related papers: Efficient Numerical Method for Models Driven by L\…
In this paper, we develop an efficient numerical solver for unsteady diffusion-type partial differential equations with random coefficients. A major computational challenge in such problems lies in repeatedly handling large-scale linear…
In this work, we study the accuracy and efficiency of hierarchical matrix ($\mathcal{H}$-matrix) based fast methods for solving dense linear systems arising from the discretization of the 3D elastodynamic Green's tensors. It is well known…
This paper detailedly discusses the locally one-dimensional numerical methods for efficiently solving the three-dimensional fractional partial differential equations, including fractional advection diffusion equation and Riesz fractional…
The fully discrete problem for convection-diffusion equation is considered. It comprises compact approximations for spatial discretization, and Crank-Nicolson scheme for temporal discretization. The expressions for the entries of inverse of…
We prove that weakly continuous solutions to martingale problems admit a canonical regular conditional probability distribution. This allows for the construction of time consistent convex dynamic procedures in a non dominated setting.…
Fractional Fokker-Planck equation plays an important role in describing anomalous dynamics. To the best of our knowledge, the existing discussions mainly focus on this kind of equation involving one diffusion operator. In this paper, we…
This paper presents a new method to approximate the time-dependent convection-diffusion equations using conforming finite element methods, ensuring that the discrete solution respects the physical bounds imposed by the differential…
We are interested in the numerical solution of nonsymmetric linear systems arising from the discretization of convection-diffusion partial differential equations with separable coefficients and dominant convection. Preconditioners based on…
This paper proposes a replica exchange preconditioned Langevin diffusion discretized by the Crank-Nicolson scheme (repCNLD) to handle high-dimensional and multi-modal distribution problems. Sampling from high-dimensional and multi-modal…
In this work, we analyze a Crank-Nicolson type time stepping scheme for the subdiffusion equation, which involves a Caputo fractional derivative of order $\alpha\in (0,1)$ in time. It hybridizes the backward Euler convolution quadrature…
In this paper, we develop and analyze a stochastic algorithm for solving space-time fractional diffusion models, which are widely used to describe anomalous diffusion dynamics. These models pose substantial numerical challenges due to the…
Several researchers have developed a rich toolbox of matrix compression techniques that exploit structure and redundancy in large matrices. Classical methods such as the block low-rank format and the Fast Multipole Method make it possible…
We develop a unified and easy to use framework to study robust fully discrete numerical methods for nonlinear degenerate diffusion equations $$ \partial_t u-\mathfrak{L}^{\sigma,\mu}[\varphi(u)]=f \quad\quad\text{in}\quad\quad…
In this paper, we study the Crank-Nicolson method for temporal dimension and the piecewise quadratic polynomial collocation method for spatial dimensions of time-dependent nonlocal problems. The new theoretical results of such…
A second-order accurate in time, positivity-preserving, and unconditionally energy stable operator splitting numerical scheme is proposed and analyzed for the system of reaction-diffusion equations with detailed balance. The scheme is…
We describe an exact and highly efficient numerical algorithm for solving a special but important class of convection-diffusion equations. These equations occur in many problems in physics, chemistry, or biology, and they are usually hard…
Along with the recent advances in scalable Markov Chain Monte Carlo methods, sampling techniques that are based on Langevin diffusions have started receiving increasing attention. These so called Langevin Monte Carlo (LMC) methods are based…
We propose and study discontinuous Galerkin methods for strongly degenerate convection-diffusion equations perturbed by a fractional diffusion (L\'evy) operator. We prove various stability estimates along with convergence results toward…
Achieving a substantial part of peak performance on todays and future high-performance computing systems is a major challenge for simulation codes. In this paper we address this question in the context of the numerical solution of partial…
In this paper, we develop a new mathematical technique which allows us to express the joint distribution of a Markov process and its running maximum (or minimum) through the marginal distribution of the process itself. This technique is an…