In this paper, a relaxed two-step modulus-based matrix synchronous multisplitting iteration method for solving the linear complementarity problems is constructed. The convergence conditions of the proposed method are analyzed with the convergence range of the relaxation parameters. Some parallel numerical experiments under OpenACC framework are given to show that the proposed method can accelerate the convergence rate of the existing method significantly.