Composing Scalable Nonlinear Algebraic Solvers

Abstract

Most efficient linear solvers use composable algorithmic components, with the most common model being the combination of a Krylov accelerator and one or more preconditioners. A similar set of concepts may be used for nonlinear algebraic systems, where nonlinear composition of different nonlinear solvers may significantly improve the time to solution. We describe the basic concepts of nonlinear composition and preconditioning and present a number of solvers applicable to nonlinear partial differential equations. We have developed a software framework in order to easily explore the possible combinations of solvers. We show that the performance gains from using composed solvers can be substantial compared with gains from standard Newton-Krylov methods.

29 pages, 14 figures, 13 tables

References in corpus (1)

• Mayavi: a package for 3D visualization of scientific data

Cited by in corpus (26)

• An augmented Lagrangian preconditioner for the 3D stationary incompressible Navier-Stokes equations at high Reynolds number
• Anderson acceleration for contractive and noncontractive operators
• Finite element interpolated neural networks for solving forward and inverse problems
• Structural Landscapes in Geometrically Frustrated Smectics
• PCPATCH: software for the topological construction of multigrid relaxation methods
• The divergence-conforming immersed boundary method: Application to vesicle and capsule dynamics
• Enhancing training of physics-informed neural networks using domain-decomposition based preconditioning strategies
• Nektar++: Design and implementation of an implicit, spectral/$hp$ element, compressible flow solver using a Jacobian-free Newton Krylov approach
• Pipelined, Flexible Krylov Subspace Methods
• Dissipation-Based Continuation Method for Multiphase Flow in Heterogeneous Porous Media
• A scalable elliptic solver with task-based parallelism for the SpECTRE numerical relativity code
• An inexact Newton-Krylov method for stochastic eigenvalue problems
• Deflation for semismooth equations
• Nonequilibrium thermomechanics of Gaussian phase packet crystals: application to the quasistatic quasicontinuum method
• Robust and Efficient Multilevel-ILU Preconditioning of Hybrid Newton-GMRES for Incompressible Navier-Stokes Equations
• A full approximation scheme multilevel method for nonlinear variational inequalities
• Anderson Acceleration as a Krylov Method with Application to Asymptotic Convergence Analysis
• Robust parallel nonlinear solvers for implicit time discretizations of the Bidomain equations
• Objective acceleration for unconstrained optimization
• An augmented Lagrangian preconditioner for implicitly-constituted non-Newtonian incompressible flow
• Numerical Analysis of Unsteady Implicitly Constituted Incompressible Fluids: Three-Field Formulation
• Nesterov Acceleration of Alternating Least Squares for Canonical Tensor Decomposition: Momentum Step Size Selection and Restart Mechanisms
• An adaptive discontinuous Petrov-Galerkin method for the Grad-Shafranov equation
• On the Asymptotic Linear Convergence Speed of Anderson Acceleration, Nesterov Acceleration, and Nonlinear GMRES
• A performance spectrum for parallel computational frameworks that solve PDEs
• Iterative convergence in phase-field brittle fracture computations: exact line search is all you need