constraints(): simplify

src/constraints.hh

@@ -42,8 +42,7 @@ namespace PPS {
   }
 
   template<typename Pred, typename GFS, typename CG>
-  void constraints(const Pred& isDirichlet, const GFS& gfs, CG& cg,
-                   const bool verbose = false)
+  void constraints(const Pred& isDirichlet, const GFS& gfs, CG& cg)
   {
     // clear global constraints
     cg.clear();
@@ -67,18 +66,15 @@ namespace PPS {
       std::set<std::size_t> local_constraints;
 
       // iterate over intersections and call metaprogram
-      unsigned int intersection_index = 0;
       for (const auto& intersection : intersections(es,element))
       {
         if(intersection.boundary())
         {
           // iterate over boundary, need intersection iterator
-          if (lfs_e.size())
-            boundary(isDirichlet, intersection,
-                     lfs_e.finiteElement().localCoefficients(),
-                     local_constraints);
+          boundary(isDirichlet, intersection,
+                   lfs_e.finiteElement().localCoefficients(),
+                   local_constraints);
         }
-        ++intersection_index;
       }
 
       for(auto i : local_constraints)