I came across this: https://github.com/otherlab/tridiagonal, but haven't tried it. Obviously it can be used for monotone LCPs as well through the convex QP formulation.