Scalable higher-order nonlinear solvers via higher-order automatic differentiation
read the original 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.
This paper has not been read by Pith yet.
Forward citations
Cited by 1 Pith paper
-
NOWS: Neural Operator Warm Starts for Accelerating Iterative Solvers
Neural operators supply warm-start guesses that cut iteration counts and runtime by up to 90% in Krylov solvers for PDEs while retaining the original methods' convergence guarantees.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.