Related papers: Algebraic Structures and Stochastic Differential E…
We introduce a general framework of low regularity integrators which allows us to approximate the time dynamics of a large class of equations, including parabolic and hyperbolic problems, as well as dispersive equations, up to arbitrary…
We discuss stochastic differential equations with a stiff linear part and their approximation by stochastic exponential integrators. Representing the exact and approximate solutions using B-series and rooted trees, we derive the order…
In this paper, we deal with a class of reflected backward stochastic differential equations associated to the subdifferential operator of a lower semi-continuous convex function driven by Teugels martingales associated with L\'{e}vy…
For sampling from a log-concave density, we study implicit integrators resulting from $\theta$-method discretization of the overdamped Langevin diffusion stochastic differential equation. Theoretical and algorithmic properties of the…
In this work, the benefits of the phase fitting technique are embedded in high order discrete Lagrangian integrators. The proposed methodology creates integrators with zero phase lag in a test Lagrangian in a similar way used in phase…
This paper investigates the parareal algorithms for solving the stochastic Maxwell equations driven by multiplicative noise, focusing on their convergence, computational efficiency and numerical performance. The algorithms use the…
The purpose of this paper is to adapt the empirical characteristic function (ECF) method to stable, but possibly not inverse stable linear stochastic system driven by the increments of a Levy-process. A remarkable property of the ECF method…
The efficient numerical integration of large-scale matrix differential equations is a topical problem in numerical analysis and of great importance in many applications. Standard numerical methods applied to such problems require an unduly…
Several numerical schemes are proposed for the solution of Nonequilibrium Langevin Dynamics (NELD), and the rate of convergence is analyzed. Due to the special deforming boundary conditions used, care must be taken when using standard…
We define the notion of effective stiffness and show that it can used to build sparsifiers, algorithms that sparsify linear systems arising from finite-element discretizations of PDEs. In particular, we show that sampling $O(n\log n)$…
We propose an efficient approach for time integration of Klein-Gordon equations with highly oscillatory in time input terms. The new methods are highly accurate in the entire range, from slowly varying up to highly oscillatory regimes. Our…
The dynamical low-rank approximation of time-dependent matrices is a low-rank factorization updating technique. It leads to differential equations for factors of the matrices, which need to be solved numerically. We propose and analyze a…
To reliably model real robot characteristics, interval linear systems of equations allow to describe families of problems that consider sets of values. This allows to easily account for typical complexities such as sets of joint states and…
We continue the investigation of the Levy processes on a q-deformed full Fock space started in a previous paper. First, we show that the vacuum vector is cyclic and separating for the algebra generated by such a process. Next, we describe a…
Derivatives play a critical role in computational statistics, examples being Bayesian inference using Hamiltonian Monte Carlo sampling and the training of neural networks. Automatic differentiation is a powerful tool to automate the…
In this paper we consider an approach to improve the performance of exponential Runge--Kutta integrators and Lawson schemes} in cases where the solution of a related, but usually much simpler, problem can be computed efficiently. While for…
We explore the connection between fractional order partial differential equations in two or more spatial dimensions with boundary integral operators to develop techniques that enable one to efficiently tackle the integral fractional…
We consider the problem of the simulation of Levy-driven stochastic differential equations. It is generally impossible to simulate the increments of a Levy-process. Thus in addition to an Euler scheme, we have to simulate approximately…
This paper is a summary of the theory of discrete embeddings introduced in [5]. A discrete embedding is an algebraic procedure associating a numerical scheme to a given ordinary differential equation. Lagrangian systems possess a…
In this work, we provide the first strong convergence result of numerical approximation of a general second order semilinear stochastic fractional order evolution equation involving a Caputo derivative in time of order $\alpha\in(\frac 34,…