Related papers: Semi-Closed Form Cubature and Applications to Fina…
In this article we present robust, efficient and accurate fully implicit time-stepping schemes and nonlinear solvers for systems of reaction-diffusion equations. The applications of reaction-diffusion systems is abundant in the literature,…
The computational cost of a Monte Carlo algorithm can only be meaningfully discussed when taking into account the magnitude of the resulting statistical error. Aiming for a fixed error per particle, we study the scaling behavior of the…
Bayesian cubature (BC) is a popular inferential perspective on the cubature of expensive integrands, wherein the integrand is emulated using a stochastic process model. Several approaches have been put forward to encode sequential…
In this thesis, I discuss the use of the Auxiliary Field Diffusion Monte Carlo method to compute the ground state of nuclear Hamiltonians, and I show several applications to interesting problems both in nuclear physics and in nuclear…
Stability and convergence of a time-weighted discrete scheme with nonuniform time steps are established for linear reaction-subdiffusion equations. The Caupto derivative is approximated at an offset point by using linear and quadratic…
In this paper, authors successfully construct a new algorithm for the new higher order scheme of weak approximation of SDEs. The algorithm presented here is based on [1][2]. Although this algorithm shares some features with the algorithm…
In this paper, we consider a numerical method for the multi-term Caputo-Fabrizio time-fractional diffusion equations (with orders $\alpha_i\in(0,1)$, $i=1,2,\cdots,n$). The proposed method employs a fast finite difference scheme to…
A recently developed self-healing diffusion Monte Carlo algorithm [PRB 79, 195117] is extended to the calculation of excited states. The formalism is based on an excited-state fixed-node approximation and the mixed estimator of the…
A method for the multifidelity Monte Carlo (MFMC) estimation of statistical quantities is proposed which is applicable to computational budgets of any size. Based on a sequence of optimization problems each with a globally minimizing…
In this paper, we are interested in the study of a problem with fractional derivatives having boundary conditions of integral types. The problem represents a Caputo type advection-diffusion equation where the fractional order derivative…
A family of semi-intrusive uncertainty propagation (UP) methods for multiscale models is introduced. The methods are semi-intrusive in the sense that inspection of the model is limited up to the level of the single scale systems, and…
In this paper, we study the variable-order (VO) time-fractional diffusion equations. For a VO function $\alpha(t)\in(0,1)$, we develop an exponential-sum-approximation (ESA) technique to approach the VO Caputo fractional derivative. The ESA…
The Auxiliary Field Diffusion Monte Carlo method has been applied to simulate droplets of 7 and 8 neutrons. Results for realistic nucleon-nucleon interactions, which include tensor, spin--orbit and three--body forces, plus a standard…
We propose a semi-discrete finite difference multiscale scheme for a concrete corrosion model consisting of a system of two-scale reaction-diffusion equations coupled with an ode. We prove energy and regularity estimates and use them to get…
We study the application of a quasi-Monte Carlo (QMC) method to a class of semi-linear parabolic reaction-diffusion partial differential equations used to model tumor growth. Mathematical models of tumor growth are largely phenomenological…
This paper concerns the numerical solution of the two-dimensional time-dependent partial integro-differential equation (PIDE) that holds for the values of European-style options under the two-asset Kou jump-diffusion model. A main feature…
Cubature formulas, asymptotically optimal with respect to accuracy, are derived for calculating multidimensional weakly singular integrals. They are used for developing a universal code for calculating capacitances of conductors of…
We have reformulated the quantum Monte Carlo (QMC) technique so that a large part of the calculation scales linearly with the number of atoms. The reformulation is related to a recent alternative proposal for achieving linear-scaling QMC,…
The subdiffusion model that involves a Caputo fractional derivative in time is widely used to describe anomalously slow diffusion processes. In this work we aim at recovering the locations of small conductivity inclusions in the model from…
We present the Continuous Empirical Cubature Method (CECM), a novel algorithm for empirically devising efficient integration rules. The CECM aims to improve existing cubature methods by producing rules that are close to the optimal,…