Merge pull request #6706 from aritorto/onlyLevelZeroNNC

Skip refined cells when computing trans and nncs for output
This commit is contained in:
Markus Blatt authored and GitHub committed 2025-12-30 09:39:56 +01:00
commit 36266401a4
1 file changed
+6
@@ -344,6 +344,9 @@ computeTrans_(const std::unordered_map<int,int>& cartesianToActive,
if (!is.neighbor())
continue; // intersection is on the domain boundary
if ( (is.inside().level()>0) || (is.outside().level()>0))
continue; // for CpGrid with LGRs, we only care about level zero cells, for now.
// Not 'const' because remapped if 'map' is non-null.
unsigned c1 = globalElemMapper.index(is.inside());
unsigned c2 = globalElemMapper.index(is.outside());
@@ -490,6 +493,9 @@ exportNncStructure_(const std::unordered_map<int,int>& cartesianToActive,
if (!is.neighbor())
continue; // intersection is on the domain boundary
if ( (is.inside().level()>0) || (is.outside().level()>0))
continue; // for CpGrid with LGRs, we only care about level zero cells, for now.
// Not 'const' because remapped if 'map' is non-null.
unsigned c1 = globalElemMapper.index(is.inside());
unsigned c2 = globalElemMapper.index(is.outside());