Clear on that. First I need some time to study ex77 carefully.
On Tue, Oct 18, 2016 at 6:19 PM, Jed Brown <j...@jedbrown.org> wrote:
> Hoang Giang Bui <hgbk2...@gmail.com> writes:
> > I do have an elasticity solver though. However it's an in-house code
> > it's not very straight forward to unroll to make a PETSc example. The
> > example can only read in the provided matrices and apply the
> > preconditioner. However, I think in general standard PETSc example
> > generates the linear system and solve it successively. In that way, ex77
> > must be extended. What do you think?
> I suggest SNES ex77 (versus writing a pure KSP example) because the
> contact problem is nonlinear and your example might as well actually
> solve a problem instead of just solving one linearized step. That said,
> VIRS isn't set up for general constraints using Lagrange multipliers, so
> it might be trying to tackle too many problems at once.