Skip to content
Snippets Groups Projects
Commit d4771bb1 authored by Dennis Gläser's avatar Dennis Gläser Committed by Timo Koch
Browse files

[facet][mpfa][mapper] make sure stencils are unique

parent fadf705c
No related branches found
No related tags found
1 merge request!1431Feature/facetcoupling with mpfa
......@@ -210,14 +210,11 @@ public:
auto& bulkCouplingData = this->couplingMap_(bulkGridId, facetGridId);
std::for_each(bulkCouplingData.begin(), bulkCouplingData.end(), makeBulkMapEntryUnique);
// coupling stencils might not be unique if box is used in lowdim domain
if (LowDimFVG::discMethod == DiscretizationMethod::box)
{
auto& bulkCouplingData = this->couplingMap_(bulkGridId, facetGridId);
std::for_each(bulkCouplingData.begin(),
bulkCouplingData.end(),
[&makeUnique] (auto& pair) { makeUnique(pair.second.couplingStencil); });
}
// coupling stencils might not be unique in lowdim domain
auto& lowDimCouplingData = this->couplingMap_(facetGridId, bulkGridId);
std::for_each(lowDimCouplingData.begin(),
lowDimCouplingData.end(),
[&makeUnique] (auto& pair) { makeUnique(pair.second.couplingStencil); });
}
private:
......
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