On Wed, Dec 14, 2011 at 20:03, Dmitry Karpeev <karpeev at mcs.anl.gov> wrote:
> If you don't care about the extra storage, you could put the conserved > variables and the "constraint variables" in separate splits, make the > Jacobian a MatNest and solve with the consevative subblock most of the time > (away from the feasible set boundary). That submatrix would retain the > nice structure and all. I'm thinking about adding "bounded variables" to > libMesh in that way. In TS, we formally add variables to reach a standard form, then rearrange and eliminate on paper to get to a computationally convenient form. I wonder if with some thought, we could do the same for SNESVI. -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20111214/c4e39495/attachment.html>