From 3c829031f02f83235878cc52bb5304187ec0dc57 Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Sat, 14 Dec 2013 13:37:26 +0100
Subject: [PATCH] [Cleanup] Move vertexCoordinates
---
src/one-body-sample.cc | 22 +++++++++++-----------
1 file changed, 11 insertions(+), 11 deletions(-)
diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc
index 2d3a781a..7d356f56 100644
--- a/src/one-body-sample.cc
+++ b/src/one-body-sample.cc
@@ -154,17 +154,6 @@ int main(int argc, char *argv[]) {
});
}
- Vector vertexCoordinates(fineVertexCount);
- {
- Dune::MultipleCodimMultipleGeomTypeMapper<
- GridView, Dune::MCMGVertexLayout> const vertexMapper(leafView);
- for (auto it = leafView.begin<dims>(); it != leafView.end<dims>(); ++it) {
- auto const geometry = it->geometry();
- assert(geometry.corners() == 1);
- vertexCoordinates[vertexMapper.map(*it)] = geometry.corner(0);
- }
- }
-
// Neumann boundary
BoundaryPatch<GridView> const neumannBoundary(leafView);
@@ -340,6 +329,17 @@ int main(int argc, char *argv[]) {
auto multigridStep = factory.getSolver();
{
+ Vector vertexCoordinates(fineVertexCount);
+ {
+ Dune::MultipleCodimMultipleGeomTypeMapper<
+ GridView, Dune::MCMGVertexLayout> const vertexMapper(leafView);
+ for (auto it = leafView.begin<dims>(); it != leafView.end<dims>();
+ ++it) {
+ auto const geometry = it->geometry();
+ assert(geometry.corners() == 1);
+ vertexCoordinates[vertexMapper.map(*it)] = geometry.corner(0);
+ }
+ }
std::fstream vertexCoordinateWriter("coordinates", std::fstream::out);
for (size_t i = 0; i < fineVertexCount; ++i)
if (frictionalNodes[i][0])
--
GitLab