[ET Trac] #2534: Evolution of grid arrays with Mol no longer works

[Yosef Zlochower]

#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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.einsteintoolkit.org/pipermail/trac/attachments/20210527/0227364b/attachment.html