From 6026cf28c1bc3d2720a38f77bab892d67364a349 Mon Sep 17 00:00:00 2001
From: Liwei Ji <jiliwei.phys@gmail.com>
Date: Wed, 20 Nov 2024 11:08:30 -0500
Subject: [PATCH] ADMconstraints: add extra SYNCs

---
 ADMconstraints/schedule.ccl           | 14 ++++++++++++--
 ADMconstraints/src/admconstraints.cxx |  4 ++++
 2 files changed, 16 insertions(+), 2 deletions(-)

diff --git a/ADMconstraints/schedule.ccl b/ADMconstraints/schedule.ccl
index dacac6c2..5cff5d09 100644
--- a/ADMconstraints/schedule.ccl
+++ b/ADMconstraints/schedule.ccl
@@ -12,7 +12,17 @@ SCHEDULE GROUP ADMconstraints_AnalysisGroup AT analysis
 
 
 if (calc_constraints) {
-  SCHEDULE ADMconstraints_CalcConstraints IN ADMconstraints_AnalysisGroup
+  SCHEDULE ADMconstraints_Sync IN ADMconstraints_AnalysisGroup
+  {
+    LANG: C
+    OPTIONS: global
+    SYNC: ADMBaseX::metric
+    SYNC: ADMBaseX::curv
+    SYNC: ADMBaseX::lapse
+    SYNC: ADMBaseX::shift
+  } "Synchronize"
+
+  SCHEDULE ADMconstraints_CalcConstraints IN ADMconstraints_AnalysisGroup AFTER ADMconstraints_Sync
   {
     LANG: C
     READS: ADMBaseX::metric(everywhere)
@@ -24,7 +34,6 @@ if (calc_constraints) {
     READS: TmunuBaseX::eTij(interior)
     WRITES: HC(interior)
     WRITES: MC(interior)
-    SYNC: HC MC
   } "Calculate ADM Constraints"
 
   if (use_lapse_mask) {
@@ -36,6 +45,7 @@ if (calc_constraints) {
       READS: MC(interior)
       WRITES: HC(interior)
       WRITES: MC(interior)
+      SYNC: HC MC
     } "Apply lapse mask"
   }
 }
diff --git a/ADMconstraints/src/admconstraints.cxx b/ADMconstraints/src/admconstraints.cxx
index e51a7eb2..18abbc2f 100644
--- a/ADMconstraints/src/admconstraints.cxx
+++ b/ADMconstraints/src/admconstraints.cxx
@@ -29,6 +29,10 @@ CCTK_DEVICE CCTK_HOST CCTK_ATTRIBUTE_ALWAYS_INLINE inline T Power(T x, int y) {
   return (y == 2) ? Arith::pow2(x) : Arith::pown(x, y);
 }
 
+extern "C" void ADMconstraints_Sync(CCTK_ARGUMENTS) {
+  // do nothing
+}
+
 extern "C" void ADMconstraints_CalcConstraints(CCTK_ARGUMENTS) {
   DECLARE_CCTK_ARGUMENTSX_ADMconstraints_CalcConstraints;
   DECLARE_CCTK_PARAMETERS;