From b99a85d12bd7d15dbc336228544437d1f5e9facf Mon Sep 17 00:00:00 2001
From: Markus Wolff <markus.wolff@twt-gmbh.de>
Date: Tue, 6 Sep 2011 12:45:02 +0000
Subject: [PATCH] switched back the way of getting the index of an intersection
in the discretization of the 2p saturation equation
- way of getting the index of an intersection was changed for
use of adaptive grids. However, this only works in the case in which the order of the iteration by
the IntersectionIterator coincides with the ordering at the
reference element (e.g. not valid for UGGrid)
- Switch back to use indexInInside() function -> Adaptive method will
not work properly until a better solution is found!
git-svn-id: svn://svn.iws.uni-stuttgart.de/DUMUX/dumux/trunk@6589 2fb0f335-1f38-0410-981e-8018bf24f1b0
---
dumux/decoupled/2p/transport/fv/fvsaturation2p.hh | 6 +++---
1 file changed, 3 insertions(+), 3 deletions(-)
diff --git a/dumux/decoupled/2p/transport/fv/fvsaturation2p.hh b/dumux/decoupled/2p/transport/fv/fvsaturation2p.hh
index ad21f09f34..c82b019dff 100644
--- a/dumux/decoupled/2p/transport/fv/fvsaturation2p.hh
+++ b/dumux/decoupled/2p/transport/fv/fvsaturation2p.hh
@@ -340,14 +340,14 @@ int FVSaturation2P<TypeTag>::update(const Scalar t, Scalar& dt,
Scalar lambdaNWI = problem_.variables().mobilityNonwetting(globalIdxI);
// run through all intersections with neighbors and boundary
- int isIndex = -1;
+// int isIndex = -1;
IntersectionIterator isItEnd = problem_.gridView().iend(*eIt);
for (IntersectionIterator isIt = problem_.gridView().ibegin(*eIt); isIt
!= isItEnd; ++isIt)
{
// local number of facet
- isIndex++;
-// int isIndex = isIt->indexInInside();
+// isIndex++;
+ int isIndex = isIt->indexInInside();
GlobalPosition unitOuterNormal = isIt->centerUnitOuterNormal();
if (switchNormals_)
--
GitLab