Related papers: Modified Hermite Integrators of Arbitrary Order
In this paper we propose a fast optimization algorithm for approximately minimizing convex quadratic functions over the intersection of affine and separable constraints (i.e., the Cartesian product of possibly nonconvex real sets). This…
Numerical algorithms for elliptic partial differential equations frequently employ error estimators and adaptive mesh refinement strategies in order to reduce the computational cost. We can extend these techniques to general vectors by…
I show that the basic structure of symplectic integrators is governed by a theorem which states {\it precisely}, how symplectic integrators with positive coefficients cannot be corrected beyond second order. All previous known results can…
The higher-order superintegrability of separable potentials is studied. It is proved that these potentials possess (in addition to the two quadratic integrals) a third integral of higher-order in the momenta that can be obtained as the…
We show that a recently discovered fourth order symplectic algorithm, which requires one evaluation of force gradient in addition to three evaluations of the force, when iterated to higher order, yielded algorithms that are far superior to…
This paper describes a regularized variant of the alternating direction method of multipliers (ADMM) for solving linearly constrained convex programs. It is shown that the pointwise iteration-complexity of the new method is better than the…
In this paper, a piecewise quadratic nonconforming finite element method on rectangular grids for a fourth-order elliptic singular perturbation problem is presented. This proposed method is robustly convergent with respect to the…
In this article we propose a new adaptive numerical quadrature procedure which includes both local subdivision of the integration domain, as well as local variation of the number of quadrature points employed on each subinterval. In this…
It is well known that in perturbation theory existence of quasi-degenerate states can rearrange the order counting. For a heavy quarkonium system, naively, enhanced effects ($l$-changing mixing effects) could contribute already to the…
In this paper, a meshless Hermite-HDMR finite difference method is proposed to solve high-dimensional Dirichlet problems. The approach is based on the local Hermite-HDMR expansion with an additional smoothing technique. First, we introduce…
Gaussian wavepacket dynamics has proven to be a useful semiclassical approximation for quantum simulations of high-dimensional systems with low anharmonicity. Compared to Heller's original local harmonic method, the variational Gaussian…
The evolution of any factorized time-reversible symplectic integrators, when applied to the harmonic oscillator, can be exactly solved in a closed form. The resulting modified Hamiltonians demonstrate the convergence of the Lie series…
Corrector estimates constitute a key ingredient in the derivation of optimal convergence rates via two-scale expansion techniques in homogenization theory of random uniformly elliptic equations. The present work follows up - in terms of…
We consider adaptive finite element methods for second-order elliptic PDEs, where the arising discrete systems are not solved exactly. For contractive iterative solvers, we formulate an adaptive algorithm which monitors and steers the…
Higher order interactions can lead to new equilibrium states and bifurcations in systems of coupled oscillators described by the Kuramoto model. However, even in the simplest case of 3-body interactions there are more than one possible…
Multiplication is an indispensable operation in most of digital signal processing systems. Recently, many systems need to execute different types of algorithms on a multiplier. Therefore, it needs complicated computation and large area…
In the last two decades, significant effort has been put in understanding and designing so-called structure-preserving numerical methods for the simulation of mechanical systems. Geometric integrators attempt to preserve the geometry…
In this work we construct an Hermite interpolant starting from basis functions that satisfy a Lagrange property. In fact, we extend and generalise an iterative approach, introduced by Cirillo and Hormann (2018) for the Floater-Hormann…
Symplectic integrators separate a problem into parts that can be solved in isolation, alternately advancing these sub-problems to approximate the evolution of the complete system. Problems with a single, dominant mass can use mixed-variable…
Exponential integrators are special time discretization methods where the traditional linear system solves used by implicit schemes are replaced with computing the action of matrix exponential-like functions on a vector. A very general…