From bbb0d78e557beade7ba06bae82ba86a141db55b2 Mon Sep 17 00:00:00 2001
From: Thomas Fetzer <thomas.fetzer@iws.uni-stuttgart.de>
Date: Fri, 8 Jan 2016 09:23:11 +0100
Subject: [PATCH] [implicit] Correct the calculation of the cornerpoints of
grid to work for all dim.
---
dumux/implicit/common/implicitproblem.hh | 11 +++++++----
1 file changed, 7 insertions(+), 4 deletions(-)
diff --git a/dumux/implicit/common/implicitproblem.hh b/dumux/implicit/common/implicitproblem.hh
index 6e74bc2cd1..2bf79b7304 100644
--- a/dumux/implicit/common/implicitproblem.hh
+++ b/dumux/implicit/common/implicitproblem.hh
@@ -689,10 +689,13 @@ public:
*/
bool isCornerPoint(const GlobalPosition &globalPos, Scalar eps = 1e-8)
{
- return ((globalPos[0] < asImp_().bBoxMin()[0] + eps && globalPos[1] < asImp_().bBoxMin()[1] + eps)
- || (globalPos[0] < asImp_().bBoxMin()[0] + eps && globalPos[1] > asImp_().bBoxMax()[1] - eps)
- || (globalPos[0] > asImp_().bBoxMax()[0] - eps && globalPos[1] < asImp_().bBoxMin()[1] + eps)
- || (globalPos[0] > asImp_().bBoxMax()[0] - eps && globalPos[1] > asImp_().bBoxMax()[1] - eps));
+ for (unsigned int dimIdx = 0; dimIdx < dimWorld; dimIdx++)
+ {
+ if (!(globalPos[dimIdx] < asImp_().bBoxMin()[dimIdx] + eps
+ || globalPos[dimIdx] > asImp_().bBoxMax()[dimIdx] - eps))
+ return false;
+ }
+ return true;
}
/*!
--
GitLab