I think we should: 1. change SNESVI to support general complementarity constraints because only doing box constraints on state variables is lame
2. have a mode to solve the system in the full space, instead of eliminating the semi-smooth variables (which produces the bad conditioning) 3. eventually allow the dimension of the constraint space to change on each SNES iteration (would be fantastic for implicit contact) -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120730/92431593/attachment.html>
