The concept of candidate particle set is introduced in the MPS gridless numerical method to generate neighboring particle setmatrix,which can reduce the CPU time to 1/11 of that before introduction.The Bi-CGSTAB (bi-conjugate gradient stabilized) algorithm is applied to solving the Poisson pressure equation,by which the solving speed is significantly accelerated.The process of solitary waves propagating over a numerical flume and interacting with a vertical wall is simulated.The simulated results of water surface elevation are in good agreement with the analytical solution as well as the measured data.The predicted maximum values of the run-up of solitary waves with various relative incident wave heights agree well with the measured results.