On Wed, Feb 20, 2019 at 2:43 PM Renato Poli <rebp...@gmail.com> wrote:
> Hi all, > > Just refreshing this one, because I am sort of stuck in inserting a "rigid" > BCs. > It seems simple, but I cannot make it work. > Hi, >From your error message, it sounds like you are introducing a coupling (through the constraint) between two dofs which would otherwise not be coupled, and this causes a new nonzero to be inserted into the system matrix after preallocation. You may be able to work around this issue by augmenting the sparsity pattern using a GhostingFunctor, similar to what is done in miscellaneous/miscellaneous_ex9/augment_sparsity_on_interface.h. -- John _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
