diff options
author | Peter Diener <diener@linux-hn3d.site> | 2009-11-12 18:15:08 -0600 |
---|---|---|
committer | Peter Diener <diener@linux-hn3d.site> | 2009-11-12 18:15:08 -0600 |
commit | 946a27631f65631e0472ae95a1fe57dc4695e4dc (patch) | |
tree | 47c53575bb8bd4554208bb27515322d5e51f3a4f /ML_BSSN_MP/src/ML_BSSN_MP_enforce.c | |
parent | a52abf35db4c68c5e28e8ce9b00b394c4ea382ee (diff) |
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 <diener@linux-hn3d.site>
Diffstat (limited to 'ML_BSSN_MP/src/ML_BSSN_MP_enforce.c')
-rw-r--r-- | ML_BSSN_MP/src/ML_BSSN_MP_enforce.c | 10 |
1 files changed, 4 insertions, 6 deletions
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; |