From e1d43a6d594dd1994ab8127469c3f8c3da54c03b Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Wed, 12 Dec 2012 15:42:15 +0100
Subject: [PATCH] Record coordinates
---
src/one-body-sample.cc | 6 ++++--
1 file changed, 4 insertions(+), 2 deletions(-)
diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc
index 3071ba35..ff9568ff 100644
--- a/src/one-body-sample.cc
+++ b/src/one-body-sample.cc
@@ -192,15 +192,17 @@ int main(int argc, char *argv[]) {
Dune::BitSetVector<dims> ignoreNodes(finestSize, false);
Dune::BitSetVector<1> neumannNodes(finestSize, false);
Dune::BitSetVector<1> frictionalNodes(finestSize, false);
+
+ VectorType coordinates(finestSize);
{
Dune::MultipleCodimMultipleGeomTypeMapper<
GridView, Dune::MCMGVertexLayout> const myVertexMapper(leafView);
for (auto it = leafView.begin<dims>(); it != leafView.end<dims>(); ++it) {
assert(it->geometry().corners() == 1);
- Dune::FieldVector<double, dims> const localCoordinates =
- it->geometry().corner(0);
size_t const id = myVertexMapper.map(*it);
+ coordinates[id] = it->geometry().corner(0);
+ auto const &localCoordinates = coordinates[id];
// Find the center of the lower face
switch (dims) {
--
GitLab