How to solve MLCP with Projected Gauss-Seidel?

Please don't post Bullet support questions here, use the above forums instead.
User avatar
Erwin Coumans
Site Admin
Posts: 4221
Joined: Sun Jun 26, 2005 6:43 pm
Location: California, USA

Re: How to solve MLCP with Projected Gauss-Seidel?

Post by Erwin Coumans »

Dirk Gregorius wrote:@Erwin/Morten: You are refering to non-linerar complementary problems. Where did we move from the LCP to the NCP? Which part is actually non-linear in the ODE/Bullet formulation?
According to Kenny, updating the low/high limits for the friction every iteration makes the problem non-linear. This is also why we cannot prove convergence of PGS anymore.
Thanks,
Erwin
silcowitz
Posts: 2
Joined: Thu Jul 05, 2007 8:25 pm

Re: How to solve MLCP with Projected Gauss-Seidel?

Post by silcowitz »

Its because if some problem is to be an LCP problem, it must be on the form
w=Ax+b>=0, x>=0, w^Tx =0.

So, when we want to model the coupling between normal force magnitudes and friction force magnitudes, the only way to put that into the form of an LCP is to introduce extra variables, like Stewart and Trinkle did (they use lambda as symbol in some of the papers if I recall). When you do that, you end up with a system matrix that cannot by solved by PGS anymore. Then, someone came up with the idea of dropping the friction cone approximation from Stewart-Trinkle, and just projecting the friction forces one by one in the PGS loop. (like its done in bullet?), That showed to work out well, but it is not on the LCP form.

/morten
Dirk Gregorius
Posts: 861
Joined: Sun Jul 03, 2005 4:06 pm
Location: Kirkland, WA

Re: How to solve MLCP with Projected Gauss-Seidel?

Post by Dirk Gregorius »

I see. Thanks for the clarification.

Cheers,
-Dirk
ngbinh
Posts: 117
Joined: Fri Aug 12, 2005 3:47 pm
Location: Newyork, USA

Re: How to solve MLCP with Projected Gauss-Seidel?

Post by ngbinh »

@Erwin: You're right! Thanks for correcting me on that.