Stable large-scale solver for Ginzburg-Landau equations for superconductors