diff options
Diffstat (limited to 'ML_BSSN/src/ML_BSSN_RHS.cc')
| -rw-r--r-- | ML_BSSN/src/ML_BSSN_RHS.cc | 68 |
1 files changed, 30 insertions, 38 deletions
diff --git a/ML_BSSN/src/ML_BSSN_RHS.cc b/ML_BSSN/src/ML_BSSN_RHS.cc index a3bb92e..1152240 100644 --- a/ML_BSSN/src/ML_BSSN_RHS.cc +++ b/ML_BSSN/src/ML_BSSN_RHS.cc @@ -2518,55 +2518,47 @@ static void ML_BSSN_RHS_Body(const cGH* restrict const cctkGH, const int dir, co } #else { - double gu[3][3] = {{ gu11, gu12, gu13 }, - { gu12, gu22, gu23 }, - { gu13, gu23, gu33 }}; - double G[3][3][3] = {{{ G111, G112, G113 }, - { G112, G122, G123 }, - { G113, G123, G133 }}, - {{ G211, G212, G213 }, - { G212, G222, G223 }, - { G213, G223, G233 }}, - {{ G311, G312, G313 }, - { G312, G322, G323 }, - { G313, G323, G333 }}}; + // double gu[3][3] = {{ gu11, gu12, gu13 }, + // { gu12, gu22, gu23 }, + // { gu13, gu23, gu33 }}; + // double G[3][3][3] = {{{ G111, G112, G113 }, + // { G112, G122, G123 }, + // { G113, G123, G133 }}, + // {{ G211, G212, G213 }, + // { G212, G222, G223 }, + // { G213, G223, G233 }}, + // {{ G311, G312, G313 }, + // { G312, G322, G323 }, + // { G313, G323, G333 }}}; - double X1 = 0.0, X3 = 0.0; + // double X1 = 0.0, X3 = 0.0; - double dux_alpha, duz_alpha; + // double dux_alpha, duz_alpha; - for (int fn = 0; fn < 3; fn++) - for (int bz = 0; bz < 3; bz++) { - X1 += gu[fn][bz] * G[0][fn][bz]; - X3 += gu[fn][bz] * G[2][fn][bz]; - } + // for (int fn = 0; fn < 3; fn++) + // for (int bz = 0; bz < 3; bz++) { + // X1 += gu[fn][bz] * G[0][fn][bz]; + // X3 += gu[fn][bz] * G[2][fn][bz]; + // } - dux_alpha = gu[0][0] * PDstandardNth1alpha + gu[0][2] * PDstandardNth3alpha; - duz_alpha = gu[2][0] * PDstandardNth1alpha + gu[2][2] * PDstandardNth3alpha; + // dux_alpha = gu[0][0] * PDstandardNth1alpha + gu[0][2] * PDstandardNth3alpha; + // duz_alpha = gu[2][0] * PDstandardNth1alpha + gu[2][2] * PDstandardNth3alpha; - beta1rhsL = alphaL * alphaL * X1 - alphaL * dux_alpha - alphaL * BetaDriver * beta1L; - beta2rhsL = 0.0; - beta3rhsL = alphaL * alphaL * X3 - alphaL * duz_alpha - alphaL * BetaDriver * beta3L; + // beta1rhsL = alphaL * alphaL * X1 - alphaL * dux_alpha - alphaL * BetaDriver * beta1L; + // beta2rhsL = 0.0; + // beta3rhsL = alphaL * alphaL * X3 - alphaL * duz_alpha - alphaL * BetaDriver * beta3L; } #endif else { - double xx = x[index]; - double zz = z[index]; - double t = cctkGH->cctk_time; - //double rr = sqrt(xx * xx + zz * zz); - //double phi = atan2(zz, xx); + beta1rhsL = theta * ToReal(ShiftGammaCoeff) * + (ToReal(ShiftBCoeff) * B1L + (kone - ToReal(ShiftBCoeff)) * (Xt1L - eta * ToReal(BetaDriver) * beta1L)); - //double gauge_r = -20.0 * (1.0 + 0.0 * cos(2 * phi) * rr * rr) * rr / (1.0 + rr * rr * rr * rr); + beta2rhsL = theta * ToReal(ShiftGammaCoeff) * + (ToReal(ShiftBCoeff) * B2L + (kone - ToReal(ShiftBCoeff)) * (Xt2L - eta * ToReal(BetaDriver) * beta2L)); - beta1rhsL = theta * ShiftGammaCoeff * - (ShiftBCoeff * B1L + (1.0 - ShiftBCoeff) * (Xt1L - eta * BetaDriver * beta1L)); - - beta2rhsL = theta * ShiftGammaCoeff * - (ShiftBCoeff * B2L + (1.0 - ShiftBCoeff) * (Xt2L - eta * BetaDriver * beta2L)); - - beta3rhsL = theta * ShiftGammaCoeff * - (ShiftBCoeff * B3L + (1.0 - ShiftBCoeff) * (Xt3L - eta * BetaDriver * beta3L)); + beta3rhsL = theta * ToReal(ShiftGammaCoeff) * + (ToReal(ShiftBCoeff) * B3L + (kone - ToReal(ShiftBCoeff)) * (Xt3L - eta * ToReal(BetaDriver) * beta3L)); } CCTK_REAL_VEC B1rhsL CCTK_ATTRIBUTE_UNUSED = |