A fast method for a generalized nonlocal elastic model

By carefully exploring the structure of the coefficient matrix of the numerical method, we develop a preconditioned fast Krylov subspace method, which reduces the computations to (Nlog⁡N) per iteration and the memory to O(N). The use of the preconditioner significantly reduces the number of iterations, and the preconditioner can be inverted in O(Nlog⁡N) computations. Numerical results show the utility of the method.