Merge branch 'master' of ssh://relativity.phys.lsu.edu/home/perturbed/gitroot/McLachlan

Conflicts:

	m/Makefile

14 files changed, 113 insertions, 54 deletions

diff --git a/ML_BSSN/src/Boundaries.c b/ML_BSSN/src/Boundaries.c
index f073be3..4220352 100644
--- a/ML_BSSN/src/Boundaries.c
+++ b/ML_BSSN/src/Boundaries.c
@@ -1,5 +1,5 @@
-/*  File produced by user eschnett */
-/*  Produced with Mathematica Version 6.0 for Mac OS X x86 (32-bit) (April 20, 2007) */
+/*  File produced by user diener */
+/*  Produced with Mathematica Version 6.0 for Linux x86 (32-bit) (April 20, 2007) */
 
 /*  Mathematica script written by Ian Hinder and Sascha Husa */
 
diff --git a/ML_BSSN/src/ML_BSSN_Minkowski.c b/ML_BSSN/src/ML_BSSN_Minkowski.c
index aaf5dd8..0c58c32 100644
--- a/ML_BSSN/src/ML_BSSN_Minkowski.c
+++ b/ML_BSSN/src/ML_BSSN_Minkowski.c
@@ -1,5 +1,5 @@
-/*  File produced by user eschnett */
-/*  Produced with Mathematica Version 6.0 for Mac OS X x86 (32-bit) (April 20, 2007) */
+/*  File produced by user diener */
+/*  Produced with Mathematica Version 6.0 for Linux x86 (32-bit) (April 20, 2007) */
 
 /*  Mathematica script written by Ian Hinder and Sascha Husa */
 
diff --git a/ML_BSSN/src/ML_BSSN_RHS.c b/ML_BSSN/src/ML_BSSN_RHS.c
index 5890a66..ca53f5e 100644
--- a/ML_BSSN/src/ML_BSSN_RHS.c
+++ b/ML_BSSN/src/ML_BSSN_RHS.c
@@ -1,5 +1,5 @@
-/*  File produced by user eschnett */
-/*  Produced with Mathematica Version 6.0 for Mac OS X x86 (32-bit) (April 20, 2007) */
+/*  File produced by user diener */
+/*  Produced with Mathematica Version 6.0 for Linux x86 (32-bit) (April 20, 2007) */
 
 /*  Mathematica script written by Ian Hinder and Sascha Husa */
 
