Abstract

This article develops a new predictor-corrector algorithm for numerical path tracking in the context of polynomial homotopy continuation. In the corrector step, it uses a newly developed Newton corrector algorithm which rejects an initial guess if it is not an approximate zero. The algorithm also uses an adaptive step size control that builds on a local understanding of the region of convergence of Newton’s method and the distance to the closest singularity following Telen, Van Barel, and Verschelde. To handle numerically challenging situations, the algorithm uses mixed precision arithmetic. The efficiency and robustness are demonstrated in several numerical examples.

Highlights

  • Systems of polynomial equations arise in many applications across the sciences including computer vision [20, 36], chemistry [28], kinematics [41], and biology [32]

  • This article introduces a new path tracking algorithm, Algorithm 3, that does not require the choice of a path tracking tolerance ε or a maximal number N of corrector iterations allowed per step

  • The implementation of the proposed algorithm significantly outperforms the available implementation of the adaptive precision algorithm [5] in the the considered examples

Read more

Summary

Introduction

Systems of polynomial equations arise in many applications across the sciences including computer vision [20, 36], chemistry [28], kinematics [41], and biology [32]. A numerical method for finding all isolated solutions of a system F of n polynomials in n unknowns is homotopy continuation [35]. This method constructs a homotopy H (x, t) : Cn × C → Cn such that H (x, t) is a polynomial system for all t, H (x, 1) = F (x) and H (x, 0) is a system whose isolated solutions are known.

75 Page 2 of 23
75 Page 4 of 23
Newton’s method: theory and computational aspects
Convergence results
75 Page 6 of 23
Computational aspects and floating point arithmetic
75 Page 8 of 23
75 Page 10 of 23
Predictors and Padeapproximants
75 Page 12 of 23
75 Page 14 of 23
Step size control and path tracking algorithm
75 Page 16 of 23
Computational experiments
Mixtures of Gaussians and the method of moments
75 Page 18 of 23
Alt’s problem
75 Page 20 of 23
Steiner’s conic problem
Findings
Conclusion and future work
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call