Abstract
Path tracking is the fundamental computational tool in homotopy continuation and is therefore key in most algorithms in the emerging field of numerical algebraic geometry. Though the basic notions of predictor-corrector methods have been known for years, there is still much to be considered, particularly in the specialized algebraic setting of solving polynomial systems. In this article, the effects of the choice of predictor method on the performance of a tracker is analyzed, and details for using Runge-Kutta methods in conjunction with adaptive precision are provided. These methods have been implemented in the Bertini software package, and several examples are described.
Similar content being viewed by others
References
Allgower, E.L., Georg, K.: Numerical continuation methods, An introduction. Springer Series in Computational Mathematics, vol. 13. Springer, Berlin (1990)
Bates, D.J., Hauenstein, J.D., Sommese, A.J., Wampler, C.W.: Bertini: Software for Numerical Algebraic Geometry. Available at http://www.nd.edu/∼sommese/bertini
Bates, D.J., Hauenstein, J.D., Sommese, A.J., Wampler, C.W.: Adaptive multiprecision path tracking. SIAM J. Numer. Anal., 46, 722–746 (2008)
Bates, D.J., Hauenstein, J.D., Sommese, A.J., Wampler, C.W.: Stepsize control for adaptive multiprecision path tracking. Contemp. Math. 496, 21–31 (2009)
Cash, J.R., Karp, A.H.: A variable order Runge-Kutta method for initial value problems with rapidly varying right-hand sides. ACM Trans. Math. Software, 16(3), 201–222 (1990)
Enright, W.H., Jackson, K.R., Nørsett, S.P., Thomsen, P.G.: Interpolants for Runge-Kutta formulas. ACM Trans. Math. Software, 12(3), 193–218 (1986)
Fehlberg, E.: Klassische Runge-Kutta-Formeln fünfter und siebenter Ordnung mit Schrittweiten-Kontrolle. Computing (Arch. Elektron. Rechnen), 4, 93–106 (1969)
Fehlberg, E.: Klassische Runge-Kutta-Formeln vierter und niedrigerer Ordnung mit Schrittweiten-Kontrolle und ihre Anwendung auf Wärmeleitungsprobleme. Computing (Arch. Elektron. Rechnen), 6, 65–71 (1970)
Kincaid, D., Cheney, W.: Numerical Analysis: Mathematics of Scientific Computing, 3rd edn. Brooks/Cole Publishing Co., Pacific Grove, CA (2002)
Li, T.Y.: Numerical solution of polynomial systems by homotopy continuation methods. In: Cucker, F. (ed.) Handbook of Numerical Analysis, Volume XI, Special Volume: Foundations of Computational Mathematics, North-Holland, pp. 209–304 (2003)
Morgan, A.P.: Solving polynomial systems using continuation for engineering and scientific problems. Prentice Hall Inc., Englewood Cliffs, NJ (1987) Reprinted as Classics in Applied Mathematics (2009) 57, SIAM
Morgan, A.P., Sommese, A.J.: Computing all solutions to polynomial systems using homotopy continuation. Appl. Math. Comput. 24(2), 115–138 (1987) Errata: Appl. Math. Comput., 51, 209 (1992)
Prince, P.J., Dormand, J.R.: High order embedded Runge-Kutta formulae. J. Comput. Appl. Math. 7(1), 67–75 1981
Shampine, L.F.: Numerical Solution of Ordinary Differential Equations. Chapman & Hall, New York (1994)
Sommese, A.J., Wampler, C.W.: The Numerical Solution of Systems of Polynomials Arising in Engineering and Science. World Scientific, Singapore (2005)
Verner, J.H.: Explicit Runge-Kutta methods with estimates of the local truncation error. SIAM J. Numer. Anal. 15(4), 772–790 (1978)
Wampler, C.W., Morgan, A., Sommese, A.J.: Complete solution of the nine-point path synthesis problem for four-bar linkages. ASME J. Mech. Des. 114(1), 153–159 (1992)
Author information
Authors and Affiliations
Corresponding author
Additional information
D. J. Bates was supported by the Institute for Mathematics and its Applications (IMA), Colorado State University, and NSF DMS-0914674.
J. D. Hauenstein was supported by the Fields Institute, Texas A&M University, the Duncan Chair of the University of Notre Dame, and NSF grant DMS-0712910.
A. J. Sommese was supported by the Duncan Chair of the University of Notre Dame and NSF grant DMS-0712910.
Rights and permissions
About this article
Cite this article
Bates, D.J., Hauenstein, J.D. & Sommese, A.J. Efficient path tracking methods. Numer Algor 58, 451–459 (2011). https://doi.org/10.1007/s11075-011-9463-8
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11075-011-9463-8
Keywords
- Path tracking
- Homotopy continuation
- Numerical algebraic geometry
- Polynomial systems
- Ordinary differential equations
- Euler’s method
- Runge-Kutta methods
- Precision
- Adaptive precision