     * Fills the vector \p di with the global degree of freedom indices
     * for the element. If no variable number is specified then all
     * variables are returned.
-@@ -979,14 +977,18 @@
+@@ -979,14 +987,18 @@
    void add_constraints_to_send_list();
  
  #endif // LIBMESH_ENABLE_CONSTRAINTS

 
       DofObject::index_buffer_t::iterator
 	it  = _idx_buf.begin(),
@@ -243,27 +247,27 @@
       _idx_buf.erase(it,end);
 
       for (unsigned int ctr=(s+1); ctr<this->n_systems(); ctr++)

 #endif // #define __variable_h__
Index: src/systems/system.C
===================================================================
--- src/systems/system.C	(revision 6517)
+++ src/systems/system.C	(working copy)
@@ -1073,53 +1073,86 @@
 			           const FEType& type,
