diff options
| author | Peter Diener <diener@cct.lsu.edu> | 2012-10-23 10:58:16 -0500 |
|---|---|---|
| committer | Peter Diener <diener@cct.lsu.edu> | 2012-10-23 10:58:16 -0500 |
| commit | 129cc3ba1a89b8336ffa0355fc74c2f1b60e8255 (patch) | |
| tree | 92e7da8e6ee4357d05cc4f8fc46037c3450e4f5f /ML_BSSN/src/ML_BSSN_enforce.cc | |
| parent | 577bb4833afad4286c8db76c6327e89c0ae0ea17 (diff) | |
Change the dafault value of apply_dissipation to "never".
Change the default value of apply_dissipation from "always" to "never".
Since the default dissipation strength is zero this should not change
results for users using default values but should save time.
Signed-off-by: Peter Diener <diener@cct.lsu.edu>
Diffstat (limited to 'ML_BSSN/src/ML_BSSN_enforce.cc')
| -rw-r--r-- | ML_BSSN/src/ML_BSSN_enforce.cc | 16 |
1 files changed, 8 insertions, 8 deletions
diff --git a/ML_BSSN/src/ML_BSSN_enforce.cc b/ML_BSSN/src/ML_BSSN_enforce.cc index d63111a..3e46718 100644 --- a/ML_BSSN/src/ML_BSSN_enforce.cc +++ b/ML_BSSN/src/ML_BSSN_enforce.cc @@ -117,12 +117,12 @@ static void ML_BSSN_enforce_Body(cGH const * restrict const cctkGH, int const di CCTK_REAL_VEC const p1o840dx = kdiv(ToReal(0.00119047619047619047619047619048),dx); CCTK_REAL_VEC const p1o840dy = kdiv(ToReal(0.00119047619047619047619047619048),dy); CCTK_REAL_VEC const p1o840dz = kdiv(ToReal(0.00119047619047619047619047619048),dz); - CCTK_REAL_VEC const p1odx = kdiv(ToReal(1.),dx); - CCTK_REAL_VEC const p1odx2 = kdiv(ToReal(1.),kmul(dx,dx)); - CCTK_REAL_VEC const p1ody = kdiv(ToReal(1.),dy); - CCTK_REAL_VEC const p1ody2 = kdiv(ToReal(1.),kmul(dy,dy)); - CCTK_REAL_VEC const p1odz = kdiv(ToReal(1.),dz); - CCTK_REAL_VEC const p1odz2 = kdiv(ToReal(1.),kmul(dz,dz)); + CCTK_REAL_VEC const p1odx = kdiv(ToReal(1),dx); + CCTK_REAL_VEC const p1odx2 = kdiv(ToReal(1),kmul(dx,dx)); + CCTK_REAL_VEC const p1ody = kdiv(ToReal(1),dy); + CCTK_REAL_VEC const p1ody2 = kdiv(ToReal(1),kmul(dy,dy)); + CCTK_REAL_VEC const p1odz = kdiv(ToReal(1),dz); + CCTK_REAL_VEC const p1odz2 = kdiv(ToReal(1),kmul(dz,dz)); CCTK_REAL_VEC const pm1o120dx = kdiv(ToReal(-0.00833333333333333333333333333333),dx); CCTK_REAL_VEC const pm1o120dy = kdiv(ToReal(-0.00833333333333333333333333333333),dy); CCTK_REAL_VEC const pm1o120dz = kdiv(ToReal(-0.00833333333333333333333333333333),dz); @@ -239,7 +239,7 @@ static void ML_BSSN_enforce_Body(cGH const * restrict const cctkGH, int const di } /* Calculate temporaries and grid functions */ - CCTK_REAL_VEC detgt = ToReal(1.); + CCTK_REAL_VEC detgt = ToReal(1); CCTK_REAL_VEC gtu11 = kdiv(kmsub(gt22L,gt33L,kmul(gt23L,gt23L)),detgt); @@ -260,7 +260,7 @@ static void ML_BSSN_enforce_Body(cGH const * restrict const cctkGH, int const di kdiv(kmsub(gt11L,gt22L,kmul(gt12L,gt12L)),detgt); CCTK_REAL_VEC trAt = - kmadd(At11L,gtu11,kmadd(At22L,gtu22,kmadd(At33L,gtu33,kmul(kmadd(At12L,gtu12,kmadd(At13L,gtu13,kmul(At23L,gtu23))),ToReal(2.))))); + kmadd(At11L,gtu11,kmadd(At22L,gtu22,kmadd(At33L,gtu33,kmul(kmadd(At12L,gtu12,kmadd(At13L,gtu13,kmul(At23L,gtu23))),ToReal(2))))); At11L = kmadd(gt11L,kmul(trAt,ToReal(-0.333333333333333333333333333333)),At11L); |