From 7d2e3815da869dc567a67e10349230e9d79e82a8 Mon Sep 17 00:00:00 2001
From: Oliver Sander <oliver.sander@tu-dresden.de>
Date: Thu, 22 Nov 2018 15:47:24 +0100
Subject: [PATCH] GlobalP2Mapper: Make sure the grid view contains only
quadrilaterals
The current implementation assign one degree of freedom to each
element, regardless of its type. That only works for quad grids:
We are implementing Lagrange FE here, and these don't assign
degrees of freedom to triangles.
---
dune/gfe/parallel/globalp2mapper.hh | 3 +++
1 file changed, 3 insertions(+)
diff --git a/dune/gfe/parallel/globalp2mapper.hh b/dune/gfe/parallel/globalp2mapper.hh
index 55c98805..1bb6a3ed 100644
--- a/dune/gfe/parallel/globalp2mapper.hh
+++ b/dune/gfe/parallel/globalp2mapper.hh
@@ -33,6 +33,9 @@ namespace Dune {
{
static_assert(GridView::dimension==2, "Only implemented for two-dimensional grids");
+ if (gridView.size(GeometryTypes::triangle)>1)
+ DUNE_THROW(NotImplemented, "GlobalP2Mapper only works for quad grids!");
+
GlobalIndexSet<GridView> globalVertexIndex(gridView,2);
GlobalIndexSet<GridView> globalEdgeIndex(gridView,1);
GlobalIndexSet<GridView> globalElementIndex(gridView,0);
--
GitLab