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?
Giang On Tue, Oct 18, 2016 at 4:36 PM, Jed Brown <j...@jedbrown.org> wrote: > Hoang Giang Bui <hgbk2...@gmail.com> writes: > > > Hi Jed > > That's a great idea. Do you have any ex* to start with? Or suggesting me > a > > starting point. > > If you don't have an elasticity solver already, you could sart with > src/snes/examples/tutorials/ex77.c which solves hyperelasticity using a > potentially unstructured mesh. >