@@ -170,6 +170,15 @@ void ML_BSSN_RHS_Body(cGH *cctkGH, CCTK_INT dir, CCTK_INT face, CCTK_REAL normal
     CCTK_REAL PDstandardNth1At33 = INITVALUE;
     CCTK_REAL PDstandardNth2At33 = INITVALUE;
     CCTK_REAL PDstandardNth3At33 = INITVALUE;
+    CCTK_REAL PDstandardNth1B1 = INITVALUE;
+    CCTK_REAL PDstandardNth2B1 = INITVALUE;
+    CCTK_REAL PDstandardNth3B1 = INITVALUE;
+    CCTK_REAL PDstandardNth1B2 = INITVALUE;
+    CCTK_REAL PDstandardNth2B2 = INITVALUE;
+    CCTK_REAL PDstandardNth3B2 = INITVALUE;
+    CCTK_REAL PDstandardNth1B3 = INITVALUE;
+    CCTK_REAL PDstandardNth2B3 = INITVALUE;
+    CCTK_REAL PDstandardNth3B3 = INITVALUE;
     CCTK_REAL PDstandardNth1beta1 = INITVALUE;
     CCTK_REAL PDstandardNth2beta1 = INITVALUE;
     CCTK_REAL PDstandardNth3beta1 = INITVALUE;
@@ -366,6 +375,15 @@ void ML_BSSN_RHS_Body(cGH *cctkGH, CCTK_INT dir, CCTK_INT face, CCTK_REAL normal
     PDstandardNth1At33 = PDstandardNth1(At33, i, j, k);
     PDstandardNth2At33 = PDstandardNth2(At33, i, j, k);
     PDstandardNth3At33 = PDstandardNth3(At33, i, j, k);
+    PDstandardNth1B1 = PDstandardNth1(B1, i, j, k);
+    PDstandardNth2B1 = PDstandardNth2(B1, i, j, k);
+    PDstandardNth3B1 = PDstandardNth3(B1, i, j, k);
+    PDstandardNth1B2 = PDstandardNth1(B2, i, j, k);
+    PDstandardNth2B2 = PDstandardNth2(B2, i, j, k);
+    PDstandardNth3B2 = PDstandardNth3(B2, i, j, k);
+    PDstandardNth1B3 = PDstandardNth1(B3, i, j, k);
+    PDstandardNth2B3 = PDstandardNth2(B3, i, j, k);
+    PDstandardNth3B3 = PDstandardNth3(B3, i, j, k);
     PDstandardNth1beta1 = PDstandardNth1(beta1, i, j, k);
     PDstandardNth2beta1 = PDstandardNth2(beta1, i, j, k);
     PDstandardNth3beta1 = PDstandardNth3(beta1, i, j, k);
@@ -1112,22 +1130,32 @@ void ML_BSSN_RHS_Body(cGH *cctkGH, CCTK_INT dir, CCTK_INT face, CCTK_REAL normal
         2*(At13L*PDstandardNth3beta1 + At23L*PDstandardNth3beta2 + At33L*PDstandardNth3beta3) + 
         em4phi*(Ats33 - g33*kthird*trAts) + alphaL*(-2*(At13L*Atm13 + At23L*Atm23 + At33L*Atm33) + At33L*trKL);
     
-    alpharhsL  =  beta1L*PDstandardNth1alpha + beta2L*PDstandardNth2alpha + beta3L*PDstandardNth3alpha - 
-        AL*harmonicF*pow(alphaL,harmonicN);
+    alpharhsL  =  LapseAdvectionCoeff*(beta1L*PDstandardNth1alpha + beta2L*PDstandardNth2alpha + 
+           beta3L*PDstandardNth3alpha) + harmonicF*(AL*(-1 + LapseAdvectionCoeff) - LapseAdvectionCoeff*trKL)*
+         pow(alphaL,harmonicN);
     
-    ArhsL  =  -(AL*AlphaDriver) + trKrhsL;
+    ArhsL  =  (-1 + LapseAdvectionCoeff)*(AL*AlphaDriver - trKrhsL);
     
-    beta1rhsL  =  B1L*ShiftGammaCoeff*pow(alphaL,ShiftAlphaPower);
+    beta1rhsL  =  (beta1L*PDstandardNth1beta1 + beta2L*PDstandardNth2beta1 + beta3L*PDstandardNth3beta1)*
+         ShiftAdvectionCoeff + B1L*ShiftGammaCoeff;
     
-    beta2rhsL  =  B2L*ShiftGammaCoeff*pow(alphaL,ShiftAlphaPower);
+    beta2rhsL  =  (beta1L*PDstandardNth1beta2 + beta2L*PDstandardNth2beta2 + beta3L*PDstandardNth3beta2)*
+         ShiftAdvectionCoeff + B2L*ShiftGammaCoeff;
     
-    beta3rhsL  =  B3L*ShiftGammaCoeff*pow(alphaL,ShiftAlphaPower);
+    beta3rhsL  =  (beta1L*PDstandardNth1beta3 + beta2L*PDstandardNth2beta3 + beta3L*PDstandardNth3beta3)*
+         ShiftAdvectionCoeff + B3L*ShiftGammaCoeff;
     
-    B1rhsL  =  -(B1L*BetaDriver) + Xt1rhsL;
+    B1rhsL  =  -(B1L*BetaDriver) + (beta1L*(PDstandardNth1B1 - PDstandardNth1Xt1) + 
+           beta2L*(PDstandardNth2B1 - PDstandardNth2Xt1) + beta3L*(PDstandardNth3B1 - PDstandardNth3Xt1))*
+         ShiftAdvectionCoeff + Xt1rhsL;
     
-    B2rhsL  =  -(B2L*BetaDriver) + Xt2rhsL;
+    B2rhsL  =  -(B2L*BetaDriver) + (beta1L*(PDstandardNth1B2 - PDstandardNth1Xt2) + 
+           beta2L*(PDstandardNth2B2 - PDstandardNth2Xt2) + beta3L*(PDstandardNth3B2 - PDstandardNth3Xt2))*
+         ShiftAdvectionCoeff + Xt2rhsL;
     
-    B3rhsL  =  -(B3L*BetaDriver) + Xt3rhsL;
+    B3rhsL  =  -(B3L*BetaDriver) + (beta1L*(PDstandardNth1B3 - PDstandardNth1Xt3) + 
+           beta2L*(PDstandardNth2B3 - PDstandardNth2Xt3) + beta3L*(PDstandardNth3B3 - PDstandardNth3Xt3))*
+         ShiftAdvectionCoeff + Xt3rhsL;
     
     
     /* Copy local copies back to grid functions */
diff --git a/ML_BSSN/src/ML_BSSN_constraints.c b/ML_BSSN/src/ML_BSSN_constraints.c
index cf6b411..c3883d6 100644
--- a/ML_BSSN/src/ML_BSSN_constraints.c
+++ b/ML_BSSN/src/ML_BSSN_constraints.c
@@ -1,5 +1,5 @@
-/*  File produced by user eschnett */
-/*  Produced with Mathematica Version 6.0 for Mac OS X x86 (32-bit) (April 20, 2007) */
+/*  File produced by user diener */
+/*  Produced with Mathematica Version 6.0 for Linux x86 (32-bit) (April 20, 2007) */
 
 /*  Mathematica script written by Ian Hinder and Sascha Husa */
 
diff --git a/ML_BSSN/src/ML_BSSN_convertFromADMBase.c b/ML_BSSN/src/ML_BSSN_convertFromADMBase.c
index 36688ce..9cefffc 100644
--- a/ML_BSSN/src/ML_BSSN_convertFromADMBase.c
+++ b/ML_BSSN/src/ML_BSSN_convertFromADMBase.c
@@ -1,5 +1,5 @@
-/*  File produced by user eschnett */
-/*  Produced with Mathematica Version 6.0 for Mac OS X x86 (32-bit) (April 20, 2007) */
+/*  File produced by user diener */
+/*  Produced with Mathematica Version 6.0 for Linux x86 (32-bit) (April 20, 2007) */
 
 /*  Mathematica script written by Ian Hinder and Sascha Husa */
 
diff --git a/ML_BSSN/src/ML_BSSN_convertFromADMBaseGamma.c b/ML_BSSN/src/ML_BSSN_convertFromADMBaseGamma.c
index af99ac2..f7fc8ad 100644
--- a/ML_BSSN/src/ML_BSSN_convertFromADMBaseGamma.c
+++ b/ML_BSSN/src/ML_BSSN_convertFromADMBaseGamma.c
@@ -1,5 +1,5 @@
-/*  File produced by user eschnett */
-/*  Produced with Mathematica Version 6.0 for Mac OS X x86 (32-bit) (April 20, 2007) */
+/*  File produced by user diener */
+/*  Produced with Mathematica Version 6.0 for Linux x86 (32-bit) (April 20, 2007) */
 
 /*  Mathematica script written by Ian Hinder and Sascha Husa */
 
@@ -114,9 +114,15 @@ void ML_BSSN_convertFromADMBaseGamma_Body(cGH *cctkGH, CCTK_INT dir, CCTK_INT fa
     /* Declare precomputed derivatives*/
     
     /* Declare derivatives */
-    CCTK_REAL PDstandardNth1alpha = INITVALUE;
-    CCTK_REAL PDstandardNth2alpha = INITVALUE;
-    CCTK_REAL PDstandardNth3alpha = INITVALUE;
+    CCTK_REAL PDstandardNth1B1 = INITVALUE;
+    CCTK_REAL PDstandardNth2B1 = INITVALUE;
+    CCTK_REAL PDstandardNth3B1 = INITVALUE;
+    CCTK_REAL PDstandardNth1B2 = INITVALUE;
+    CCTK_REAL PDstandardNth2B2 = INITVALUE;
+    CCTK_REAL PDstandardNth3B2 = INITVALUE;
+    CCTK_REAL PDstandardNth1B3 = INITVALUE;
+    CCTK_REAL PDstandardNth2B3 = INITVALUE;
+    CCTK_REAL PDstandardNth3B3 = INITVALUE;
     CCTK_REAL PDstandardNth1gt11 = INITVALUE;
     CCTK_REAL PDstandardNth2gt11 = INITVALUE;
     CCTK_REAL PDstandardNth3gt11 = INITVALUE;
@@ -138,6 +144,9 @@ void ML_BSSN_convertFromADMBaseGamma_Body(cGH *cctkGH, CCTK_INT dir, CCTK_INT fa
     
     /* Assign local copies of grid functions */
     alphaL = alpha[index];
+    B1L = B1[index];
+    B2L = B2[index];
+    B3L = B3[index];
     beta1L = beta1[index];
     beta2L = beta2[index];
     beta3L = beta3[index];
@@ -157,9 +166,15 @@ void ML_BSSN_convertFromADMBaseGamma_Body(cGH *cctkGH, CCTK_INT dir, CCTK_INT fa
     /* Include user supplied include files */
     
     /* Precompute derivatives (new style) */
-    PDstandardNth1alpha = PDstandardNth1(alpha, i, j, k);
-    PDstandardNth2alpha = PDstandardNth2(alpha, i, j, k);
-    PDstandardNth3alpha = PDstandardNth3(alpha, i, j, k);
+    PDstandardNth1B1 = PDstandardNth1(B1, i, j, k);
+    PDstandardNth2B1 = PDstandardNth2(B1, i, j, k);
+    PDstandardNth3B1 = PDstandardNth3(B1, i, j, k);
+    PDstandardNth1B2 = PDstandardNth1(B2, i, j, k);
+    PDstandardNth2B2 = PDstandardNth2(B2, i, j, k);
+    PDstandardNth3B2 = PDstandardNth3(B2, i, j, k);
+    PDstandardNth1B3 = PDstandardNth1(B3, i, j, k);
+    PDstandardNth2B3 = PDstandardNth2(B3, i, j, k);
+    PDstandardNth3B3 = PDstandardNth3(B3, i, j, k);
     PDstandardNth1gt11 = PDstandardNth1(gt11, i, j, k);
     PDstandardNth2gt11 = PDstandardNth2(gt11, i, j, k);
     PDstandardNth3gt11 = PDstandardNth3(gt11, i, j, k);
@@ -256,14 +271,16 @@ void ML_BSSN_convertFromADMBaseGamma_Body(cGH *cctkGH, CCTK_INT dir, CCTK_INT fa
     
     Xt3L  =  Gt311*gtu11 + Gt322*gtu22 + 2*(Gt312*gtu21 + Gt313*gtu31 + Gt323*gtu32) + Gt333*gtu33;
     
-    AL  =  (-dtalpL + beta1L*PDstandardNth1alpha + beta2L*PDstandardNth2alpha + beta3L*PDstandardNth3alpha)*INV(harmonicF)*
-        pow(alphaL,-harmonicN);
+    AL  =  -1.*dtalpL*(-1. + LapseAdvectionCoeff)*INV(harmonicF)*pow(alphaL,-harmonicN);
     
-    B1L  =  dtbetaxL*ShiftGammaCoeff*INV(1.e-100 + SQR(ShiftGammaCoeff))*pow(alphaL,-ShiftAlphaPower);
+    B1L  =  (dtbetaxL - (beta1L*PDstandardNth1B1 + beta2L*PDstandardNth2B1 + beta3L*PDstandardNth3B1)*ShiftAdvectionCoeff)*
+        INV(ShiftGammaCoeff);
     
-    B2L  =  dtbetayL*ShiftGammaCoeff*INV(1.e-100 + SQR(ShiftGammaCoeff))*pow(alphaL,-ShiftAlphaPower);
+    B2L  =  (dtbetayL - (beta1L*PDstandardNth1B2 + beta2L*PDstandardNth2B2 + beta3L*PDstandardNth3B2)*ShiftAdvectionCoeff)*
+        INV(ShiftGammaCoeff);
     
-    B3L  =  dtbetazL*ShiftGammaCoeff*INV(1.e-100 + SQR(ShiftGammaCoeff))*pow(alphaL,-ShiftAlphaPower);
+    B3L  =  (dtbetazL - (beta1L*PDstandardNth1B3 + beta2L*PDstandardNth2B3 + beta3L*PDstandardNth3B3)*ShiftAdvectionCoeff)*
+        INV(ShiftGammaCoeff);
     
     
     /* Copy local copies back to grid functions */
diff --git a/ML_BSSN/src/ML_BSSN_convertToADMBase.c b/ML_BSSN/src/ML_BSSN_convertToADMBase.c
index c6071ee..14359b2 100644
--- a/ML_BSSN/src/ML_BSSN_convertToADMBase.c
+++ b/ML_BSSN/src/ML_BSSN_convertToADMBase.c
@@ -1,5 +1,5 @@
-/*  File produced by user eschnett */
-/*  Produced with Mathematica Version 6.0 for Mac OS X x86 (32-bit) (April 20, 2007) */
+/*  File produced by user diener */
+/*  Produced with Mathematica Version 6.0 for Linux x86 (32-bit) (April 20, 2007) */
 
 /*  Mathematica script written by Ian Hinder and Sascha Husa */
 
@@ -133,6 +133,9 @@ void ML_BSSN_convertToADMBase_Body(cGH *cctkGH, CCTK_INT dir, CCTK_INT face, CCT
     CCTK_REAL PDstandardNth1alpha = INITVALUE;
     CCTK_REAL PDstandardNth2alpha = INITVALUE;
     CCTK_REAL PDstandardNth3alpha = INITVALUE;
+    CCTK_REAL PDstandardNth1beta1 = INITVALUE;
+    CCTK_REAL PDstandardNth2beta1 = INITVALUE;
+    CCTK_REAL PDstandardNth3beta1 = INITVALUE;
     
     /* Assign local copies of grid functions */
     AL = A[index];
@@ -166,6 +169,9 @@ void ML_BSSN_convertToADMBase_Body(cGH *cctkGH, CCTK_INT dir, CCTK_INT face, CCT
     PDstandardNth1alpha = PDstandardNth1(alpha, i, j, k);
     PDstandardNth2alpha = PDstandardNth2(alpha, i, j, k);
     PDstandardNth3alpha = PDstandardNth3(alpha, i, j, k);
+    PDstandardNth1beta1 = PDstandardNth1(beta1, i, j, k);
+    PDstandardNth2beta1 = PDstandardNth2(beta1, i, j, k);
+    PDstandardNth3beta1 = PDstandardNth3(beta1, i, j, k);
     
     /* Precompute derivatives (old style) */
     
@@ -228,14 +234,17 @@ void ML_BSSN_convertToADMBase_Body(cGH *cctkGH, CCTK_INT dir, CCTK_INT face, CCT
     
     betazL  =  beta3L;
     
-    dtalpL  =  beta1L*PDstandardNth1alpha + beta2L*PDstandardNth2alpha + beta3L*PDstandardNth3alpha - 
-        AL*harmonicF*pow(alphaL,harmonicN);
+    dtalpL  =  LapseAdvectionCoeff*(beta1L*PDstandardNth1alpha + beta2L*PDstandardNth2alpha + beta3L*PDstandardNth3alpha) + 
+        harmonicF*(AL*(-1 + LapseAdvectionCoeff) - LapseAdvectionCoeff*trKL)*pow(alphaL,harmonicN);
     
-    dtbetaxL  =  B1L*ShiftGammaCoeff*pow(alphaL,ShiftAlphaPower);
+    dtbetaxL  =  (beta1L*PDstandardNth1beta1 + beta2L*PDstandardNth2beta1 + beta3L*PDstandardNth3beta1)*
+         ShiftAdvectionCoeff + B1L*ShiftGammaCoeff;
     
-    dtbetayL  =  B2L*ShiftGammaCoeff*pow(alphaL,ShiftAlphaPower);
+    dtbetayL  =  (beta1L*PDstandardNth1beta1 + beta2L*PDstandardNth2beta1 + beta3L*PDstandardNth3beta1)*
+         ShiftAdvectionCoeff + B2L*ShiftGammaCoeff;
     
-    dtbetazL  =  B3L*ShiftGammaCoeff*pow(alphaL,ShiftAlphaPower);
+    dtbetazL  =  (beta1L*PDstandardNth1beta1 + beta2L*PDstandardNth2beta1 + beta3L*PDstandardNth3beta1)*
+         ShiftAdvectionCoeff + B3L*ShiftGammaCoeff;
     
     
     /* Copy local copies back to grid functions */
@@ -270,5 +279,5 @@ void ML_BSSN_convertToADMBase(CCTK_ARGUMENTS)
   DECLARE_CCTK_ARGUMENTS
   DECLARE_CCTK_PARAMETERS
   
-  GenericFD_LoopOverInterior(cctkGH, &ML_BSSN_convertToADMBase_Body);
+  GenericFD_LoopOverEverything(cctkGH, &ML_BSSN_convertToADMBase_Body);
 }
diff --git a/ML_BSSN/src/ML_BSSN_enforce.c b/ML_BSSN/src/ML_BSSN_enforce.c
index 5a24546..8c69f82 100644
--- a/ML_BSSN/src/ML_BSSN_enforce.c
+++ b/ML_BSSN/src/ML_BSSN_enforce.c
@@ -1,5 +1,5 @@
-/*  File produced by user eschnett */
-/*  Produced with Mathematica Version 6.0 for Mac OS X x86 (32-bit) (April 20, 2007) */
+/*  File produced by user diener */
+/*  Produced with Mathematica Version 6.0 for Linux x86 (32-bit) (April 20, 2007) */
 
 /*  Mathematica script written by Ian Hinder and Sascha Husa */
 
@@ -99,6 +99,7 @@ void ML_BSSN_enforce_Body(cGH *cctkGH, CCTK_INT dir, CCTK_INT face, CCTK_REAL no
     CCTK_REAL trA = INITVALUE;
     
     /* Declare local copies of grid functions */
+    CCTK_REAL alphaL = INITVALUE;
     CCTK_REAL At11L = INITVALUE, At12L = INITVALUE, At13L = INITVALUE, At22L = INITVALUE, At23L = INITVALUE, At33L = INITVALUE;
     CCTK_REAL gt11L = INITVALUE, gt12L = INITVALUE, gt13L = INITVALUE, gt22L = INITVALUE, gt23L = INITVALUE, gt33L = INITVALUE;
     /* Declare precomputed derivatives*/
@@ -106,6 +107,7 @@ void ML_BSSN_enforce_Body(cGH *cctkGH, CCTK_INT dir, CCTK_INT face, CCTK_REAL no
     /* Declare derivatives */
     
     /* Assign local copies of grid functions */
+    alphaL = alpha[index];
     At11L = At11[index];
     At12L = At12[index];
     At13L = At13[index];
@@ -156,8 +158,11 @@ void ML_BSSN_enforce_Body(cGH *cctkGH, CCTK_INT dir, CCTK_INT face, CCTK_REAL no
     
     At33L  =  At33L - gt33L*kthird*trA;
     
+    alphaL  =  fmax(1.e-10,alphaL);
+    
     
     /* Copy local copies back to grid functions */
+    alpha[index] = alphaL;
     At11[index] = At11L;
     At12[index] = At12L;
     At13[index] = At13L;
diff --git a/ML_BSSN/src/ML_BSSN_matter.c b/ML_BSSN/src/ML_BSSN_matter.c
index 3834a1f..2da5337 100644
--- a/ML_BSSN/src/ML_BSSN_matter.c
+++ b/ML_BSSN/src/ML_BSSN_matter.c
@@ -1,5 +1,5 @@
-/*  File produced by user eschnett */
-/*  Produced with Mathematica Version 6.0 for Mac OS X x86 (32-bit) (April 20, 2007) */
+/*  File produced by user diener */
+/*  Produced with Mathematica Version 6.0 for Linux x86 (32-bit) (April 20, 2007) */
 
 /*  Mathematica script written by Ian Hinder and Sascha Husa */
 
diff --git a/ML_BSSN/src/ML_BSSN_matter_constraints.c b/ML_BSSN/src/ML_BSSN_matter_constraints.c
index b8b41a3..2043e8e 100644
--- a/ML_BSSN/src/ML_BSSN_matter_constraints.c
+++ b/ML_BSSN/src/ML_BSSN_matter_constraints.c
@@ -1,5 +1,5 @@
-/*  File produced by user eschnett */
-/*  Produced with Mathematica Version 6.0 for Mac OS X x86 (32-bit) (April 20, 2007) */
+/*  File produced by user diener */
+/*  Produced with Mathematica Version 6.0 for Linux x86 (32-bit) (April 20, 2007) */
 
 /*  Mathematica script written by Ian Hinder and Sascha Husa */
 
diff --git a/ML_BSSN/src/RegisterMoL.c b/ML_BSSN/src/RegisterMoL.c
index 0a5e78e..0bdc194 100644
--- a/ML_BSSN/src/RegisterMoL.c
+++ b/ML_BSSN/src/RegisterMoL.c
@@ -1,5 +1,5 @@
-/*  File produced by user eschnett */
-/*  Produced with Mathematica Version 6.0 for Mac OS X x86 (32-bit) (April 20, 2007) */
+/*  File produced by user diener */
+/*  Produced with Mathematica Version 6.0 for Linux x86 (32-bit) (April 20, 2007) */
 
 /*  Mathematica script written by Ian Hinder and Sascha Husa */
 
diff --git a/ML_BSSN/src/RegisterSymmetries.c b/ML_BSSN/src/RegisterSymmetries.c
index 5bcde48..44f11de 100644
--- a/ML_BSSN/src/RegisterSymmetries.c
+++ b/ML_BSSN/src/RegisterSymmetries.c
@@ -1,5 +1,5 @@
-/*  File produced by user eschnett */
-/*  Produced with Mathematica Version 6.0 for Mac OS X x86 (32-bit) (April 20, 2007) */
+/*  File produced by user diener */
+/*  Produced with Mathematica Version 6.0 for Linux x86 (32-bit) (April 20, 2007) */
 
 /*  Mathematica script written by Ian Hinder and Sascha Husa */
 
diff --git a/ML_BSSN/src/Startup.c b/ML_BSSN/src/Startup.c
index 2bb6a00..3871ecb 100644
--- a/ML_BSSN/src/Startup.c
+++ b/ML_BSSN/src/Startup.c
@@ -1,5 +1,5 @@
-/*  File produced by user eschnett */
-/*  Produced with Mathematica Version 6.0 for Mac OS X x86 (32-bit) (April 20, 2007) */
+/*  File produced by user diener */
+/*  Produced with Mathematica Version 6.0 for Linux x86 (32-bit) (April 20, 2007) */
 
 /*  Mathematica script written by Ian Hinder and Sascha Husa */
 
diff --git a/ML_BSSN/src/make.code.defn b/ML_BSSN/src/make.code.defn
index 2d1b2a0..42f1a1b 100644
--- a/ML_BSSN/src/make.code.defn
+++ b/ML_BSSN/src/make.code.defn
@@ -1,5 +1,5 @@
-# File produced by user eschnett
-# Produced with Mathematica Version 6.0 for Mac OS X x86 (32-bit) (April 20, 2007)
+# File produced by user diener
+# Produced with Mathematica Version 6.0 for Linux x86 (32-bit) (April 20, 2007)
 
 # Mathematica script written by Ian Hinder and Sascha Husa