From b74dbfba10dba09c3049fa66c6cb53af189da8c9 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Maximilian=20H=C3=B6rl?=
<maximilian.hoerl@mathematik.uni-stuttgart.de>
Date: Thu, 7 May 2020 17:45:09 +0200
Subject: [PATCH] [bugfix] fix subentity codimension in getDiameter()
---
dune/mmdg/dg.hh | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/dune/mmdg/dg.hh b/dune/mmdg/dg.hh
index c842b35..623c825 100644
--- a/dune/mmdg/dg.hh
+++ b/dune/mmdg/dg.hh
@@ -521,9 +521,9 @@ protected:
{
Scalar diameter(0.0);
- for (int i = 0; i < ((dim + 1) * dim) / 2; i++)
+ for (int i = 0; i < elem.subEntities(dim-1); i++)
{
- const auto& edge = elem.template subEntity<1>(i);
+ const auto& edge = elem.template subEntity<dim-1>(i);
diameter = std::max(diameter, edge.geometry().volume());
}
--
GitLab