The carpet patch timed Sat Mar 4 18:38:10 CET 2006, not allowing any longer the use of CCTK_DELTA_SPACE in global mode, had broken our way of computing volume integrals. Here is the upgrade.

Further changes are the transformation of divisions into moltiplications and the boundary corrections in the computation of the integral when PUGH is used. This last part is not tested (but will be), anyway who uses PUGH nowadays?


git-svn-id: http://svn.aei.mpg.de/numrel/AEIThorns/ADMMass/trunk@9 54511f98-0e4f-0410-826e-eb8b393f5a1e

1 files changed, 9 insertions, 1 deletions

diff --git a/schedule.ccl b/schedule.ccl
index 42ebdd0..fc54cb1 100644
--- a/schedule.ccl
+++ b/schedule.ccl
@@ -22,11 +22,19 @@ schedule ADMMass_Loop IN ADMMass BEFORE ADMMass_Local
   OPTIONS: global
 } "Decrement loop counter"
 
-
+#########################################################################################
+# We must schedule the local routines to compute the integrals in global-loop-local 
+# (as opposed to local) mode, in order to make sure that the scheduling condition
+# "AFTER", which describes the ADMMass group, is respected. This group may depend 
+# on local routines (excision, emask) and must be run after all local routines. 
+# If we had no such a dependence, we could have scheduled the integral computations 
+# simply in local mode.
+##########################################################################################
 schedule GROUP ADMMass AT POSTSTEP AFTER ADMMass_InitLoopCounter WHILE ADMMass::ADMMass_LoopCounter
 {
   STORAGE:ADMMass_GFs[3]
   STORAGE:ADMMass_box
+  STORAGE:grid_spacing_product
 } "ADMMass loop"
 
 schedule ADMMass_Surface IN ADMMass