An efficient matrix splitting method for the second-order cone complementarity problem

Given a symmetric and positive (semi)definite $n$-by-$n$ matrix $M$ and a vector, in this paper, we consider the matrix splitting method for solving the second-order cone linear complementarity problem (SOCLCP). The matrix splitting method is among the most widely used approaches for large scale and sparse classical linear complementarity problems (LCP), and its linear convergence is proved by [Luo and Tseng, \textit{SIAM J. Control Optim.}, 30, 408-425 (1992)]. Our first contribution is to prove that, when the general matrix splitting algorithm is applied to SOCLCP with $M$ symmetric and positive definite, it also converges at least linearly. Numerically, our second contribution is to propose a special and efficient matrix splitting algorithm, the block successive over-relaxation (BSOR) method, for solving SOCLCP. The algorithm makes good use of the underlying geometry of SOCLCP and each iteration only involves solving triangular linear systems and requires $O(n^2)$ flops; moreover, the algorithm does not destroy the sparse structure of $M$ and is able to exploit the sparsity effectively. Our code \textsf{BSOR\_BN\_L} based on this method is tested against other four state-of-the-art methods on problems with dense, ill-conditioned symmetric and positive definite matrices, as well as on problems with large scale sparse, symmetric and positive (semi)definite matrices. The numerical results are quite encouraging and \textsf{BSOR\_BN\_L} exhibits a very good performance in terms of its speed, the accuracy and its efficiency in exploiting the sparsity of $M$.

Citation

SIAM Journal on Optimization, to appear.

Article

Download

View An efficient matrix splitting method for the second-order cone complementarity problem