An Efficient Matrix Splitting Method for the Second-Order Cone Complementarity Problem
Lei-Hong Zhang,Wei Hong Yang
DOI: https://doi.org/10.1137/13090938x
IF: 2.763
2014-01-01
SIAM Journal on Optimization
Abstract:Given a symmetric and positive (semi)definite $n$-by-$n$ matrix $M$and a vector, in this paper, we consider the matrix splittingmethod for solving the second-order cone linear complementarityproblem (SOCLCP). The matrix splitting method is among the mostwidely used approaches for large scale and sparse classical linearcomplementarity problems, and its linear convergence isproved by [Luo and Tseng, SIAM J. Control Optim. , 30 (1992),pp. 408--425]. Our first contribution is to prove that, when thegeneral matrix splitting algorithm is applied to SOCLCP with$M$ symmetric and positive definite, it also converges at leastlinearly. Numerically, our second contribution is to propose aspecial and efficient matrix splitting algorithm, the blocksuccessive overrelaxation method, for solving the SOCLCP.The algorithm makes good use of the underlying geometry ofthe SOCLCP and each iteration only involves solving triangular linearsystems and requires $O(n^2)$ flops; moreover, the algorithm doesnot destroy the sparse structure of $M$ and is able to exploit thesparsity effectively. Our code BSOR_BN_L based on thismethod is tested against four other state-of-the-art methods on problems with dense, ill-conditionedsymmetric and positive definite matrices, as well as on problemswith large scale sparse, symmetric and positive (semi)definite matrices. The numericalresults are quite encouraging and BSOR_BN_L exhibits avery good performance in terms of its speed, accuracy, and efficiency in exploiting the sparsity of $M$.