From 946a27631f65631e0472ae95a1fe57dc4695e4dc Mon Sep 17 00:00:00 2001 From: Peter Diener Date: Thu, 12 Nov 2009 18:15:08 -0600 Subject: Add support for the W-method. Choose between phi and W-methods by setting the integer parameter conformalmethod to either 0 or 1. Signed-off-by: Peter Diener --- ML_BSSN_MP/src/ML_BSSN_MP_enforce.c | 10 ++++------ 1 file changed, 4 insertions(+), 6 deletions(-) (limited to 'ML_BSSN_MP/src/ML_BSSN_MP_enforce.c') diff --git a/ML_BSSN_MP/src/ML_BSSN_MP_enforce.c b/ML_BSSN_MP/src/ML_BSSN_MP_enforce.c index 3b8364e..b33a381 100644 --- a/ML_BSSN_MP/src/ML_BSSN_MP_enforce.c +++ b/ML_BSSN_MP/src/ML_BSSN_MP_enforce.c @@ -1,4 +1,7 @@ -/* File produced by Kranc */ +/* File produced by user diener */ +/* Produced with Mathematica Version 7.0 for Linux x86 (64-bit) (February 18, 2009) */ + +/* Mathematica script written by Ian Hinder and Sascha Husa */ #define KRANC_C @@ -106,7 +109,6 @@ void ML_BSSN_MP_enforce_Body(cGH const * const cctkGH, CCTK_INT const dir, CCTK_ CCTK_REAL trAt = 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*/ @@ -114,7 +116,6 @@ void ML_BSSN_MP_enforce_Body(cGH const * const cctkGH, CCTK_INT const dir, CCTK_ /* Declare derivatives */ /* Assign local copies of grid functions */ - alphaL = alpha[index]; At11L = At11[index]; At12L = At12[index]; At13L = At13[index]; @@ -165,11 +166,8 @@ void ML_BSSN_MP_enforce_Body(cGH const * const cctkGH, CCTK_INT const dir, CCTK_ At33L = At33L - gt33L*kthird*trAt; - alphaL = fmax(alphaL,MinimumLapse); - /* Copy local copies back to grid functions */ - alpha[index] = alphaL; At11[index] = At11L; At12[index] = At12L; At13[index] = At13L; -- cgit v1.2.3