Related papers: A Second-Order Symplectic Integrator for Guiding-C…
A new form of the Bohr-Sommerfeld quantization is given of the motion of guiding center in strong magnetic field. This is obtained by the effective action for the degree of guiding center which is deduced from the coherent state path…
Many applications in computational physics that use numerical integrators based on splitting and composition can benefit from the development of optimized algorithms and from choosing the best ordering of terms. The cost in programming and…
Asymptotic formula is derived for the behavior of the fundamental solution of the second-order elliptic self-adjoint operator with a piecewise-smooth coefficient in front of the senior derivatives near the discontinuity surface of the…
We apply a heuristic method based on counting points over finite fields to the Poincar\'e center problem. We show that this method gives the correct results for homogeneous non linearities of degree 2 and 3. Also we obtain new evidence for…
In this article, we introduce a novel concept for second-order information of a nonsmooth function inspired by the Goldstein eps-subdifferential. It comprises the coefficients of all existing second-order Taylor expansions in an eps-ball…
We report on the use of the recently-developed Mathematica package \emph{VEST} (Vector Einstein Summation Tools) to automatically derive the guiding center transformation. Our Mathematica code employs a recursive procedure to derive the…
In [Solving second order ordinary differential equations by extending the Prelle-Singer method, J. Phys. A: Math.Gen., 34, 3015-3024 (2001)] we defined a function (we called S) associated to a rational second order ordinary differential…
In this work, we study two first-order primal-dual based algorithms, the Gradient Primal-Dual Algorithm (GPDA) and the Gradient Alternating Direction Method of Multipliers (GADMM), for solving a class of linearly constrained non-convex…
We consider the problem of the estimation of the mean function of an inhomogeneous Poisson process when its intensity function is periodic. For the mean integrated squared error (MISE) there is a classical lower bound for all estimators and…
This work introduces the high-order Boris-SDC method for integrating the equations of motion for electrically charged particles in an electric and magnetic field. Boris-SDC relies on a combination of the Boris-integrator with spectral…
Long-term stability studies of nonlinear Hamiltonian systems require symplectic integration algorithms which are both fast and accurate. In this paper, we study a symplectic integration method wherein the symplectic map representing the…
In this brief, we discuss the implementation of a third order semi-implicit differentiator as a complement of the recent work by the author that proposes an interconnected semi-implicit Euler double differentiators algorithm through Taylor…
An explicit second-order numerical method to integrate the isokinetic equations of motion is derived by fitting circular arcs through every three consecutive points of the discretized trajectory, so that the tangent and the curvature…
In this work, we address the problem of Hessian inversion bias in distributed second-order optimization algorithms. We introduce a novel shrinkage-based estimator for the resolvent of gram matrices which is asymptotically unbiased, and…
An integrator for a class of stochastic Lie-Poisson systems driven by Stratonovich noise is developed. The integrator is suited for Lie-Poisson systems that also admit an isospectral formulation, which enables scalability to…
This paper contributes to the exploration of a recently introduced computational paradigm known as second-order flows, which are characterized by novel dissipative hyperbolic partial differential equations extending accelerated gradient…
The dynamic equation of mass point in rotating coordinates is governed by Coriolis and centrifugal force, besides a corotating potential relative to frame. Such a system is no longer a canonical Hamiltonian system so that the construction…
We develop a second order primal-dual method for optimization problems in which the objective function is given by the sum of a strongly convex twice differentiable term and a possibly nondifferentiable convex regularizer. After introducing…
We revisit the problem of introducing an a priori control for devices that can be modeled via a symplectic map in a neighborhood of an elliptic equilibrium. Using a technique based on Lie transform methods we produce a normal form algorithm…
Recently a new class of numerical integration methods -- ``mixed variable symplectic integrators'' -- has been introduced for studying long-term evolution in the conservative gravitational few-body problem. These integrators are an order of…