Use face hierarchy instead of intersection iterator

dune/fufem/functionspacebases/dofconstraints.hh

@@ -12,6 +12,7 @@
 
 
 #include "dune/fufem/functions/alienelementlocalbasisfunction.hh"
+//#include "dune/fufem/facehierarchy.hh"
 
 /** \brief Class to handle linear constraints on a subset of a set of DOFs
  *
@@ -128,6 +129,9 @@ class DOFConstraints {
                         continue;
 
                     NeighborPointer np = nIt->outside();
+//                for(int k=0; k<it->template count<1>(); ++k)
+//                {
+//                    NeighborPointer np = FaceHierarchy<typename B::GridView::Grid>::finestCoarseNeighbor(gridview.grid(), *it, k);
                     auto& neighbor = *np;
 
                     if (not(neighbor.level()<level))