From: John P. <jwp...@gm...> - 2019-02-20 20:53:21
|
On Wed, Feb 20, 2019 at 2:43 PM Renato Poli <re...@gm...> 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 |