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... 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? --- Roy ------------------------------------------------------------------------------ _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
