diff options
| author | Erik Schnetter <schnetter@gmail.com> | 2012-01-04 12:20:56 -0500 |
|---|---|---|
| committer | Erik Schnetter <schnetter@gmail.com> | 2012-01-04 12:20:56 -0500 |
| commit | 8cb581f64d8116593bbc0b046571ac17e932b612 (patch) | |
| tree | e501909be38ff04f5debef51e0634bd49a6b5acc /m | |
| parent | 76afdc03ad5341cd043b5740f06fe0a30388fc00 (diff) | |
Apply advection terms to A or B^i only if they are evolved
Diffstat (limited to 'm')
| -rw-r--r-- | m/McLachlan_BSSN.m | 13 |
1 files changed, 7 insertions, 6 deletions
diff --git a/m/McLachlan_BSSN.m b/m/McLachlan_BSSN.m index 4df6573..a8f6984 100644 --- a/m/McLachlan_BSSN.m +++ b/m/McLachlan_BSSN.m @@ -705,17 +705,18 @@ advectCalc = + LapseAdvectionCoeff Upwind[beta[ua], alpha, la], dot[A] -> dot[A] - + LapseAdvectionCoeff Upwind[beta[ua], A, la] - + (LapseACoeff - LapseAdvectionCoeff) - Upwind[beta[ua], trK, la], + + LapseACoeff ( + + LapseAdvectionCoeff Upwind[beta[ua], A, la] + + (1 - LapseAdvectionCoeff) Upwind[beta[ua], trK, la]), dot[beta[ua]] -> dot[beta[ua]] + ShiftAdvectionCoeff Upwind[beta[ub], beta[ua], lb], dot[B[ua]] -> dot[B[ua]] - + ShiftAdvectionCoeff Upwind[beta[ub], B[ua], lb] - + (ShiftBCoeff - ShiftAdvectionCoeff) - Upwind[beta[ub], Xt[ua], lb] + + ShiftBCoeff ( + + ShiftAdvectionCoeff Upwind[beta[ub], B[ua], lb] + + ((1 - ShiftAdvectionCoeff) + Upwind[beta[ub], Xt[ua], lb])) } }; |