aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorAnton Khirnov <anton@khirnov.net>2020-06-26 15:08:25 +0200
committerAnton Khirnov <anton@khirnov.net>2020-06-26 15:08:25 +0200
commit073fc3f8982842a5dc91438fb9b351afff6b179e (patch)
tree9090b9cc2011382b35b40ba3bd7f7207e1431001
parent9faa628ffb96d1f5e9e2a97e0e150f7534ce7003 (diff)
apply_dissipation: fix sign for 3rd order dissipationcartoon
-rw-r--r--src/apply_dissipation.c2
1 files changed, 1 insertions, 1 deletions
diff --git a/src/apply_dissipation.c b/src/apply_dissipation.c
index f550223..e5ff6ec 100644
--- a/src/apply_dissipation.c
+++ b/src/apply_dissipation.c
@@ -38,7 +38,7 @@ void apply_dissipation(const cGH *gh, const double *var, double *rhs, double dx[
for (int k = 2; k < nk - 2; k++)
for (int i = 2; i < nk - 2; i++) {
int idx = CCTK_GFINDEX3D(gh, i, y_idx, k);
- rhs[idx] += epsdis[idx] / 16.0 *
+ rhs[idx] -= epsdis[idx] / 16.0 *
((var[idx - 2] - 4.0 * var[idx - 1] + 6.0 * var[idx] - 4.0 * var[idx + 1] + var[idx + 2]) * inv_dx[0] +
(var[idx - 2 * ni * nj] - 4.0 * var[idx - 1 * ni * nj] + 6.0 * var[idx] - 4.0 * var[idx + 1 * ni * nj] + var[idx + 2 * ni * nj]) * inv_dx[2]);
}