Register symmetries for grid functions.

2 files changed, 34 insertions, 0 deletions

diff --git a/schedule.ccl b/schedule.ccl
index 244300f..8072632 100644
--- a/schedule.ccl
+++ b/schedule.ccl
@@ -40,6 +40,17 @@ if (CCTK_EQUALS(lapse_source, "QMS_MG")) {
         LANG: C
     } ""
 
+    SCHEDULE qms_mg_register_symmetries in SymmetryRegister {
+        LANG: C
+    } "register symmetry properties"
+
+    schedule qms_mg_select_bc in MoL_PostStep
+    {
+      LANG: C
+      OPTIONS: level
+      SYNC: W_val
+    } "select boundary conditions"
+
     STORAGE: W_pred0
     STORAGE: W_pred1
     STORAGE: W_pred0_time
diff --git a/src/qms.c b/src/qms.c
index 63fbce1..047e4b6 100644
--- a/src/qms.c
+++ b/src/qms.c
@@ -1691,3 +1691,26 @@ void qms_mg_terminate_print_stats(CCTK_ARGUMENTS)
         context_free(&qms_context);
     }
 }
+
+void qms_mg_register_symmetries(CCTK_ARGUMENTS)
+{
+  DECLARE_CCTK_ARGUMENTS;
+  DECLARE_CCTK_PARAMETERS;
+
+  int sym[3] = { 1, 1, 1 };
+  int ret;
+
+  ret = SetCartSymVN(cctkGH, sym, "QuasiMaximalSlicingMG::W_val");
+  if (ret != 0)
+      CCTK_WARN(0, "Error registering symmetries");
+}
+
+void qms_mg_select_bc(CCTK_ARGUMENTS)
+{
+  DECLARE_CCTK_ARGUMENTS;
+  DECLARE_CCTK_PARAMETERS;
+
+  int ret = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1, "QuasiMaximalSlicingMG::W_val", "none");
+  if (ret != 0)
+      CCTK_WARN(0, "Error registering boundaries");
+}