"Salazar De Troya, Miguel" <salazardet...@llnl.gov> writes:
> I write the boundary conditions using their DirichletBoundary class, not > the penalty term. > The options I¹m using are ones that I found in the libMesh mail list > from a user who suggested them for elasticity problems. The idea he > mentioned was to use field split to separate each field of the > displacement vector solution. I honestly do not know the role of > -pc_fieldsplit_type symmetric_multiplicative, but it was working for > me. I would use GAMG or ML (without fieldsplit; set a near null space) instead of Hypre, because the algorithm is usually better for elasticity. Fieldsplit can work for this, but the performance degrades for higher Poisson ratio. In any case, "-pc_fieldsplit_0 0,1" is ignored (I think) and "-pc_fieldsplit_0_fields 0,1" is likely not what you want.
signature.asc
Description: PGP signature