PABLO: mark octants as "new for refinement" only if they are generate…

…d by a refinement
optimad · Jun 12, 2024 · d2dbdaf · d2dbdaf
1 parent 0bd8c75
commit d2dbdaf
Show file tree

Hide file tree

Showing 2 changed files with 5 additions and 2 deletions.
diff --git a/src/PABLO/LocalTree.cpp b/src/PABLO/LocalTree.cpp
@@ -524,6 +524,11 @@ namespace bitpit {
                 uint32_t firstChildIdx = futureIdx - (nChildren - 1);
                 fatherOctant.buildChildren(m_octants.data() + firstChildIdx);
 
+                // Set children information
+                for (int i = 0; i < nChildren; ++i) {
+                    m_octants[firstChildIdx + i].m_info[Octant::INFO_NEW4REFINEMENT] = true;
+                }
+
                 // Update the mapping
                 if(!mapidx.empty()){
                     for (uint8_t i=0; i<nChildren; i++){

diff --git a/src/PABLO/Octant.cpp b/src/PABLO/Octant.cpp
@@ -844,8 +844,6 @@ void	Octant::buildChildren(Octant *children) const {
 		uint32_t dh = oct.getLogicalSize();
 		oct.m_morton = PABLO::computeMorton(m_dim, coords[0] + dh * dx, coords[1] + dh * dy, coords[2] + dh * dz);
 
-		oct.m_info[OctantInfo::INFO_NEW4REFINEMENT] = true;
-
 		oct.m_info[INFO_BOUNDFACE0 + xf] = false;
 		oct.m_info[INFO_BOUNDFACE0 + yf] = false;
 		oct.m_info[INFO_BOUNDFACE0 + zf] = false;