On Fri, Feb 22, 2019 at 2:13 PM Renato Poli <rebp...@gmail.com> wrote:
> The good news is that, if I "reinit" the equation systems in the first > timestep only, it works beautifully. > However, as I rely on this call to change the boundary conditions, that > means I am not able to change the BCs during the run. > > If it is a bug, it looks like workaroundable? > Perhaps forcing reinitialization of some stuff? > What do you think? > If it's a bug, I think we should try to fix it. For that it will be helpful to have a debug-mode stack trace or even better a minimum working example that demonstrates the error. It's possible we are forgetting to close some vector after constraints are applied, I really don't know... -- John _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
