[petsc-users] preconditioner for contact / mesh tying problem

[Jed Brown]

Hoang Giang Bui <hgbk2008 at gmail.com> writes:

> I do have an elasticity solver though. However it's an in-house code hence
> 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.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 800 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20161018/804e51aa/attachment.pgp>