[boxdfm] store and expose indexInIntersection in scv(f)

dumux/porousmediumflow/boxdfm/subcontrolvolume.hh

@@ -93,6 +93,7 @@ public:
     , elemLocalScvIdx_(scvIdx)
     , dofIndex_(dofIndex)
     , facetIdx_(0)
+    , indexInIntersection_(0)
     {
         auto corners = geometryHelper.getScvCorners(scvIdx);
         dofPosition_ = corners[0];
@@ -132,6 +133,7 @@ public:
     , elemLocalScvIdx_(elemLocalScvIdx)
     , dofIndex_(dofIndex)
     , facetIdx_(elemLocalFacetIdx)
+    , indexInIntersection_(indexInIntersection)
     {
         auto corners = geometryHelper.getBoundaryScvfCorners(intersection.indexInInside(), indexInIntersection);
         dofPosition_ = corners[0];
@@ -166,6 +168,10 @@ public:
     LocalIndexType facetIndexInElement() const
     { assert(isFractureScv_); return facetIdx_; }
 
+    //! The local vertex index in the intersection
+    LocalIndexType indexInsideIntersection() const
+    { assert(isFractureScv_); return indexInIntersection_; }
+
     //! The index of the dof this scv is embedded in
     GridIndexType dofIndex() const
     { return dofIndex_; }
@@ -194,6 +200,7 @@ private:
 
     // for fracture scvs only!
     LocalIndexType facetIdx_;
+    LocalIndexType indexInIntersection_;
 };
 
 } // end namespace

dumux/porousmediumflow/boxdfm/subcontrolvolumeface.hh

@@ -95,6 +95,7 @@ public:
     , isFractureScvf_(false)
     , boundaryFlag_{}
     , facetIdx_(0)
+    , indexInIntersection_(0)
     {
         auto corners = geometryHelper.getScvfCorners(scvfIndex);
         unitOuterNormal_ = geometryHelper.normal(corners, scvIndices_);
@@ -123,6 +124,7 @@ public:
     , isFractureScvf_(false)
     , boundaryFlag_{intersection}
     , facetIdx_(0)
+    , indexInIntersection_(0)
     {
         auto corners = geometryHelper.getBoundaryScvfCorners(intersection.indexInInside(), indexInIntersection);
         area_ = Dumux::convexPolytopeVolume<T::dim-1>(
@@ -149,6 +151,7 @@ public:
     , isFractureScvf_(true)
     , boundaryFlag_{intersection}
     , facetIdx_(intersection.indexInInside())
+    , indexInIntersection_(indexInIntersection)
     {
         auto corners = geometryHelper.getFractureScvfCorners(intersection.indexInInside(), indexInIntersection);
         // The area here is given in meters. In order to
@@ -200,6 +203,10 @@ public:
     LocalIndexType facetIndexInElement() const
     { assert(isFractureScvf_); return facetIdx_; }
 
+    //! The local edge index inside the intersection
+    LocalIndexType indexInIntersection() const
+    { assert(isFractureScvf_); return indexInIntersection_; }
+
     //! Returns the boundary flag
     typename BoundaryFlag::value_type boundaryFlag() const
     { return boundaryFlag_.get(); }
@@ -232,6 +239,7 @@ private:
     bool isFractureScvf_;
     BoundaryFlag boundaryFlag_;
     LocalIndexType facetIdx_;
+    LocalIndexType indexInIntersection_;
 };
 
 } // end namespace Dumux