File Exchange

image thumbnail

LCP / MCP solver (Newton-based)

version 1.0.0.0 (2.79 KB) by Yuval
Solve the Linear Complementarity Problem.

6 Downloads

Updated 04 Aug 2008

View License

x = LCP(M,q) solves the LCP

x >= 0
Mx + q >= 0
x'(Mx + q) = 0

x = LCP(M,q,l,u) solves the generalized LCP (a.k.a MCP)

l < x < u => Mx + q = 0
x = u => Mx + q < 0
l = x => Mx + q > 0

x = LCP(M,q,l,u,x0,display) allows the optional initial value 'x0' and
a binary flag 'display' which controls the display of iteration data.

Cite As

Yuval (2021). LCP / MCP solver (Newton-based) (https://www.mathworks.com/matlabcentral/fileexchange/20952-lcp-mcp-solver-newton-based), MATLAB Central File Exchange. Retrieved .

Comments and Ratings (3)

Sara s

pluton schmidt

This does not seem to work?? Considering a simple one-dimensional system with M=1 and q=1 return x=1 and Mx+q=2 (complementarity condition not satisfied) while it should be x=0 and Mx+q=1. I've played with the tolerances without success...

Erez

Thanks!

MATLAB Release Compatibility
Created with R2008a
Compatible with any release
Platform Compatibility
Windows macOS Linux

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!