feature partitioning: consider newly_refined patch at end of skip_local

hannesbrandt · Aug 8, 2024 · f63af68 · f63af68
1 parent c15a2a0
commit f63af68
Showing 1 changed file with 6 additions and 1 deletion.
diff --git a/src/forestclaw2d.c b/src/forestclaw2d.c
@@ -1798,9 +1798,14 @@ fclaw2d_domain_iterate_pack (fclaw2d_domain_t * domain, size_t data_size,
                 {
                     break;      /* jumped to end of array */
                 }
+                /* skip newly refined patches that stay local */
                 while (skip_refined && next_refined < i)
                 {
-                    /* skip newly refined patches that stay local */
+                    if (next_refined + P4EST_CHILDREN > i) {
+                        /* i is a newly refined sibling patch */
+                        num_src -= next_refined + P4EST_CHILDREN - i;
+                        i = next_refined + P4EST_CHILDREN;
+                    }
                     next_refined = (nri == num_nr) ? old_lnp :
                         *(int *) sc_array_index_int (wrap->newly_refined,
                                                      nri++);