An optimal linear solver for the Jacobian system of the extreme type-II Ginzburg-Landau problem
Abstract
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).
- Publication:
-
Journal of Computational Physics
- Pub Date:
- February 2013
- DOI:
- 10.1016/j.jcp.2012.10.013
- arXiv:
- arXiv:1204.5125
- Bibcode:
- 2013JCoPh.234..560S
- Keywords:
-
- Physics - Computational Physics;
- Mathematics - Numerical Analysis;
- 35Q56;
- G.1.5;
- G.1.8;
- J.2
- E-Print:
- doi:10.1016/j.jcp.2012.10.013