Skip to content
Snippets Groups Projects
Commit b99a85d1 authored by Markus Wolff's avatar Markus Wolff
Browse files

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
parent 6a8365e1
No related branches found
No related tags found
No related merge requests found
......@@ -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_)
......
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