Towards Scalable Automatic Exploration of Bifurcation Diagrams for Large-Scale Applications

Jonas Thies*, Michiel Wouters, Rebekka Sarah Hennig, Wim Vanroose

*Corresponding author for this work

Research output: Chapter in Book/Conference proceedings/Edited volumeConference contributionScientificpeer-review


The Trilinos library 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 ( ). 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 ( ), 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.

Original languageEnglish
Title of host publicationNumerical Mathematics and Advanced Applications, ENUMATH 2019 - European Conference
EditorsFred J. Vermolen, Cornelis Vuik
Number of pages9
ISBN (Print)9783030558734
Publication statusPublished - 2021
Externally publishedYes
EventEuropean Conference on Numerical Mathematics and Advanced Applications, ENUMATH 2019 - Egmond aan Zee, Netherlands
Duration: 30 Sep 20194 Oct 2019

Publication series

NameLecture Notes in Computational Science and Engineering
ISSN (Print)1439-7358
ISSN (Electronic)2197-7100


ConferenceEuropean Conference on Numerical Mathematics and Advanced Applications, ENUMATH 2019
CityEgmond aan Zee


Dive into the research topics of 'Towards Scalable Automatic Exploration of Bifurcation Diagrams for Large-Scale Applications'. Together they form a unique fingerprint.

Cite this