Schedule constraint calculation slighly differently when

constraints_persist=yes:

Calculate constraints at POSTINITIAL (after MoL_PostInitial and
MoL_PostStep) and at EVOL (after MoL_Evolution).

Re-apply the constraint boundary conditions at POSTREGRID and
POSTRESTRICT.

(This method is closer to the way in which time-evolved quantities are
handled.)

Introduce a new schedule group ADMConstraintsBoundariesGroup which
applies only the constraint boundaries.


git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/ADMConstraints/trunk@129 b7a48df3-cbbf-4440-997f-b4b717c9f7fc

1 files changed, 16 insertions, 4 deletions

diff --git a/schedule.ccl b/schedule.ccl
index 2004588..8ae3cee 100644
--- a/schedule.ccl
+++ b/schedule.ccl
@@ -34,13 +34,21 @@ if (constraints_persist)
     STORAGE: hamiltonian[3], normalized_hamiltonian[3], momentum[3]
   }
   
-  schedule GROUP ADMConstraintsGroup at CCTK_EVOL after MoL_Evolution
+  schedule GROUP ADMConstraintsGroup at CCTK_POSTINITIAL after (MoL_PostStep MoL_PostInitial)
   {
   } "Evaluate ADM constraints, and perform symmetry boundary conditions"
   
-  schedule GROUP ADMConstraintsGroup at CCTK_POSTREGRID
+  schedule GROUP ADMConstraintsGroup at CCTK_EVOL after MoL_Evolution
   {
   } "Evaluate ADM constraints, and perform symmetry boundary conditions"
+  
+  schedule GROUP ADMConstraintsBoundaryGroup at CCTK_POSTREGRID
+  {
+  } "Set ADM constraints on the boundary"
+  
+  schedule GROUP ADMConstraintsBoundaryGroup at CCTK_POSTRESTRICT
+  {
+  } "Set ADM constraints on the boundary"
 }
 else
 {
@@ -56,13 +64,17 @@ schedule ADMConstraints in ADMConstraintsGroup
   LANG: Fortran
 } "Evaluate ADM constraints"
 
-schedule ADMConstraints_Boundaries in ADMConstraintsGroup after ADMConstraints
+schedule GROUP ADMConstraintsBoundariesGroup in ADMConstraintsGroup after ADMConstraints
+{
+} "Set ADM constraints on the boundary"
+
+schedule ADMConstraints_Boundaries in ADMConstraintsBoundariesGroup
 {
   LANG: Fortran
   OPTIONS: level
   SYNC: hamiltonian, normalized_hamiltonian, momentum
 } "Select boundary conditions for the ADM constraints"
 
-schedule GROUP ApplyBCs as ADMConstraints_ApplyBCs in ADMConstraintsGroup after ADMConstraints_Boundaries
+schedule GROUP ApplyBCs as ADMConstraints_ApplyBCs in ADMConstraintsBoundariesGroup after ADMConstraints_Boundaries
 {
 } "Apply boundary conditions to the ADM constraints"