Full linear multistep methods as root-finders

Root-finders based on full linear multistep methods (LMMs) use previous function values, derivatives and root estimates to iteratively find a root of a nonlinear function. As ODE solvers, full LMMs are typically not zero-stable. However, used as root-finders, the interpolation points are convergent...

Full description

Saved in:
Bibliographic Details
Published in:arXiv.org 2017-09
Main Authors: van Lith, Bart S, Jan H M ten Thije Boonkkamp, IJzerman, Wilbert L
Format: Article
Language:English
Subjects:
Convergence
Interpolation
Mathematical analysis
Polynomials
Robustness (mathematics)
Solvers
Stability analysis
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:Root-finders based on full linear multistep methods (LMMs) use previous function values, derivatives and root estimates to iteratively find a root of a nonlinear function. As ODE solvers, full LMMs are typically not zero-stable. However, used as root-finders, the interpolation points are convergent so that such stability issues are circumvented. A general analysis is provided based on inverse polynomial interpolation, which is used to prove a fundamental barrier on the convergence rate of any LMM-based method. We show, using numerical examples, that full LMM-based methods perform excellently. Finally, we also provide a robust implementation based on Brent's method that is guaranteed to converge.
ISSN:2331-8422