Re-generate code

4 files changed, 7 insertions, 24 deletions

diff --git a/ML_ADMConstraints/param.ccl b/ML_ADMConstraints/param.ccl
index d35f62d..460ea74 100644
--- a/ML_ADMConstraints/param.ccl
+++ b/ML_ADMConstraints/param.ccl
@@ -3,11 +3,6 @@
 
 shares: GenericFD
 
-USES CCTK_INT stencil_width
-USES CCTK_INT stencil_width_x
-USES CCTK_INT stencil_width_y
-USES CCTK_INT stencil_width_z
-USES CCTK_INT boundary_width
 
 
 shares: MethodOfLines
@@ -15,37 +10,37 @@ shares: MethodOfLines
 USES CCTK_INT MoL_Num_Evolved_Vars
 
 restricted:
-CCTK_INT verbose "verbose"
+CCTK_INT verbose "verbose" STEERABLE=ALWAYS
 {
   *:* :: ""
 } 0
 
 restricted:
-CCTK_INT ML_ADMConstraints_MaxNumEvolvedVars "Number of evolved variables used by this thorn" ACCUMULATOR-BASE=MethodofLines::MoL_Num_Evolved_Vars
+CCTK_INT ML_ADMConstraints_MaxNumEvolvedVars "Number of evolved variables used by this thorn" ACCUMULATOR-BASE=MethodofLines::MoL_Num_Evolved_Vars STEERABLE=RECOVER
 {
   0:0 :: "Number of evolved variables used by this thorn"
 } 0
 
 restricted:
-CCTK_INT timelevels "Number of active timelevels"
+CCTK_INT timelevels "Number of active timelevels" STEERABLE=RECOVER
 {
   0:3 :: ""
 } 3
 
 restricted:
-CCTK_INT rhs_timelevels "Number of active RHS timelevels"
+CCTK_INT rhs_timelevels "Number of active RHS timelevels" STEERABLE=RECOVER
 {
   0:3 :: ""
 } 1
 
 restricted:
-CCTK_INT ML_ADMConstraints_calc_every "ML_ADMConstraints_calc_every"
+CCTK_INT ML_ADMConstraints_calc_every "ML_ADMConstraints_calc_every" STEERABLE=ALWAYS
 {
   *:* :: ""
 } 1
 
 restricted:
-CCTK_INT ML_ADMConstraints_calc_offset "ML_ADMConstraints_calc_offset"
+CCTK_INT ML_ADMConstraints_calc_offset "ML_ADMConstraints_calc_offset" STEERABLE=ALWAYS
 {
   *:* :: ""
 } 0
diff --git a/ML_ADMConstraints/schedule.ccl b/ML_ADMConstraints/schedule.ccl
index 74b3d66..fa471af 100644
--- a/ML_ADMConstraints/schedule.ccl
+++ b/ML_ADMConstraints/schedule.ccl
@@ -73,7 +73,7 @@ schedule ML_ADMConstraints_CheckBoundaries at BASEGRID
   OPTIONS: meta
 } "check boundaries treatment"
 
-schedule group ApplyBCs as ML_ADMConstraints_ApplyBCs in MoL_PseudoEvolution after ML_ADMConstraints_SelectBoundConds
+schedule group ApplyBCs as ML_ADMConstraints_ApplyBCs in MoL_PseudoEvolution  after ML_ADMConstraints_SelectBoundConds
 {
   # no language specified
 } "Apply boundary conditions controlled by thorn Boundary"
diff --git a/ML_ADMConstraints/src/Differencing.h b/ML_ADMConstraints/src/Differencing.h
index 1cdce99..aa21060 100644
--- a/ML_ADMConstraints/src/Differencing.h
+++ b/ML_ADMConstraints/src/Differencing.h
@@ -10,9 +10,3 @@
 #define PDstandardNth23(u,i,j,k) (p1o144dydz*(-64*((u)[index+di*(0)+dj*(-1)+dk*(1)] + (u)[index+di*(0)+dj*(1)+dk*(-1)]) + 64*((u)[index+di*(0)+dj*(-1)+dk*(-1)] + (u)[index+di*(0)+dj*(1)+dk*(1)]) + 8*((u)[index+di*(0)+dj*(-1)+dk*(2)] + (u)[index+di*(0)+dj*(1)+dk*(-2)] + (u)[index+di*(0)+dj*(-2)+dk*(1)] + (u)[index+di*(0)+dj*(2)+dk*(-1)]) - 8*((u)[index+di*(0)+dj*(-1)+dk*(-2)] + (u)[index+di*(0)+dj*(1)+dk*(2)] + (u)[index+di*(0)+dj*(-2)+dk*(-1)] + (u)[index+di*(0)+dj*(2)+dk*(1)]) + (u)[index+di*(0)+dj*(-2)+dk*(-2)] - (u)[index+di*(0)+dj*(-2)+dk*(2)] - (u)[index+di*(0)+dj*(2)+dk*(-2)] + (u)[index+di*(0)+dj*(2)+dk*(2)]))
 #define PDstandardNth31(u,i,j,k) (p1o144dxdz*(-64*((u)[index+di*(-1)+dj*(0)+dk*(1)] + (u)[index+di*(1)+dj*(0)+dk*(-1)]) + 64*((u)[index+di*(-1)+dj*(0)+dk*(-1)] + (u)[index+di*(1)+dj*(0)+dk*(1)]) + 8*((u)[index+di*(-1)+dj*(0)+dk*(2)] + (u)[index+di*(1)+dj*(0)+dk*(-2)] + (u)[index+di*(-2)+dj*(0)+dk*(1)] + (u)[index+di*(2)+dj*(0)+dk*(-1)]) - 8*((u)[index+di*(-1)+dj*(0)+dk*(-2)] + (u)[index+di*(1)+dj*(0)+dk*(2)] + (u)[index+di*(-2)+dj*(0)+dk*(-1)] + (u)[index+di*(2)+dj*(0)+dk*(1)]) + (u)[index+di*(-2)+dj*(0)+dk*(-2)] - (u)[index+di*(-2)+dj*(0)+dk*(2)] - (u)[index+di*(2)+dj*(0)+dk*(-2)] + (u)[index+di*(2)+dj*(0)+dk*(2)]))
 #define PDstandardNth32(u,i,j,k) (p1o144dydz*(-64*((u)[index+di*(0)+dj*(-1)+dk*(1)] + (u)[index+di*(0)+dj*(1)+dk*(-1)]) + 64*((u)[index+di*(0)+dj*(-1)+dk*(-1)] + (u)[index+di*(0)+dj*(1)+dk*(1)]) + 8*((u)[index+di*(0)+dj*(-1)+dk*(2)] + (u)[index+di*(0)+dj*(1)+dk*(-2)] + (u)[index+di*(0)+dj*(-2)+dk*(1)] + (u)[index+di*(0)+dj*(2)+dk*(-1)]) - 8*((u)[index+di*(0)+dj*(-1)+dk*(-2)] + (u)[index+di*(0)+dj*(1)+dk*(2)] + (u)[index+di*(0)+dj*(-2)+dk*(-1)] + (u)[index+di*(0)+dj*(2)+dk*(1)]) + (u)[index+di*(0)+dj*(-2)+dk*(-2)] - (u)[index+di*(0)+dj*(-2)+dk*(2)] - (u)[index+di*(0)+dj*(2)+dk*(-2)] + (u)[index+di*(0)+dj*(2)+dk*(2)]))
