Skip to content
Snippets Groups Projects
Commit 62cd7529 authored by akbib's avatar akbib Committed by akbib@FU-BERLIN.DE
Browse files

same bugfix as in revision 352, this time for the case of a nonconforming hostgrid

[[Imported from SVN: r353]]
parent f74b6f0b
No related branches found
No related tags found
No related merge requests found
......@@ -133,10 +133,18 @@ public:
{
// top to bottom
FatherFaceIterator fatherIt(subGrid_,hostIterator_,FatherFaceIterator::begin);
FatherFaceIterator fatherEndIt(subGrid_,hostIterator_,FatherFaceIterator::end);
// get first internal face
while (true) {
if (fatherIt->inside()->level() == 0 || (fatherIt->neighbor() && subGrid_->template contains<0>(*(fatherIt->outside()))))
break;
else if (fatherIt == fatherEndIt) {
// if no father found the intersection is in the inside of an element that is not completly covered by the subgrid children,
// This can happen when "insertPartial" is used. In this case this is a boundary intersection
hierarchicFaceIterator_ = HierarchicFaceIterator(subGrid_, hostIterator_, hostIterator_->inside(), 12345);
return;
}
else
++fatherIt;
}
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment