Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userSolving Keplers equation with CORDIC double iterations
56
0
0.0
(
0
)
Added by
Mathias Zechmeister
Publication date
2020
fields
Physics
and research's language is
English
Authors
Mathias Zechmeister
Ask ChatGPT about the research
No Arabic abstract
In a previous work, we developed the idea to solve Keplers equation with a CORDIC-like algorithm, which does not require any division, but still multiplications in each iteration. Here we overcome this major shortcoming and solve Keplers equation using only bitshifts, additions, and one initial multiplication. We prescale the initial vector with the eccentricity and the scale correction factor. The rotation direction is decided without correction for the changing scale. We find that double CORDIC iterations are self-correcting and compensate possible wrong rotations in subsequent iterations. The algorithm needs 75% more iterations and delivers the eccentric anomaly and its sine and cosine terms times the eccentricity. The algorithm can be adopted for the hyperbolic case, too. The new shift-and-add algorithm brings Keplers equation close to hardware and allows to solve it with cheap and simple hardware components.
rate research
Read More
We investigate the use of parabolic equation (PE) methods for solving radio-wave propagation in polar ice. PE methods provide an approximate solution to Maxwells equations, in contrast to full-field solutions such as finite-difference-time-domain (FDTD) methods, yet provide a more complete model of propagation than simple geometric ray-tracing (RT) methods that are the current state of the art for simulating in-ice radio detection of neutrino-induced cascades. PE are more computationally efficient than FDTD methods, and more flexible than RT methods, allowing for the inclusion of diffractive effects, and modeling of propagation in regions that cannot be modeled with geometric methods. We present a new PE approximation suited to the in-ice case. We conclude that current ray-tracing methods may be too simplistic in their treatment of ice properties, and their continued use could overestimate experimental sensitivity for in-ice neutrino detection experiments. We discuss the implications for current in-ice Askaryan-type detectors and for the upcoming Radar Echo Telescope; two families of experiments for which these results are most relevant. We suggest that PE methods be investigated further for in-ice radio applications.
In this paper we present a framework which provides an analytical (i.e., infinitely differentiable) transformation between spatial coordinates and orbital elements for the solution of the gravitational two-body problem. The formalism omits all singular variables which otherwise would yield discontinuities. This method is based on two simple real functions for which the derivative rules are only required to be known, all other applications -- e.g., calculating the orbital velocities, obtaining the partial derivatives of radial velocity curves with respect to the orbital elements -- are thereafter straightforward. As it is shown, the presented formalism can be applied to find optimal instants for radial velocity measurements in transiting exoplanetary systems to constrain the orbital eccentricity as well as to detect secular variations in the eccentricity or in the longitude of periastron.
Recently, there has been significant progress in solving quantum many-particle problem via machine learning based on the restricted Boltzmann machine. However, it is still highly challenging to solve frustrated models via machine learning, which has not been demonstrated so far. In this work, we design a brand new convolutional neural network (CNN) to solve such quantum many-particle problems. We demonstrate, for the first time, of solving the highly frustrated spin-1/2 J$_1$-J$_2$ antiferromagnetic Heisenberg model on square lattices via CNN. The energy per site achieved by the CNN is even better than previous string-bond-state calculations. Our work therefore opens up a new routine to solve challenging frustrated quantum many-particle problems using machine learning.
Calculation of band structure of three dimensional photonic crystals amounts to solving large-scale Maxwell eigenvalue problems, which are notoriously challenging due to high multiplicity of zero eigenvalue. In this paper, we try to address this problem in such a broad context that band structure of three dimensional isotropic photonic crystals with all 14 Bravais lattices can be efficiently computed in a unified framework. We uncover the delicate machinery behind several key results of our work and on the basis of this new understanding we drastically simplify the derivations, proofs and arguments in our framework. In this work particular effort is made on reformulating the Bloch boundary condition for all 14 Bravais lattices in the redefined orthogonal coordinate system, and establishing eigen-decomposition of discrete partial derivative operators by systematic use of commutativity among them, which has been overlooked previously, and reducing eigen-decomposition of double-curl operator to the canonical form of a 3x3 complex skew-symmetric matrix under unitary congruence. With the validity of the novel nullspace free method in the broad context, we perform some calculations on one benchmark system to demonstrate the accuracy and efficiency of our algorithm.
The Bethe-Salpeter equation plays a crucial role in understanding the physics of correlated fermions, relating to optical excitations in solids as well as resonances in high-energy physics. Yet, it is notoriously difficult to control numerically, typically requiring an effort that scales polynomially with energy scales and accuracy. This puts many interesting systems out of computational reach. Using the intermediate representation and sparse modelling for two-particle objects on the Matsubara axis, we develop an algorithm that solves the Bethe-Salpeter equation in $O(L^8)$ time with $O(L^4)$ memory, where $L$ grows only logarithmically with inverse temperature, bandwidth, and desired accuracy, This opens the door for computations in hitherto inaccessible regimes. We benchmark the method on the Hubbard atom and on the multi-orbital weak-coupling limit, where we observe the expected exponential convergence to the analytical results. We then showcase the method for a realistic impurity problem.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
