diff options
author | Anton Khirnov <anton@khirnov.net> | 2020-08-12 18:09:58 +0200 |
---|---|---|
committer | Anton Khirnov <anton@khirnov.net> | 2020-08-12 18:09:58 +0200 |
commit | 3da7478f9c5852563ee706b991e40ce4b6b6656f (patch) | |
tree | 5443b67a39becac360f267ccea7ee0d23c05ae7b | |
parent | 779de67ad5a8ffc90fd5b908941a492461473ebc (diff) |
Vectorize BSSN
-rw-r--r-- | ML_BSSN/src/ML_BSSN_lapse_evol.cc | 2 |
1 files changed, 1 insertions, 1 deletions
diff --git a/ML_BSSN/src/ML_BSSN_lapse_evol.cc b/ML_BSSN/src/ML_BSSN_lapse_evol.cc index 0388510..3cef863 100644 --- a/ML_BSSN/src/ML_BSSN_lapse_evol.cc +++ b/ML_BSSN/src/ML_BSSN_lapse_evol.cc @@ -219,7 +219,7 @@ static void ML_BSSN_lapse_evol_Body(const cGH* restrict const cctkGH, const int #if 1 if (lapse_gh) { - CCTK_REAL_VEC e4phi = kifthen(conformalMethod == 1, ToReal(1.0) / (phiL * phiL), kexp(ToReal(4.0) * phiL)); + CCTK_REAL_VEC e4phi = conformalMethod == 1 ? ToReal(1.0) / (phiL * phiL) : kexp(ToReal(4.0) * phiL); CCTK_REAL_VEC det_gamma = kpow(e4phi, 3.0); alpharhsL = -(alphaL * alphaL) * trKL + ToReal(gh_eta_l_bar) * kpow(alphaL, gh_q) * (alphaL * alphaL) * klog(kpow(det_gamma, 0.5 * gh_p) / alphaL) |