[petsc-dev] controlling vector values while doing matrix free operations

[Jed Brown]

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>