Related papers: Variable stepsize commutator free Lie group integr…
We develop the equations of motion for full body models that describe the dynamics of rigid bodies, acting under their mutual gravity. The equations are derived using a variational approach where variations are defined on the Lie group of…
We show that existing Runge-Kutta methods for ordinary differential equations (odes) can be modified to solve stochastic differential equations (sdes) with strong solutions provided that appropriate changes are made to the way stepsizes are…
We present initial results on Hessian-free force-gradient integrators for lattice field theories. Integrators of this framework promise to provide substantial performance enhancements, particularly for larger lattice volumes where…
Matrix differential Riccati equation (DRE) typically exhibits transient and steady-state phases, posing challenges for fixed-step time integration methods, which may lack accuracy during transients or oversample in steady regimes. In this…
We give a short and elementary introduction to Lie group methods. A selection of applications of Lie group integrators are discussed. Finally, a family of symplectic integrators on cotangent bundles of Lie groups is presented and the notion…
By using commutator methods, we show uniform resolvent estimates and obtain globally smooth operators for self-adjoint injective homogeneous operators $H$ on graded groups, including Rockland operators, sublaplacians and many others. Left…
Multirate integration is an increasingly relevant tool that enables scientists to simulate multiphysics systems. Existing multirate methods are designed for equations whose fast and slow variables can be linearly separated using additive or…
There exist many Runge-Kutta methods (explicit or implicit), more or less adapted to specific problems. Some of them have interesting properties, such as stability for stiff problems or symplectic capability for problems with energy…
In this paper we derive and analyze the properties of explicit singly diagonal implicit Runge-Kutta (ESDIRK) integration methods. We discuss the principles for construction of Runge-Kutta methods with embedded methods of different order for…
Exponential Runge-Kutta methods are a well-established tool for the numerical integration of parabolic evolution equations. However, these schemes are typically developed under the assumption of homogeneous boundary conditions. In this…
In this paper, we extend the popular integral control technique to systems evolving on Lie groups. More explicitly, we provide an alternative definition of "integral action" for proportional(-derivative)-controlled systems whose…
In this paper we present two new methods based on an implicit Runge-Kutta method Gauss which is of algebraic order fourth and has two stages: the first one has zero dispersion and the second one has zero dispersion and zero dissipation. The…
We endow the group of invertible Fourier integral operators on an open}manifold with the structure of an ILH Lie group. This is done by establishing such structures for the groups of invertible pseudodifferential operators and contact…
We explore two classes of exponential integrators in this letter to design nonlinear Fourier transform (NFT) algorithms with a desired accuracy-complexity trade-off and a convergence order of $4$ on an equispaced grid. The integrating…
The goal of this project is to compare the performance of exponential time integrators with traditional methods such as diagonally implicit Runge-Kutta methods in the context of solving the system of reduced magnetohydrodynamics (RMHD). In…
This study introduces new time-stepping strategies with built-in global error estimators. The new methods propagate the defect along with the numerical solution much like solving for the correction or Zadunaisky's procedure; however, the…
Magnetic quadrupoles are essential components of particle accelerators like the Large Hadron Collider. In order to study numerically the stability of the particle beam crossing a quadrupole, a large number of particle revolutions in the…
Time integration of Fourier pseudo-spectral DNS is usually performed using the classical fourth-order accurate Runge--Kutta method, or other methods of second or third order, with a fixed step size. We investigate the use of higher-order…
We examine various integration schemes for the time-dependent Kohn-Sham equations. Contrary to the time-dependent Schr\"odinger's equation, this set of equations is non-linear, due to the dependence of the Hamiltonian on the electronic…
We propose a practical implementation of high-order fully implicit Runge-Kutta(IRK) methods in a multiple precision floating-point environment. Although implementations based on IRK methods in an IEEE754 double precision environment have…