
<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Hmmmm..sorry guys - the entire thread got put away in a folder by my mail reader and I just discovered all your emails :-)<div><br></div><div>SNESVI looks interesting but as someone pointed out I am using matrix free. And as I understand from Barry a matrix free version of SNESVI is not implemented. The idea of flipping the sign seems to be a poor mans approach that might work for me though I would still need some mechanisms in SNES to put the hooks in.</div><div><br></div><div>Bobby</div><div><br><div><div>On Dec 14, 2011, at 10:45 PM, Jed Brown wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div class="gmail_quote">On Wed, Dec 14, 2011 at 19:40, Dmitry Karpeev <span dir="ltr"><<a href="mailto:karpeev@mcs.anl.gov">karpeev@mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

The trouble is that the constraints can get "eliminated" only when they become active. </blockquote></div><br><div>That's not what I meant. I meant to transform the algebraic system so that those extra variables were eliminated. The point is that we tend to put a lot of effort into designing effective preconditioners for a standard formulation (e.g. conservative variables), and that is partly lost of we have this other system in which the evaluations of constitutive relations are added explicitly.</div>

</blockquote></div><br></div></body></html>