On Wed, 2 Mar 2011, Roman Vetter wrote: > Thank you for the detailed description, I have been able to implement > it as you suggested. Both System::ProjectVector::operator() and > System::BuildProjectionList::operator() needed an additional continue > in the case that the old_dof_object doesn't exist for a given element > or it has no dofs for the current system (this is relevant for > equation systems with more than one system).
The patch is committed now, but I'm afraid I don't yet understand under what circumstances the second half of those if statements will trigger. Surely if you add a new element to an existing mesh then it won't have an old_dof_object at all, but if it does have an old_dof_object then there should be valid dofs on it for every System. Or are you looking at adding a new System to a problem as well as new elements? --- Roy ------------------------------------------------------------------------------ What You Don't Know About Data Connectivity CAN Hurt You This paper provides an overview of data connectivity, details its effect on application quality, and explores various alternative solutions. http://p.sf.net/sfu/progress-d2d _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
