diff options
| author | Peter Diener <diener@diener-3.lsu.edu> | 2008-02-05 10:44:41 -0600 |
|---|---|---|
| committer | Peter Diener <diener@diener-3.lsu.edu> | 2008-02-05 10:44:41 -0600 |
| commit | 3e597749893c2ae7b5025801e70ea7a2efc752fa (patch) | |
| tree | fea30b3d2cb4ced32b1564da593f962af0a47aa9 /ML_BSSN/src | |
| parent | f6fc18f1a4087b106e84f3113b5c8af3e1314ac1 (diff) | |
Don't enforce lapse>0.
Comment out enforcement of lapse>0. Doesn't seem to be necessary and it
fails on santaka (fmax seems to be bad).
Signed-off-by: Peter Diener <diener@diener-3.lsu.edu>
Diffstat (limited to 'ML_BSSN/src')
| -rw-r--r-- | ML_BSSN/src/ML_BSSN_enforce.c | 5 |
1 files changed, 0 insertions, 5 deletions
diff --git a/ML_BSSN/src/ML_BSSN_enforce.c b/ML_BSSN/src/ML_BSSN_enforce.c index 57b9ed6..792adb4 100644 --- a/ML_BSSN/src/ML_BSSN_enforce.c +++ b/ML_BSSN/src/ML_BSSN_enforce.c @@ -100,7 +100,6 @@ void ML_BSSN_enforce_Body(cGH *cctkGH, CCTK_INT dir, CCTK_INT face, CCTK_REAL no 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*/ @@ -108,7 +107,6 @@ 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]; @@ -159,11 +157,8 @@ void ML_BSSN_enforce_Body(cGH *cctkGH, CCTK_INT dir, CCTK_INT face, CCTK_REAL no At33L = At33L - gt33L*kthird*trAt; - 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; |