Scalable higher-order nonlinear solvers via higher-order automatic differentiation
Abstract
This paper demonstrates new methods and implementations of nonlinear solvers with higher-order of convergence, which is achieved by efficiently computing higher-order derivatives. Instead of computing full derivatives, which could be expensive, we compute directional derivatives with Taylor-mode automatic differentiation. We first implement Householder's method with arbitrary order for one variable, and investigate the trade-off between computational cost and convergence order. We find that the second-order variant, i.e., Halley's method, to be the most valuable, and further generalize Halley's method to systems of nonlinear equations and demonstrate that it can scale efficiently to large-scale problems. We further apply Halley's method on solving large-scale ill-conditioned nonlinear problems, as well as solving nonlinear equations inside stiff ODE solvers, and demonstrate that it could outperform Newton's method.
Cite
@article{arxiv.2501.16895,
title = {Scalable higher-order nonlinear solvers via higher-order automatic differentiation},
author = {Songchen Tan and Keming Miao and Alan Edelman and Christopher Rackauckas},
journal= {arXiv preprint arXiv:2501.16895},
year = {2025}
}