Regenerate code

This was regenerated with commit 66668c409754acc1f76ed50118d8831fc3992c59 of Kranc and passes all existing test suites.

1 files changed, 8 insertions, 1 deletions

diff --git a/ML_BSSN/src/ML_BSSN_enforce.cc b/ML_BSSN/src/ML_BSSN_enforce.cc
index 54d04bc..485b834 100644
--- a/ML_BSSN/src/ML_BSSN_enforce.cc
+++ b/ML_BSSN/src/ML_BSSN_enforce.cc
@@ -6,6 +6,7 @@
 #include <math.h>
 #include <stdio.h>
 #include <stdlib.h>
+#include <string.h>
 #include "cctk.h"
 #include "cctk_Arguments.h"
 #include "cctk_Parameters.h"
@@ -41,15 +42,20 @@ static void ML_BSSN_enforce_Body(cGH const * restrict const cctkGH, int const di
   const char *groups[] = {"ML_BSSN::ML_curv","ML_BSSN::ML_lapse","ML_BSSN::ML_metric"};
   GenericFD_AssertGroupStorage(cctkGH, "ML_BSSN_enforce", 3, groups);
   
+  
   /* Include user-supplied include files */
   
   /* Initialise finite differencing variables */
   ptrdiff_t const di = 1;
   ptrdiff_t const dj = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0);
   ptrdiff_t const dk = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0);
+  ptrdiff_t const cdi = sizeof(CCTK_REAL) * di;
+  ptrdiff_t const cdj = sizeof(CCTK_REAL) * dj;
+  ptrdiff_t const cdk = sizeof(CCTK_REAL) * dk;
   CCTK_REAL const dx = ToReal(CCTK_DELTA_SPACE(0));
   CCTK_REAL const dy = ToReal(CCTK_DELTA_SPACE(1));
   CCTK_REAL const dz = ToReal(CCTK_DELTA_SPACE(2));
+  CCTK_REAL const dt = ToReal(CCTK_DELTA_TIME);
   CCTK_REAL const dxi = INV(dx);
   CCTK_REAL const dyi = INV(dy);
   CCTK_REAL const dzi = INV(dz);
@@ -91,6 +97,7 @@ static void ML_BSSN_enforce_Body(cGH const * restrict const cctkGH, int const di
     ptrdiff_t const index = di*i + dj*j + dk*k;
     
     /* Assign local copies of grid functions */
+    
     CCTK_REAL alphaL = alpha[index];
     CCTK_REAL At11L = At11[index];
     CCTK_REAL At12L = At12[index];
@@ -105,6 +112,7 @@ static void ML_BSSN_enforce_Body(cGH const * restrict const cctkGH, int const di
     CCTK_REAL gt23L = gt23[index];
     CCTK_REAL gt33L = gt33[index];
     
+    
     /* Include user supplied include files */
     
     /* Precompute derivatives */
@@ -141,7 +149,6 @@ static void ML_BSSN_enforce_Body(cGH const * restrict const cctkGH, int const di
     
     alphaL = fmax(alphaL,ToReal(MinimumLapse));
     
-    
     /* Copy local copies back to grid functions */
     alpha[index] = alphaL;
     At11[index] = At11L;