Related papers: A systematic method to enforce conservativity on s…
Dimensionally split advection schemes are attractive for atmospheric modelling due to their efficiency and accuracy in each spatial dimension. Accurate long time-steps can be achieved without significant cost using the flux-form…
Clustering may be the most fundamental problem in unsupervised learning which is still active in machine learning research because its importance in many applications. Popular methods like K-means, may suffer from instability as they are…
This paper develops maintenance policies for a system under condition monitoring. We assume that a number of defects may develop and the degradation process of each defect follows a gamma process, respectively. The system is inspected…
Semi-implicit semi-Lagrangian (SISL) methods are commonly used for the shallow water equations (SWE) because they allow for larger time steps than those permitted by the Courant-Friedrichs-Lewy (CFL) stability condition in Eulerian schemes.…
In this paper, we consider a class of stochastic midpoint and trapezoidal Lawson schemes for the numerical discretization of highly oscillatory stochastic differential equations. These Lawson schemes incorporate both the linear drift and…
The linear Support Vector Machine (SVM) is a classic classification technique in machine learning. Motivated by applications in modern high dimensional statistics, we consider penalized SVM problems involving the minimization of a…
A numerical method is proposed in order to track field lines of three-dimensional divergence free fields. Field lines are computed by a locally valid Hamiltonian mapping, which is computed using a symplectic scheme. The method is…
Safety is essential for reinforcement learning (RL) applied in the real world. Adding chance constraints (or probabilistic constraints) is a suitable way to enhance RL safety under uncertainty. Existing chance-constrained RL methods like…
Mass-conservative reaction-diffusion systems have recently been proposed as a general framework to describe intracellular pattern formation. These systems have been used to model the conformational switching of proteins as they cycle from…
An explicit high order semi-Lagrangian method is developed for solving Lagrangian transport equations in Eulerian-Lagrangian formulations. To ensure a semi-Lagrangian approximation that is consistent with an explicit Eulerian, discontinuous…
In this paper, we consider a model reduction technique for stabilizable and detectable stochastic systems. It is based on a pair of Gramians that we analyze in terms of well-posedness. Subsequently, dominant subspaces of the stochastic…
Subgradient methods are the natural extension to the non-smooth case of the classical gradient descent for regular convex optimization problems. However, in general, they are characterized by slow convergence rates, and they require…
In this exploratory study, we apply shock-capturing schemes within the framework of the Particles on Demand kinetic model to simulate compressible flows with mild and strong shock waves and discontinuities. The model is based on the…
There are several well-established approaches to constructing finite difference schemes that preserve global invariants of a given partial differential equation. However, few of these methods preserve more than one conservation law locally.…
We present a method for constructing numerical schemes with up to 3rd strong convergence order for solution of a class of stochastic differential equations, including equations of the Langevin type. The construction proceeds in two stages.…
We propose a novel, highly efficient, mean-reverting-SAV-BDF2-based, long-time unconditionally stable numerical scheme for a class of finite-dimensional nonlinear models important in geophysical fluid dynamics. The scheme is highly…
The variational free-Lagrange (VFL) method for shallow water is a free-Lagrange method with the additional property that it preserves the variational structure of shallow water. The VFL method was first derived in this context by…
Within recent years, considerable progress has been made regarding high-performance solvers for Partial Differential Equations (PDEs), yielding potential gains in efficiency compared to industry standard tools. However, the latter largely…
Dynamical low-rank approximation, as has been demonstrated recently, can be extremely efficient in solving kinetic equations. However, a major deficiency is that they do not preserve the structure of the underlying physical problem. For…
The aim of this paper is to develop some tools in order to obtain the weak consistency of (in other words, analogues of the Lax-Wendroff theorem for) finite volume schemes for balance laws in the multi-dimensional case and under minimal…