GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems
SIAM Journal on Scientific and Statistical Computing
SIAM Journal on Scientific and Statistical Computing
Ground-state solution of Bose--Einstein condensate by directly minimizing the energy functional
Journal of Computational Physics
Iterative Methods for Sparse Linear Systems
Iterative Methods for Sparse Linear Systems
Computing the Ground State Solution of Bose--Einstein Condensates by a Normalized Gradient Flow
SIAM Journal on Scientific Computing
A Relaxation Scheme for the Nonlinear Schrödinger Equation
SIAM Journal on Numerical Analysis
Journal of Computational Physics
Journal of Computational Physics
A minimisation approach for computing the ground state of Gross-Pitaevskii systems
Journal of Computational Physics
Journal of Computational Physics
Journal of Computational Physics
Journal of Computational and Applied Mathematics
A New Sobolev Gradient Method for Direct Minimization of the Gross-Pitaevskii Energy with Rotation
SIAM Journal on Scientific Computing
Hi-index | 31.45 |
We consider the Backward Euler SPectral (BESP) scheme proposed in [10] for computing the stationary states of Bose-Einstein Condensates (BECs) through the Gross-Pitaevskii equation. We show that the fixed point approach introduced in [10] fails to converge for fast rotating BECs. A simple alternative approach based on Krylov subspace solvers with a Laplace or Thomas-Fermi preconditioner is given. Numerical simulations (obtained with the associated freely available Matlab toolbox GPELab) for complex configurations show that the method is accurate, fast and robust for 2D/3D problems and multi-components BECs.