Projected Splitting Methods for Vertical Linear Complementarity Problems

In this paper, we generalize the projected Jacobi and the projected Gauss–Seidel methods to vertical linear complementarity problems (VLCPs) characterized by matrices with positive diagonal entries. First, we formulate the methods and show that the subproblems that must be solved at each iteration have an explicit solution, which is easy to compute. Then, we prove the convergence of the proposed procedures when the matrices of the problem satisfy some assumptions of strict or irreducible diagonal dominance. In this context, for simplicity, we first analyze the convergence in the special case of VLCPs of dimension $$2n\times n$$ , and we then generalize the results to VLCPs of an arbitrary dimension $$\ell n\times n$$ . Finally, we provide several numerical experiments (involving both full and sparse matrices) that show the effectiveness of the proposed approaches. In this context, our methods are compared with existing solution methods for VLCPs. A parallel implementation of the projected Jacobi method in CUDA is also presented and analyzed.