To reduce the communication among processors and improve the computing time for solving linear complementarity problems, we present a two-step modulus-based synchronous multisplitting iteration method and the corresponding symmetric modulus-based multisplitting relaxation methods. The convergence theorems are established when the system matrix is an H+-matrix, which improve the existing convergence theory. Numerical results show that the symmetric modulus-based multisplitting relaxation methods are effective in actual implementation.Mathematics subject classification: 65F10, 68W10, 90C33. 101 the communication by making full use of the previous iteration and communication. To this end, we present the two-step modulus-based synchronous multisplitting iteration methods as well as their relaxed variants in this paper, which consist of two sweeps at each iteration step. We remark that these two-step methods are different from the two-stage methods presented in [8,27], which are inner/outer iteration methods aimed to solve the outer iteration efficiently.The remaining part of this paper is organized as follows: In Section 2, we introduce some notations and briefly review the MSM iteration methods. In Section 3, we propose the two-step modulus-based synchronous multisplitting iteration methods as well as their relaxed variants. In Section 4, we prove their convergence when the system matrix is an H + -matrix. Numerical results are given in Section 5. Finally, we make a conclusion in Section 6.