Journal
APPLIED MATHEMATICS AND COMPUTATION
Volume 338, Issue -, Pages 758-773Publisher
ELSEVIER SCIENCE INC
DOI: 10.1016/j.amc.2018.06.010
Keywords
Nonlinear fractional Schrodinger equation; Finite element method; Krylov subspace method; Toeplitz matrix; Fast Fourier transform; Circulant preconditioner
Categories
Funding
- NSF of China [11771163]
- China Postdoctoral Science Foundation [2018M632791]
Ask authors/readers for more resources
The main aim of this paper is to apply the Galerkin finite element method to numerically solve the nonlinear fractional Schrodinger equation with wave operator. We first construct a fully discrete scheme combining the Crank-Nicolson method with the Galerkin finite element method. Two conserved quantities of the discrete system are shown. Meanwhile, the prior bound of the discrete solutions are proved. Then, we prove that the discrete scheme is unconditionally convergent in the senses of L-2-norm and H-alpha/2-norm. Moreover, by the proposed iterative algorithm, some numerical examples are given to verify the theoretical results and show the effectiveness of the numerical scheme. Finally, a fast Krylov subspace solver with suitable circulant preconditioner is designed to solve above Toeplitz-like linear system. In each iterative step, this method can effectively reduce the memory requirement of the proposed iterative finite element scheme from O(M-2) to O(M), and the computational complexity from O(M-3) to O(MlogM), where M is the number of grid nodes. Several numerical tests are carried out to show that this fast algorithm is more practical than the traditional backslash and LU factorization methods, in terms of memory requirement and computational cost. (C) 2018 Elsevier Inc. All rights reserved.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available