{
  "id": "1108.5773",
  "version": "v3",
  "published": "2011-08-29T22:37:31.000Z",
  "updated": "2016-10-08T16:25:51.000Z",
  "title": "On the Efficient Global Dynamics of Newton's Method for Complex Polynomials",
  "authors": [
    "Dierk Schleicher"
  ],
  "comment": "Updated and improved; 26 pages, 4 figures, 2 tables. Main change (besides improved exposition): description of practical experiments that lead to complete factorization of polynomials of degrees greater than 100 million",
  "categories": [
    "math.DS",
    "math.NA"
  ],
  "abstract": "We investigate Newton's method as a root finder for complex polynomials of arbitrary degree. While polynomial root finding continues to be one of the fundamental tasks of computing, with essential use in all areas of theoretical mathematics, numerics, computer graphics and physics, known methods have either excellent theoretical complexity but cannot be used in practice, or are practically efficient but are a lacking a successful theory behind them. In this manuscript we describe the theoretical complexity of Newton's method for finding all roots of polynomials of given degree and show that it is near-optimal for the known set of starting points that find all roots. This theoretical result is complemented by a recent implementation of Newton's method that finds all roots of various polynomials of degree more than a million, significantly faster than our upper bounds on the complexity indicate, and often much faster than established fast root finders. In some experiments, it was possible to find all roots even with complexity $O(d\\log d)$ for degrees exceeding 100 million. Newton's method thus stands out as a method that has merits both from the theoretical and from the practical point of view. Our study is based on the known explicit set of universal starting points that is guaranteed to find all roots of polynomials of degree $d$ (appropriately normalized). We show that this set contains $d$ points that converge very quickly to the $d$ roots: the expected total number of Newton iterations required to find all roots with precision $\\eps$ is $O(d^3\\log^3d+d\\log|\\log\\eps|)$, which can be further improved to $O(d^2\\log^4d+d\\log|\\log\\eps|)$; in the worst case allowing near-multiple roots, the complexity is $O(d^4\\log^2d+d^3\\log^2d|\\log\\eps|)$. The arithmetic complexity for all these estimates is the same as the number of Newton iterations steps, up to a factor of $\\log^2 d$.",
  "revisions": [
    {
      "version": "v2",
      "updated": "2013-12-16T11:58:16.000Z",
      "abstract": "We investigate Newton's method for complex polynomials of arbitrary degree $d$, normalized so that all their roots are in the unit disk. We specify an explicit universal set of starting points, consisting of $O(d\\log^2d)$ points and depending only on $d$, so that among them there are $d$ points that converge very quickly to the $d$ roots: we prove that the expected total number of Newton iterations required to find all $d$ roots with precision $\\eps$ is $O(d^3\\log^3d+d\\log|\\log\\eps|)$, which can be further improved to $O(d^2\\log^4d+d\\log|\\log\\eps|)$; in the worst case possibly with near-multiple roots, the complexity is $O(d^3\\log^2d(d+|\\log\\eps|))$. The arithmetic complexity for all these Newton iterations is the same as the number of Newton iterations, up to a factor of $\\log^2 d$.",
      "comment": "Updated and improved; 21 pages, 3 figures",
      "journal": null,
      "doi": null
    },
    {
      "version": "v3",
      "updated": "2016-10-08T16:25:51.000Z"
    }
  ],
  "analyses": {
    "subjects": [
      "37F10",
      "49M15"
    ],
    "keywords": [
      "efficient global dynamics",
      "complex polynomials",
      "newtons method",
      "newton iterations",
      "explicit universal set"
    ],
    "note": {
      "typesetting": "TeX",
      "pages": 26,
      "language": "en",
      "license": "arXiv",
      "status": "editable",
      "adsabs": "2011arXiv1108.5773S"
    }
  }
}