diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc index cc5ece24d241498c30ef7c021212c79593569985..5f105be4b4947495301d46945cc6908dd0e357c5 100644 --- a/src/one-body-sample.cc +++ b/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