This paper presents the localized collocation solver based on fundamental solutions to 3D elastic wave propagation analysis. In the proposed collocation solver, the approximated solution at the considered node is represented by a linear combination of the solutions at a few nearest nodes in the stencil support of the considered node instead of the whole discretization nodes. The related weighting functions are determined by solving a system of linear equations, which is constructed as a linear combination of fundamental solutions in the related stencil support. Therefore, the proposed collocation solver produces sparse resultant matrix, which makes it possible to perform large-scale elastic wave propagation simulations on a desktop computer. Furthermore, it avoids highly ill-conditioned dense matrix encountered in most of collocation methods. The efficiency and accuracy of the proposed method have been verified under several benchmark examples.