Abstract
The Trilinos library LOCA (http://www.cs.sandia.gov/LOCA/) allows computing branches of steady states of large-scale dynamical systems like (discretized) nonlinear PDEs. The core algorithms typically are (pseudo-)arclength continuation, Newton–Krylov methods and (sparse) eigenvalue solvers. While LOCA includes some basic techniques for computing bifurcation points and switching branches, the exploration of a complete bifurcation diagram still takes a lot of programming effort and manual interference.
On the other hand, recent developments in algorithms for fully automatic exploration are condensed in PyNCT (https://pypi.org/project/PyNCT/). The scope of this algorithmically versatile software is, however, limited to relatively small (e.g. 2D) problems because it relies on linear algebra from Python libraries like NumPy. Furthermore, PyNCT currently does not support problems with a non-Hermitian Jacobian matrix, which rules out interesting applications in chemistry and fluid dynamics.
In this paper we aim to combine the best of both worlds: a high-level implementation of algorithms in PyNCT with parallel models and linear algebra implemented in Trilinos. PyNCT is extended to non-symmetric systems and its complete backend is replaced by the PHIST library (https://bitbucket.org/essex/phist), which allows us to use the same underlying HPC libraries as LOCA does.
We then apply the new code to a reaction-diffusion model to demonstrate its potential of enabling fully automatic bifurcation analysis on parallel computers.
This is a preview of subscription content, log in via an institution to check access.
Access this chapter
Tax calculation will be finalised at checkout
Purchases are for personal use only
Similar content being viewed by others
References
S. Balay, S. Abhyankar, M. F. Adams, J. Brown, P. Brune, K. Buschelman, L. Dalcin, V. Eijkhout, W. D. Gropp, D. Kaushik, M. G. Knepley, L. C. McInnes, K. Rupp, B. F. Smith, S. Zampini, and H. Zhang. PETSc Web page. http://www.mcs.anl.gov/petsc, 2016.
H. A. Dijkstra, F. W. Wubs, A. K. Cliffe, E. Doedel, I. F. Dragomirescu, B. Eckhardt, A. Y. Gelfgat, A. L. Hazel, V. Lucarini, A. G. Salinger, E. T. Phipps, J. Sanchez-Umbria, H. Schuttelaars, L. S. Tuckerman, and U. Thiele. Numerical bifurcation methods and their application to fluid dynamics: Analysis beyond simulation. Communications in Computational Physics, 15(01):1–45, Jan 2014.
D. Draelants, D. Avitabile, and W. Vanroose. Localized auxin peaks in concentration-based transport models of the shoot apical meristem. J. R. Soc. Interface, 12, May 2015.
D. Draelants, J. Broeckhove, G. Beemster, and W. Vanroose. Numerical bifurcation analysis of pattern formation in a cell based auxin transport model. J. Math. Biol., 67, Sep 2012.
D. Draelants, P. Klosiewicz, J. Broeckhove, and W. Vanroose. Solving general auxin transport models with a numerical continuation toolbox in Python: PyNCT. Lect. Notes Comput. Sci., 9271:211–225, 2015.
M. Golubitsky and I. Stewart. The symmetry perspective. Birkhäuser, 2002.
M. Heroux, R. Bartlett, V. H. R. Hoekstra, J. Hu, T. Kolda, R. Lehoucq, K. Long, R. Pawlowski, E. Phipps, A. Salinger, H. Thornquist, R. Tuminaro, J. Willenbring, and A. Williams. An overview of Trilinos. Technical Report SAND2003-2927, Sandia National Laboratories, 2003.
H. B. Keller. Numerical solution of bifurcation and nonlinear eigenvalue problems. In P. H. Rabinowitz, editor, Applications of Bifurcation Theory, pages 359–384. Academic Press, New York, U.S.A., 1977.
T. Leppänen. Computational studies of pattern formation in Turing systems, 2004-11-27.
W. Song, F. Wubs, J. Thies, and S. Baars. Numerical bifurcation analysis of a 3D Turing-type reaction-diffusion model. Communications in nonlinear science and numerical simulation, 60:145–164, 7 2018.
J. Thies, M. Röhrig-Zöllner, N. Overmars, A. Basermann, D. Ernst, G. Hager, and G. Wellein. PHIST: a pipelined, hybrid-parallel iterative solver toolkit. ACM Transactions on Mathematical Software, 46(4), 2020.
M. Wouters. Automatic exploration techniques for the numerical continuation of large–scale nonlinear systems. PhD thesis, Universiteit Antwerpen, 2019.
M. Wouters and W. Vanroose. Automatic exploration techniques of numerical bifurcation diagrams illustrated by the Ginzburg–Landau equation. SIAM Journal on Applied Dynamical Systems, 18(4):2047–2098, 2019.
Author information
Authors and Affiliations
Corresponding author
Editor information
Editors and Affiliations
Rights and permissions
Copyright information
© 2021 Springer Nature Switzerland AG
About this paper
Cite this paper
Thies, J., Wouters, M., Hennig, RS., Vanroose, W. (2021). Towards Scalable Automatic Exploration of Bifurcation Diagrams for Large-Scale Applications. In: Vermolen, F.J., Vuik, C. (eds) Numerical Mathematics and Advanced Applications ENUMATH 2019. Lecture Notes in Computational Science and Engineering, vol 139. Springer, Cham. https://doi.org/10.1007/978-3-030-55874-1_97
Download citation
DOI: https://doi.org/10.1007/978-3-030-55874-1_97
Published:
Publisher Name: Springer, Cham
Print ISBN: 978-3-030-55873-4
Online ISBN: 978-3-030-55874-1
eBook Packages: Mathematics and StatisticsMathematics and Statistics (R0)