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