“…In the end, we implemented the difference scheme through two numerical tests, which showed a perfect consistency with our theoretical findings. In addition, our numerical method and the analysis technique of the optimal pointwise error estimates can be easily extended to the cases with spatial fourth-order accuracy [54], the 2D coupled SFNSE [40], 2D space fractional Ginzburg-Landau equation [37,53] and some other space fractional diffusion equation in 2D and 3D [15,51,52,56], which will be our future work. In addition, for the resulting systems of algebraic equations, the coefficient matrices have Toeplitz structure, which can be solved by adopting a super-fast solver with preconditioner [23,[27][28][29] or multigrid methods [36,38] to reduce the CPU time and storage requirement in future.…”