Count nodes

src/one-body-sample.cc

@@ -105,6 +105,7 @@ int main() {
       size_t neumann_nodes = 0;
       size_t bounding_nodes = 0;
       size_t extremal_nodes = 0;
+      size_t overall_nodes = 0;
       for (VertexLeafIterator it = leafView.begin<dim>();
            it != leafView.end<dim>(); ++it) {
         Dune::GeometryType gt = it->type();
@@ -112,6 +113,7 @@ int main() {
         SmallVector coordinates = it->geometry().corner(0);
         bool bounding(false);
         bool extremal(true);
+        ++overall_nodes;
         for (int i = 0; i < dim; ++i) {
           if (coordinates[i] == end_points[i] || coordinates[i] == 0) {
             bounding = true;
@@ -152,6 +154,7 @@ int main() {
                 << std::endl;
       std::cout << "Number of bounding nodes: " << bounding_nodes << std::endl;
       std::cout << "Number of extremal nodes: " << extremal_nodes << std::endl;
+      std::cout << "Number of nodes: " << overall_nodes << std::endl;
     }
 
     // TODO: populate on S_F