Automatic exploration techniques for the numerical bifurcation study of the Ginzburg-Landau equation

This paper considers the extreme type-II Ginzburg-Landau equations, a nonlinear PDE model that describes the states of a wide range of superconductors. For two-dimensional grids, a robust method is developed that performs a numerical continuation of the equations, automatically exploring the whole solution landscape. The strength of the applied magnetic field is used as the bifurcation parameter. Our branch switching algorithm is based on Lyapunov-Schmidt reduction, but we will show that for an important class of grids an alternative method based on the equivariant branching lemma can be applied as well. The complete algorithm has been implemented in Python and tested for multiple examples. For each example a complete solution landscape was constructed, showing the robustness of the algorithm.