-#define PDupwindNth1(u,i,j,k) (p1o12dx*(-10*(u)[index+di*(0)+dj*(0)+dk*(0)] - 6*(u)[index+di*(2*dir1)+dj*(0)+dk*(0)] + (u)[index+di*(3*dir1)+dj*(0)+dk*(0)] - 3*(u)[index+di*(-dir1)+dj*(0)+dk*(0)] + 18*(u)[index+di*(dir1)+dj*(0)+dk*(0)])*dir1)
-#define PDupwindNth2(u,i,j,k) (p1o12dy*(-10*(u)[index+di*(0)+dj*(0)+dk*(0)] - 6*(u)[index+di*(0)+dj*(2*dir2)+dk*(0)] + (u)[index+di*(0)+dj*(3*dir2)+dk*(0)] - 3*(u)[index+di*(0)+dj*(-dir2)+dk*(0)] + 18*(u)[index+di*(0)+dj*(dir2)+dk*(0)])*dir2)
-#define PDupwindNth3(u,i,j,k) (p1o12dz*(-10*(u)[index+di*(0)+dj*(0)+dk*(0)] - 6*(u)[index+di*(0)+dj*(0)+dk*(2*dir3)] + (u)[index+di*(0)+dj*(0)+dk*(3*dir3)] - 3*(u)[index+di*(0)+dj*(0)+dk*(-dir3)] + 18*(u)[index+di*(0)+dj*(0)+dk*(dir3)])*dir3)
-#define PDonesided1(u,i,j,k) (p1odx*(-(u)[index+di*(0)+dj*(0)+dk*(0)] + (u)[index+di*(dir1)+dj*(0)+dk*(0)])*dir1)
-#define PDonesided2(u,i,j,k) (p1ody*(-(u)[index+di*(0)+dj*(0)+dk*(0)] + (u)[index+di*(0)+dj*(dir2)+dk*(0)])*dir2)
-#define PDonesided3(u,i,j,k) (p1odz*(-(u)[index+di*(0)+dj*(0)+dk*(0)] + (u)[index+di*(0)+dj*(0)+dk*(dir3)])*dir3)
diff --git a/ML_ADMConstraints/src/ML_ADMConstraints.c b/ML_ADMConstraints/src/ML_ADMConstraints.c
index 54035b8..44990c1 100644
--- a/ML_ADMConstraints/src/ML_ADMConstraints.c
+++ b/ML_ADMConstraints/src/ML_ADMConstraints.c
@@ -53,9 +53,6 @@ void ML_ADMConstraints_Body(cGH const * restrict const cctkGH, int const dir, in
     return;
   }
   
-  const char *groups[] = {"ADMBase::curv","ADMBase::lapse","ADMBase::metric","ADMBase::shift","ML_ADMConstraints::ML_Ham","ML_ADMConstraints::ML_mom"};
-  GenericFD_AssertGroupStorage(cctkGH, "ML_ADMConstraints", 6, groups);
-  
   /* Include user-supplied include files */
   
   /* Initialise finite differencing variables */
@@ -84,9 +81,6 @@ void ML_ADMConstraints_Body(cGH const * restrict const cctkGH, int const dir, in
   CCTK_REAL const p1o144dxdy = (INV(dx)*INV(dy))/144.;
   CCTK_REAL const p1o144dxdz = (INV(dx)*INV(dz))/144.;
   CCTK_REAL const p1o144dydz = (INV(dy)*INV(dz))/144.;
-  CCTK_REAL const p1odx = INV(dx);
-  CCTK_REAL const p1ody = INV(dy);
-  CCTK_REAL const p1odz = INV(dz);
   CCTK_REAL const pm1o12dx2 = -pow(dx,-2)/12.;
   CCTK_REAL const pm1o12dy2 = -pow(dy,-2)/12.;
   CCTK_REAL const pm1o12dz2 = -pow(dz,-2)/12.;