Nico Schlömer

Nico Schlömer, Wim Vanroose

This paper considers the extreme type-II Ginzburg--Landau equations, a nonlinear PDE model for describing the states of a wide range of superconductors. Based on properties of the Jacobian operator and an AMG strategy, a preconditioned Newton--Krylov method is constructed. After a finite-volume-type discretization, numerical experiments are done for representative two- and three-dimensional domains. Strong numerical evidence is provided that the number of Krylov iterations is independent of the dimension $n$ of the solution space, yielding an overall solver complexity of O(n).

André Gaul, Nico Schlömer

The authors propose a recycling Krylov subspace method for the solution of a sequence of self-adjoint linear systems. Such problems appear, for example, in the Newton process for solving nonlinear equations. Ritz vectors are automatically extracted from one MINRES run and then used for self-adjoint deflation in the next. The method is designed to work with arbitrary inner products and arbitrary self-adjoint positive-definite preconditioners whose inverse can be computed with high accuracy. Numerical experiments with nonlinear Schrödinger equations indicate a substantial decrease in computation time when recycling is used.