Related papers: Solving Kepler's equation with CORDIC double itera…
We consider an overdetermined problem of Serrin-type with respect to an operator in divergence form with piecewise constant coefficients. We give sufficient condition for unique solvability near radially symmetric configurations by means of…
Numerically obtaining the inverse of a function is a common task for many scientific problems, often solved using a Newton iteration method. Here we describe an alternative scheme, based on switching variables followed by spline…
In this paper we introduce the algorithm and the fixed point hardware to calculate the normalized singular value decomposition of a non-symmetric matrices using Givens fast (approximate) rotations. This algorithm only uses the basic…
Compact binaries consisting of neutron stars / black holes on eccentric orbit undergo a perturbed Keplerian motion. The perturbations are either of relativistic origin or are related to the spin, mass quadrupole and magnetic dipole moments…
We give a self-contained presentation and comparison of two different algorithms to explicitly solve quantum many body models of indistinguishable particles moving on a circle and interacting with two-body potentials of $1/\sin^2$-type. The…
We provide an elegant way of solving analytically the third post-Newtonian (3PN) accurate Kepler equation, associated with the 3PN-accurate generalized quasi-Keplerian parametrization for compact binaries in eccentric orbits. An additional…
We investigate the Kepler problem using a symplectic structure consistent with the commutation rules of the noncommutative quantum mechanics. We show that a noncommutative parameter of the order of $10^{-58} \text m^2$ gives observable…
This paper studies multistep methods for the integration of reversible dynamical systems, with particular emphasis on the planar Kepler problem. It has previously been shown by Cano & Sanz-Serna that reversible linear multisteps for…
We show that when time-reversible symplectic algorithms are used to solve periodic motions, the energy error after one period is generally two orders higher than that of the algorithm. By use of correctable algorithms, we show that the…
We consider a single particle which is bound by a central potential and obeys the Dirac equation in d dimensions. We first apply the asymptotic iteration method to recover the known exact solutions for the pure Coulomb case. For a…
A new, second-order solution in curvilinear coordinates is introduced for the relative motion of two spacecraft on eccentric orbits. The second-order equations for unperturbed orbits are derived in spherical coordinates with true anomaly as…
Finding roots of equations is at the heart of most computational science. A well-known and widely used iterative algorithm is the Newton's method. However, its convergence depends heavily on the initial guess, with poor choices often…
An unconventional approach is applied to solve the one-dimensional Burgers' equation. It is based on spline polynomial interpolations and Hopf-Cole transformation. Taylor expansion is used to approximate the exponential term in the…
We investigate a method to compute a finite set of preliminary orbits for solar system bodies using the first integrals of the Kepler problem. This method is thought for the applications to the modern sets of astrometric observations, where…
The sum of elliptic integrals simultaneously determines orbits in thr Kepler problem and the addition of divisors on elliptic curves. Periodic motion of a body in physical space is defined by symmetries, whereas periodic motion of divisors…
There are thousands of papers on rootfinding for nonlinear scalar equations. Here is one more, to talk about an apparently new method, which I call ``Inverse Cubic Iteration'' (ICI) in analogy to the Inverse Quadratic Iteration in Richard…
There has been growing interest in high-order tensor methods for nonconvex optimization, with adaptive regularization, as they possess better/optimal worst-case evaluation complexity globally and faster convergence asymptotically. These…
We present a family of analytic solutions for the nearly-equatorial motion of a test particle with precessing spin in Kerr spacetime. We solve the equations of motion up to linear order in the small body's spin for periodic and homoclinic…
Quantum algorithm, as compared to classical algorithm, plays a notable role in solving linear systems of equations with an exponential speedup. Here, we demonstrate a method for solving a particular system of equations by using the concept…
In this paper, a class of optimization problems with nonlinear inequality constraints is discussed. Based on the ideas of sequential quadratic programming algorithm and the method of strongly sub-feasible directions, a new superlinearly…