Commit c8be3d69 authored by Kilian Weishaupt's avatar Kilian Weishaupt
Browse files

[discretization][boxfvgeometry] Construct correct number of scv faces in 3d

* number of inner scvfs is element.subEntities(dim-1), not element.subEntities(1)
parent 005fd4a6
......@@ -299,7 +299,7 @@ private:
}
// construct the sub control volume faces
const auto numInnerScvf = (dim==1) ? 1 : element.subEntities(1);
const auto numInnerScvf = (dim==1) ? 1 : element.subEntities(dim-1);
scvfs_.resize(numInnerScvf);
unsigned int scvfLocalIdx = 0;
......
......@@ -141,8 +141,8 @@ public:
// construct the sub control volume faces
unsigned int scvfLocalIdx = 0;
scvfs_[eIdx].resize(element.subEntities(1));
for (; scvfLocalIdx < element.subEntities(1); ++scvfLocalIdx)
scvfs_[eIdx].resize(element.subEntities(dim-1));
for (; scvfLocalIdx < element.subEntities(dim-1); ++scvfLocalIdx)
{
// find the global and local scv indices this scvf is belonging to
std::vector<IndexType> localScvIndices({static_cast<IndexType>(referenceElement.subEntity(scvfLocalIdx, dim-1, 0, dim)),
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment