Abstract: In order to save memory in the stiffness matrix and improve the global search efficiency for nodes and integral points, in this paper we propose an improved algorithm for the element free Galerkin method by using the idea of crossing interacting node pairs to assemble a stiffness matrix, and using the triangular integral mesh to carry out integral calculation. Meanwhile, the CSR format is employed to store stiffness matricies, and a local search method is presented for searching nodes and integral points. The memory requirements of the stiffness matrix and the search efficiency of nodes and integral points for different sized meshes are compared by numerical examples. The results obtained show that the improved algorithm can save memory effectively and improve the search efficiency of nodes and integral points, and is well adapted to complex geometry models.