Skip to content
Snippets Groups Projects
Commit 4a97dcf2 authored by Elias Pipping's avatar Elias Pipping
Browse files

[Extern] Adjust to changes in the MCMG mapper

parent e703313e
No related branches found
No related tags found
No related merge requests found
...@@ -37,7 +37,7 @@ class GlobalRateStateFriction : public GlobalFriction<Matrix, Vector> { ...@@ -37,7 +37,7 @@ class GlobalRateStateFriction : public GlobalFriction<Matrix, Vector> {
GridView, Dune::MCMGVertexLayout> const vertexMapper(gridView); GridView, Dune::MCMGVertexLayout> const vertexMapper(gridView);
for (auto it = gridView.template begin<block_size>(); for (auto it = gridView.template begin<block_size>();
it != gridView.template end<block_size>(); ++it) { it != gridView.template end<block_size>(); ++it) {
auto const i = vertexMapper.map(*it); auto const i = vertexMapper.index(*it);
auto const coordinate = it->geometry().corner(0); auto const coordinate = it->geometry().corner(0);
if (not frictionalBoundary.containsVertex(i)) { if (not frictionalBoundary.containsVertex(i)) {
restrictions[i] = trivialNonlinearity; restrictions[i] = trivialNonlinearity;
......
...@@ -297,7 +297,7 @@ int main(int argc, char *argv[]) { ...@@ -297,7 +297,7 @@ int main(int argc, char *argv[]) {
for (auto it = leafView.begin<dims>(); it != leafView.end<dims>(); ++it) { for (auto it = leafView.begin<dims>(); it != leafView.end<dims>(); ++it) {
auto const geometry = it->geometry(); auto const geometry = it->geometry();
assert(geometry.corners() == 1); assert(geometry.corners() == 1);
vertexCoordinates[vertexMapper.map(*it)] = geometry.corner(0); vertexCoordinates[vertexMapper.index(*it)] = geometry.corner(0);
} }
} }
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment