diff --git a/test/discretization/box/test_boxfvgeometry.cc b/test/discretization/box/test_boxfvgeometry.cc index 227191bc55f96971dc174bc99da585526df45b2f..b2130b25dacfd7e010c31bf630bb5a9e47ebba95 100644 --- a/test/discretization/box/test_boxfvgeometry.cc +++ b/test/discretization/box/test_boxfvgeometry.cc @@ -98,12 +98,21 @@ int main (int argc, char *argv[]) try if(0 != testForwardIterator(range2.begin(), range2.end(), op2)) DUNE_THROW(Dune::Exception, "Iterator does not fulfill the forward iterator concept"); + std::size_t boundaryCount = 0; for (auto&& scvf : scvfs(fvGeometry)) { std::cout << "-- scvf " << scvf.index() << " ip at: " << scvf.ipGlobal() << " normal: " << scvf.unitOuterNormal(); - if (scvf.boundary()) std::cout << " (on boundary)."; + if (scvf.boundary()) + { + ++boundaryCount; + std::cout << " (on boundary)."; + } std::cout << std::endl; } + + if ((boundaryCount>0) != fvGeometry.hasBoundaryScvf()) + DUNE_THROW(Dune::InvalidStateException, "fvGeometry.hasBoundaryScvf() reports " << fvGeometry.hasBoundaryScvf() + << " but the number of boundary scvfs is " << boundaryCount); } } // ////////////////////////////////// diff --git a/test/discretization/cellcentered/tpfa/test_tpfafvgeometry.cc b/test/discretization/cellcentered/tpfa/test_tpfafvgeometry.cc index e7223e64bbef2895fc380130cb67e121f565ba7e..49d715f3ff45ce05bd0c132f6695d3a6fc85fa52 100644 --- a/test/discretization/cellcentered/tpfa/test_tpfafvgeometry.cc +++ b/test/discretization/cellcentered/tpfa/test_tpfafvgeometry.cc @@ -99,12 +99,21 @@ int main (int argc, char *argv[]) try if(0 != testForwardIterator(range2.begin(), range2.end(), op2)) DUNE_THROW(Dune::Exception, "Iterator does not fulfill the forward iterator concept"); + std::size_t boundaryCount = 0; for (auto&& scvf : scvfs(fvGeometry)) { std::cout << "-- scvf " << scvf.index() << " ip at: " << scvf.ipGlobal(); - if (scvf.boundary()) std::cout << " (on boundary)."; + if (scvf.boundary()) + { + ++boundaryCount; + std::cout << " (on boundary)."; + } std::cout << std::endl; } + + if ((boundaryCount>0) != fvGeometry.hasBoundaryScvf()) + DUNE_THROW(Dune::InvalidStateException, "fvGeometry.hasBoundaryScvf() reports " << fvGeometry.hasBoundaryScvf() + << " but the number of boundary scvfs is " << boundaryCount); } } // ////////////////////////////////// diff --git a/test/discretization/staggered/test_staggeredfvgeometry.cc b/test/discretization/staggered/test_staggeredfvgeometry.cc index d3d45349e8bac2efb59c80a166d11b4253ac3793..130e0be2df6cc3d2028dfeaaaf67b4e1c9427413 100644 --- a/test/discretization/staggered/test_staggeredfvgeometry.cc +++ b/test/discretization/staggered/test_staggeredfvgeometry.cc @@ -133,12 +133,21 @@ int main (int argc, char *argv[]) try if(0 != testForwardIterator(range2.begin(), range2.end(), op2)) DUNE_THROW(Dune::Exception, "Iterator does not fulfill the forward iterator concept"); + std::size_t boundaryCount = 0; for (auto&& scvf : scvfs(fvGeometry)) { std::cout << "-- scvf " << scvf.index() << " ip at: " << scvf.ipGlobal() << " normal: " << scvf.unitOuterNormal(); - if (scvf.boundary()) std::cout << " (on boundary)."; + if (scvf.boundary()) + { + ++boundaryCount; + std::cout << " (on boundary)."; + } std::cout << std::endl; } + + if ((boundaryCount>0) != fvGeometry.hasBoundaryScvf()) + DUNE_THROW(Dune::InvalidStateException, "fvGeometry.hasBoundaryScvf() reports " << fvGeometry.hasBoundaryScvf() + << " but the number of boundary scvfs is " << boundaryCount); } } // //////////////////////////////////