Inform the Accelerator thorn that past timelevels have been filled on the host

2 files changed, 13 insertions, 0 deletions

diff --git a/Carpet/Carpet/interface.ccl b/Carpet/Carpet/interface.ccl
index 506f9e8d4..4f53054f6 100644
--- a/Carpet/Carpet/interface.ccl
+++ b/Carpet/Carpet/interface.ccl
@@ -231,6 +231,12 @@ CCTK_INT FUNCTION Accelerator_PostSync          \
    CCTK_INT              IN ngroups)
 USES FUNCTION Accelerator_PostSync
 
+CCTK_INT FUNCTION Accelerator_NotifyVariableWritten       \
+  (CCTK_POINTER_TO_CONST IN cctkGH,                       \
+   CCTK_INT              IN vi,                           \
+   CCTK_INT              IN tl,                           \
+   CCTK_INT              IN onhost)
+USES FUNCTION Accelerator_NotifyVariableWritten
 
 
 # Access coordinate information (on the coarse level)
diff --git a/Carpet/Carpet/src/Cycle.cc b/Carpet/Carpet/src/Cycle.cc
index bd0acac0a..64b7ff25b 100644
--- a/Carpet/Carpet/src/Cycle.cc
+++ b/Carpet/Carpet/src/Cycle.cc
@@ -238,6 +238,13 @@ namespace Carpet {
           assert (reflevel>=0 and reflevel<reflevels);
 	  for (int m=0; m<(int)arrdata.AT(group).size(); ++m) {
 	    for (int var=0; var<CCTK_NumVarsInGroupI(group); ++var) {
+
+              if (CCTK_IsFunctionAliased("Accelerator_NotifyVariableWritten")) {
+                for (int tl = 1; tl < arrdata.AT(group).AT(m).data.AT(var)->timelevels(mglevel,reflevel); tl++) {
+                  Accelerator_NotifyVariableWritten(cctkGH, CCTK_FirstVarIndexI(group)+var, tl, 1 /* on host */);
+                }
+              }
+
               arrdata.AT(group).AT(m).data.AT(var)->
                 fill_all (reflevel, mglevel);
             }