Related papers: Continuous and Discrete Homotopy Operators with Ap…
We develop a monotone, two-scale discretization for a class of integrodifferential operators of order $2s$, $s \in (0,1)$. We apply it to develop numerical schemes, and derive pointwise convergence rates, for linear and obstacle problems…
In this paper, a backward Euler method combined with finite element discretization in spatial direction is discussed for the equations of motion arising in the $2D$ Oldroyd model of viscoelastic fluids of order one with the forcing term…
We introduce a new algorithm to parallelise the computation of persistent homology of 2D alpha complexes. Our algorithm distributes the input point cloud among the cores which then compute a cover based on a rectilinear grid. We show how to…
This paper presents a novel structure-preserving scheme for Euler equations, focusing on the numerical conservation of entropy and kinetic energy. Explicit flux functions engineered to conserve entropy are introduced within the…
In this article, we prove the convergence of a semi-discrete numerical method applied to a general class of nonlocal nonlinear wave equations where the nonlocality is introduced through the convolution operator in space. The most important…
We prove in this work that the well-known lasso problem can be solved exactly without homotopy using novel differential inclusions techniques. Specifically, we show that a selection principle from the theory of differential inclusions…
In this article we propose and investigate a hierarchy of mathematical models based on partial differential equations (PDE) and ordinary differential equations (ODE) for the simulation of the biophysical phenomena occurring in the…
Noether's theorem connects symmetries to invariants in continuous systems, however its extension to discrete systems has remained elusive. Recognizing the lowest-order finite difference as the foundation of local continuity, a viable method…
Entropy stable schemes can be constructed with a specific choice of the numerical flux function. First, an entropy conserving flux is constructed. Secondly, an entropy stable dissipation term is added to this flux to guarantee dissipation…
A discrete rotation algorithm can be apprehended as a parametric application $f\_\alpha$ from $\ZZ[i]$ to $\ZZ[i]$, whose resulting permutation ``looks like'' the map induced by an Euclidean rotation. For this kind of algorithm, to be…
A novel implicit integration scheme for the Discrete Element Method (DEM) based on the variational integrator approach is presented. The numerical solver provides a fully dynamical description that, notably, reduces to an energy…
Integration of nonlinear partial differential equations with the help of the non-commutative integration over octonions is studied. An apparatus permitting to take into account symmetry properties of PDOs is developed. For this purpose…
Many simulated complex systems that support persistent self-organizing patterns, i.e. gliders, have a 'state-plus-update' paradigm. This approach can be found in computational models of physics, continuous and neural cellular automata,…
Algorithms for the computation of the real zeros of hypergeometric functions which are solutions of second order ODEs are described. The algorithms are based on global fixed point iterations which apply to families of functions satisfying…
Numerical methods that preserve geometric invariants of the system, such as energy, momentum or the symplectic form, are called geometric integrators. Variational integrators are an important class of geometric integrators. The general idea…
In this work we have presented a rather general and easy-to-apply method for discrete Hilbert space representation of quantum mechanical Green's operators. We have shown that if in some discrete Hilbert space basis representation the…
The rapid advancements in high-dimensional statistics and machine learning have increased the use of first-order methods. Many of these methods can be regarded as instances of the proximal point algorithm. Given the importance of the…
We compute fundamental solutions of homogeneous elliptic differential operators, with constant coefficients, on $\mathbb{R}^n$ by mean of analytic continuation of distributions. The result obtained is valid in any dimension, for any degree…
This work aims to construct an efficient and highly accurate numerical method to address the time singularity at $t=0$ involved in a class of time-fractional parabolic integro-partial differential equations in one and two dimensions. The…
We present a new class of exponential integrators for ordinary differential equations: locally exact modifications of known numerical schemes. Local exactness means that they preserve the linearization of the original system at every point.…