#2534: Evolution of grid arrays with Mol no longer works
| Reporter: | |
| Status: | new |
| Milestone: | |
| Version: | development version |
| Type: | bug |
| Priority: | major |
| Component: | Carpet |
Comment (by Yosef Zlochower):
The following patch seems to “fix” the problem.
diff --git a/Carpet/src/Poison.cc b/Carpet/src/Poison.cc index 12fc3d3..39a989a 100644 --- a/Carpet/src/Poison.cc +++ b/Carpet/src/Poison.cc @@ -131,12 +131,16 @@ void PoisonGroup(cGH const *const cctkGH, int const group, END_LOCAL_MAP_LOOP; // mark as invalid - for (int var = 0; var < nvar; ++var) { - int const map0 = 0; - ggf *const ff = arrdata.AT(group).AT(map0).data.AT(var); - assert(ff); - for (int tl = min_tl; tl <= max_tl; ++tl) { - ff->set_valid(mglevel, reflevel, tl, CCTK_VALID_NOWHERE); + // Exclude CCTK_ARRAY from this + if (CCTK_GroupTypeI(group) != CCTK_ARRAY) + { + for (int var = 0; var < nvar; ++var) { + int const map0 = 0; + ggf *const ff = arrdata.AT(group).AT(map0).data.AT(var); + assert(ff); + for (int tl = min_tl; tl <= max_tl; ++tl) { + ff->set_valid(mglevel, reflevel, tl, CCTK_VALID_NOWHERE); + } } }
--
Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2534/evolution-of-grid-arrays-with-mol-no