Jump to content

Lehmer–Schur algorithm

From Wikipedia, the free encyclopedia
This is an old revision of this page, as edited by LutzL (talk | contribs) at 13:05, 22 December 2013 (External links). The present address (URL) is a permanent link to this revision, which may differ significantly from the current revision.

In mathematics, the Lehmer–Schur algorithm (named after Derrick Henry Lehmer and Issai Schur) is a root-finding algorithm extending the one-dimensional bracketing used by the bisection method to find the roots of a function of one complex variable inside any rectangular region of the function's holomorphicity (i.e., analyticity).

The rectangle in question is quadrisected into four, congruent quarter rectangles. The argument principle is then applied to the boundary of each quarter to find the winding number (about the origin) for the boundary. Given that the function is analytic within each of these quarters, a nonzero winding number N (always an integer) identifies N zeros of the function inside the quarter in question, each zero counted as many times as its multiplicity.

Analogously with the bisection method, the algorithm is then applied recursively to any quarter whose boundary has nonzero winding number to further refine the estimates of the zeros. The recursion is repeated until the zero-containing rectangles are either small enough that their centres give sufficiently accurate zero estimates or, alternatively, that another root-finding algorithm can be applied to the estimates to further refine them.

References

  • D. H. Lehmer, (1961), "A Machine Method for Solving Polynomial Equations", Journal of the ACM, 8, (2): 151–162, doi:10.1145/321062.321064 {{citation}}: Unknown parameter |month= ignored (help)CS1 maint: extra punctuation (link)
  • Herbert S. Wilf (1978), "A Global Bisection Algorithm for Computing the Zeros of Polynomials in the Complex Plane"", Journal of the ACM, 25 (3), doi:10.1145/322077.322084