I'm using GhostingFunctor for a contact solve, in which I consider a 1/4
domain with partial Dirichlet boundary conditions that impose a symmetry
condition (i.e. displacement normal to the symmetry boundary is clamped to
zero, and tangential displacement is unconstrained). This means that I have
Dirichlet constraints that affect the dofs on the contact surface.
What I find is that the solve works fine in serial, but in parallel the
nonlinear convergence fails, presumably because of an incorrect Jacobian. I
have actually run into this exact issue before (when I was augmenting the
sparsity pattern "manually", prior to GhostingFunctor) and I found that the
issue was that the dof constraints on the contact surface were not being
communicated in parallel, which caused the incorrect Jacobian in parallel.
I fixed it by adding artificial Edge2 elements into the mesh (like in
systems_of_equations_ex8) to ensure that the dof constraints are
communicated properly in parallel.
So, my question is, is there a way to achieve the necessary dof constraint
communication with the new GhostingFunctor API? I had hoped that using
"add_algebraic_ghosting_functor" would do the job, but it apparently
doesn't.
Thanks,
David
------------------------------------------------------------------------------
Check out the vibrant tech community on one of the world's most
engaging tech sites, SlashDot.org! http://sdm.link/slashdot
_______________________________________________
Libmesh-devel mailing list
Libmesh-devel@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-devel
