Paper detail

Fast Solutions to Projective Monotone Linear Complementarity Problems

We present a new interior-point potential-reduction algorithm for solving monotone linear complementarity problems (LCPs) that have a particular special structure: their matrix $M\in{\mathbb R}^{n\times n}$ can be decomposed as $M=ΦU + Π_0$, where the rank of $Φ$ is $k<n$, and $Π_0$ denotes Euclidean projection onto the nullspace of $Φ^\top$. We call such LCPs projective. Our algorithm solves a monotone projective LCP to relative accuracy $ε$ in $O(\sqrt n \ln(1/ε))$ iterations, with each iteration requiring $O(nk^2)$ flops. This complexity compares favorably with interior-point algorithms for general monotone LCPs: these algorithms also require $O(\sqrt n \ln(1/ε))$ iterations, but each iteration needs to solve an $n\times n$ system of linear equations, a much higher cost than our algorithm when $k\ll n$. Our algorithm works even though the solution to a projective LCP is not restricted to lie in any low-rank subspace.