On 6/1/10 7:19 PM, Roy Stogner wrote: > > On Tue, 1 Jun 2010, Boyce Griffith wrote: > >> I guess there must be some other necessary DOFs in the default send_list? > > Not for enforce_constraints_exactly(), there shouldn't...
Well, using the full send list seems to work correctly, at least in the sense that the assertion does not fail. > In the library it's a little complicated because we have to be sure to > expand out recursive constraints before adding them to the send_list, > but that gets done in the DofConstraints structure and > EquationSystems::reinit() should have handled it for you already. > It's possible to break this system (adding constraints manually > outside of System::user_constrain(), for instance) but not easily, and > I can't figure out how it might have been broken in such a way as to > trip an assert on your special vectors but not on the default vectors. > Your vectors just have non-standard ghost DoFs, right, not > non-standard local DoF ranges? Only non-standard ghost DOFs --- I am leaving the local DOFs alone (at least for now). This is for a problem which uses AMR, but this assertion failure occurs before any mesh refinement has taken place, so the mesh is a simple 3D uniform grid. Are there methods I could call to try to ensure that all of the DofConstraints stuff has been (re-)initialized? -- Boyce ------------------------------------------------------------------------------ _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
