Algorithm 956: PAMPAC, a parallel adaptive method for pseudo-arclength continuation
From MaRDI portal
Abstract: Pseudo-arclength continuation is a well-established method for generating a numerical curve approximating the solution of an underdetermined system of nonlinear equations. It is an inherently sequential predictor-corrector method in which new approximate solutions are extrapolated from previously converged results and then iteratively refined. Convergence of the iterative corrections is guaranteed only for sufficiently small prediction steps. In high-dimensional systems, corrector steps are extremely costly to compute and the prediction step-length must be adapted carefully to avoid failed steps or unnecessarily slow progress. We describe a parallel method for adapting the step-length employing several predictor-corrector sequences of different step lengths computed concurrently. In addition, the algorithm permits intermediate results of unconverged correction sequences to seed new predictions. This strategy results in an aggressive optimization of the step length at the cost of redundancy in the concurrent computation. We present two examples of convoluted solution curves of high-dimensional systems showing that speed-up by a factor of two can be attained on a multi-core CPU while a factor of three is attainable on a small cluster.
Recommendations
- Algorithms for Simultaneous Padé Approximations
- Algorithm 986: A suite of compact finite difference schemes
- scientific article; zbMATH DE number 1183041
- A Parallel Implementation of the p-Version of the Finite Element Method
- A Pseudo-Arclength Continuation Method for Nonlinear Eigenvalue Problems
- A parallel multisplitting iterative algorithm for a partial differential equation
- Parallel algorithms for a singularly perturbed parabolic problem
- scientific article; zbMATH DE number 873752
- scientific article; zbMATH DE number 761757
- Parallel adaptive and time-stabilizing schemes for constant-coefficient parabolic PDE's
Cites work
- scientific article; zbMATH DE number 5305088 (Why is no real title available?)
- scientific article; zbMATH DE number 1178976 (Why is no real title available?)
- A Direct Method for Computing Higher Order Folds
- Condition Estimates for Pseudo‐Arclength Continuation
- Elements of applied bifurcation theory.
- GMRES: A Generalized Minimal Residual Algorithm for Solving Nonsymmetric Linear Systems
- Introduction to Numerical Continuation Methods
- Jacobian-free Newton-Krylov methods: a survey of approaches and applications.
- LAPACK Users' Guide
- Newton-Krylov continuation of periodic orbits for Navier-Stokes flows
- Nonlinear analysis of hydrodynamic instability in laminar flames—I. Derivation of basic equations
- Numerical Methods for Bifurcations of Dynamical Equilibria
- Periodic motion representing isotropic turbulence
- Solving Nonlinear Equations with Newton's Method
- The significance of simple invariant solutions in turbulent flows
Cited in
(5)- A \(J\)-\textit{integral}-based arc-length solver for brittle and ductile crack propagation in finite deformation-finite strain hyperelastic solids with an application to graphene kirigami
- Parallel numerical continuation of periodic responses of local nonlinear systems
- PAMPAC
- Solution space exploration of low-thrust minimum-time trajectory optimization by combining two homotopies
- A new paradigm for multi-fidelity continuation using parallel model refinement
This page was built for publication: Algorithm 956: PAMPAC, a parallel adaptive method for pseudo-arclength continuation
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2828171)
