Lehmer–Schur algorithm
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
 |
- Acton, Foreman S. (1970), Numerical Methods That Work (corrected ed.), Washington: Mathematical Association of America (published 1990), p. Chapter 7, ISBN 0-88385-450-3.
- W. H. Press, B. P. Flannery, S. A. Teukolsky, W. T. Vetterling, Numerical Recipes in C: The Art of Scientific Computing, Cambridge University Press, 1992. ISBN 0-521-43108-5 (available free online, with code samples: [1]), section 9.5.