Fourier pseudospectral methods for the spatial variable-order fractional wave equations

In this paper, we propose Fourier pseudospectral methods to solve variable-order fractional viscoacoustic wave equations. Our approach involves a Fourier pseudospectral method for spatial discretization and an accelerated matrix-free technique for efficient computation and storage costs, with a computational cost of $\mathcal{O}(MN\log N)$ and storage cost $\mathcal{O}(MN)$ where $M\ll N$. For temporal discretization, we employ the Crank-Nicolson, leap-frog, and time-splitting schemes. Numerical experiments are conducted to assess their performance. The results demonstrate the advantages of our fast method, particularly in computational and storage costs, and its feasibility in high dimensions. The numerical findings reveal that all three temporal discretization methods exhibit second-order accuracy, while the Fourier pseudospectral spatial discretization showcases spectral accuracy.