From c25dcdd8783a9d0e543c0fe97f9124dcebdee6ab Mon Sep 17 00:00:00 2001
From: Timo Koch <timo.koch@iws.uni-stuttgart.de>
Date: Mon, 4 Sep 2017 10:39:51 +0200
Subject: [PATCH] [test] Small readability improvement in 2d growth test
---
dune/foamgrid/test/growth-test-2d.cc | 10 +++++++---
1 file changed, 7 insertions(+), 3 deletions(-)
diff --git a/dune/foamgrid/test/growth-test-2d.cc b/dune/foamgrid/test/growth-test-2d.cc
index 5e91ec3..ac53f15 100644
--- a/dune/foamgrid/test/growth-test-2d.cc
+++ b/dune/foamgrid/test/growth-test-2d.cc
@@ -20,19 +20,21 @@ void checkGridElementGrowth(Grid& grid)
enum { dim = Grid::dimension };
// the vertex mapper
- Dune::LeafMultipleCodimMultipleGeomTypeMapper<Grid,Dune::MCMGVertexLayout> mapper(grid);
+ LeafMultipleCodimMultipleGeomTypeMapper<Grid, MCMGVertexLayout> mapper(grid);
// let the top element grow
for (const auto& element : elements(grid.leafGridView()))
+ {
+ const auto refElement = ReferenceElements<double,dim>::general(element.type());
for (const auto& intersection : intersections(grid.leafGridView(), element))
+ {
if(intersection.boundary() && intersection.geometry().center()[1] > 0.5)
{
- const auto refElement = ReferenceElements<double,dim>::general(element.type());
auto&& v0 = element.template subEntity<dim>(refElement.subEntity(intersection.indexInInside(), dim-1, 0, dim));
auto&& v1 = element.template subEntity<dim>(refElement.subEntity(intersection.indexInInside(), dim-1, 1, dim));
// calculate new vertex position
- Dune::FieldVector<double, dimworld> newVertex = v0.geometry().center();
+ auto newVertex = v0.geometry().center();
newVertex += v1.geometry().center(); newVertex /= 2.0;
newVertex[1] += 0.5;
@@ -42,6 +44,8 @@ void checkGridElementGrowth(Grid& grid)
// insert the element
grid.insertElement(element.type(), {mapper.index(v0), mapper.index(v1), vIdx});
}
+ }
+ }
std::size_t numBoundarySegments = grid.numBoundarySegments();
std::cout << std::endl<< "numBoundarySegments before growth: " << numBoundarySegments << std::endl;
--
GitLab