diff options
Diffstat (limited to 'ML_BSSN/src/ML_BSSN_lapse_evol.cc')
| -rw-r--r-- | ML_BSSN/src/ML_BSSN_lapse_evol.cc | 16 |
1 files changed, 13 insertions, 3 deletions
diff --git a/ML_BSSN/src/ML_BSSN_lapse_evol.cc b/ML_BSSN/src/ML_BSSN_lapse_evol.cc index 040cbb4..988b586 100644 --- a/ML_BSSN/src/ML_BSSN_lapse_evol.cc +++ b/ML_BSSN/src/ML_BSSN_lapse_evol.cc @@ -159,6 +159,7 @@ static void ML_BSSN_lapse_evol_Body(const cGH* restrict const cctkGH, const int CCTK_REAL_VEC beta1L CCTK_ATTRIBUTE_UNUSED = vec_load(beta1[index]); CCTK_REAL_VEC beta2L CCTK_ATTRIBUTE_UNUSED = vec_load(beta2[index]); CCTK_REAL_VEC beta3L CCTK_ATTRIBUTE_UNUSED = vec_load(beta3[index]); + CCTK_REAL_VEC phiL CCTK_ATTRIBUTE_UNUSED = vec_load(phi[index]); CCTK_REAL_VEC PDupwindNthAnti1alpha CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDupwindNthSymm1alpha CCTK_ATTRIBUTE_UNUSED; @@ -213,9 +214,18 @@ static void ML_BSSN_lapse_evol_Body(const cGH* restrict const cctkGH, const int PDupwindNthAnti2alpha = 0.0; PDupwindNthSymm2alpha = 0.0; - alpharhs[index] = -harmonicF * pow(alphaL, harmonicN) * trKL + WFactor * WL - + LapseAdvectionCoeff * (beta1L * PDupwindNthAnti1alpha + beta3L * PDupwindNthAnti3alpha + - fabs(beta1L) * PDupwindNthSymm1alpha + fabs(beta3L) * PDupwindNthSymm3alpha); + if (lapse_gh) { + const double e4phi = conformalMethod == 1 ? 1.0 / (phiL * phiL) : exp(4.0 * phiL); + const double det_gamma = pow(e4phi, 3.0); + + alpharhs[index] = -(alphaL * alphaL) * trKL + gh_eta_l_bar * pow(alphaL, gh_q) * (alphaL * alphaL) * log(pow(det_gamma, 0.5 * gh_p) / alphaL) + + LapseAdvectionCoeff * (beta1L * PDupwindNthAnti1alpha + beta3L * PDupwindNthAnti3alpha + + fabs(beta1L) * PDupwindNthSymm1alpha + fabs(beta3L) * PDupwindNthSymm3alpha); + } else { + alpharhs[index] = -harmonicF * pow(alphaL, harmonicN) * trKL + WFactor * WL + + LapseAdvectionCoeff * (beta1L * PDupwindNthAnti1alpha + beta3L * PDupwindNthAnti3alpha + + fabs(beta1L) * PDupwindNthSymm1alpha + fabs(beta3L) * PDupwindNthSymm3alpha); + } } CCTK_ENDLOOP3STR(ML_BSSN_lapse_evol); |