this->_dof_coupling,
     this->_coupling_functors,
     implicit_neighbor_dofs,
     _need_full_sparsity_pattern,
     calculate_constrained,
     sc);

  _extra_send_list_context(nullptr),
  _default_coupling(std::make_unique<DefaultCoupling>()),
  _default_evaluating(std::make_unique<DefaultCoupling>()),
  _need_full_sparsity_pattern(false),
  _need_ghost_constraints(false),
  _n_SCALAR_dofs(0)
#ifdef LIBMESH_ENABLE_AMR
  , _first_old_scalar_df()

  this->update_sparsity_pattern(matrix);

  if (matrix.need_full_sparsity_pattern())
    _need_full_sparsity_pattern = true;
}

        {
          // We'd better have already computed the full sparsity
          // pattern if we need it here
          libmesh_assert(_need_full_sparsity_pattern);

          matrix.update_sparsity_pattern (_sp->get_sparsity_pattern());
        }

  _first_scalar_df.clear();
  this->clear_send_list();
  this->clear_sparsity();
  _need_full_sparsity_pattern = false;
  _need_ghost_constraints = false;

#ifdef LIBMESH_ENABLE_AMR

  for (const auto & mat : _matrices)
    {
      mat->attach_sparsity_pattern (*_sp);
      if (_need_full_sparsity_pattern)
        mat->update_sparsity_pattern (_sp->get_sparsity_pattern());
    }
  // If we don't need the full sparsity pattern anymore, free the
  // parts of it we don't need.
  if (!_need_full_sparsity_pattern)
    _sp->clear_full_sparsity();
}

  // Now that we have our root constraint dependencies sorted out, add
  // them to the send_list
  this->add_constraints_to_send_list(mesh);
}

    }
}

void DofMap::add_constraints_to_send_list (const MeshBase & mesh)
{
  // This function must be run on all processors at once
  parallel_object_only();

  // If we only need constraint DoFs constraining DoFs which are
  // algebraically local, we're done.
  if (!this->has_static_condensation() && !_need_ghost_constraints)
    return;

  // If we use StaticCondensation, though, we may need constraint DoFs

  //
  // We'll potentially be hitting the same constrained DoFs from
  // multiple directions.
  for (auto & elem : mesh.active_local_element_ptr_range())
    {
      std::vector<dof_id_type> di;
      this->dof_indices (elem, di);