Skip to content
Snippets Groups Projects
Commit d326aa0c authored by alkaemper's avatar alkaemper
Browse files

avoid checking the artificial faces for periodic boundaries

parent d23ed7af
No related branches found
No related tags found
No related merge requests found
...@@ -975,6 +975,17 @@ namespace Dune ...@@ -975,6 +975,17 @@ namespace Dune
const FaceMapIterator fend = faceMap.end(); const FaceMapIterator fend = faceMap.end();
for( FaceMapIterator fit = faceMap.begin(); fit != fend; ++fit ) for( FaceMapIterator fit = faceMap.begin(); fit != fend; ++fit )
{ {
//for dimension == 2 we do not want to search
// the artificially introduced faces
if(dimension == 2)
{
if(elementType == hexa)
if(fit->second.second > 3)
continue;
if(elementType == tetra)
if(fit->second.second > 2)
continue;
}
FaceType key2; FaceType key2;
generateFace( fit->second, key2 ); generateFace( fit->second, key2 );
...@@ -1064,6 +1075,17 @@ namespace Dune ...@@ -1064,6 +1075,17 @@ namespace Dune
{ {
for(auto it = faceMap.begin(); it!=faceMap.end(); ++it) for(auto it = faceMap.begin(); it!=faceMap.end(); ++it)
{ {
//for dimension == 2 we do not want to search
// the artificially introduced faces
if(dimension == 2)
{
if(elementType == hexa)
if(it->second.second > 3)
continue;
if(elementType == tetra)
if(it->second.second > 2)
continue;
}
searchPeriodicNeighbor( faceMap, it, defaultId ); searchPeriodicNeighbor( faceMap, it, defaultId );
} }
} }
......
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