diff options
author | Anton Khirnov <anton@khirnov.net> | 2019-06-02 17:42:45 +0200 |
---|---|---|
committer | Anton Khirnov <anton@khirnov.net> | 2019-06-02 17:42:45 +0200 |
commit | 9faa628ffb96d1f5e9e2a97e0e150f7534ce7003 (patch) | |
tree | f538e93e512ebf75568f60bd2a821488b5823f3c | |
parent | b212d1e5048f48d3e5146ccdac8386f075e8aee7 (diff) |
Add 1st order dissipation.
-rw-r--r-- | src/apply_dissipation.c | 10 |
1 files changed, 10 insertions, 0 deletions
diff --git a/src/apply_dissipation.c b/src/apply_dissipation.c index ce6e223..f550223 100644 --- a/src/apply_dissipation.c +++ b/src/apply_dissipation.c @@ -23,6 +23,16 @@ void apply_dissipation(const cGH *gh, const double *var, double *rhs, double dx[ } switch (order) { + case 1: +#pragma omp parallel for + for (int k = 1; k < nk - 1; k++) + for (int i = 1; i < nk - 1; i++) { + int idx = CCTK_GFINDEX3D(gh, i, y_idx, k); + rhs[idx] += epsdis[idx] * + ((var[idx - 1] - 2.0 * var[idx] + var[idx + 1]) * inv_dx[0] + + (var[idx - 1 * ni * nj] - 2.0 * var[idx] + var[idx + 1 * ni * nj]) * inv_dx[2]); + } + break; case 3: #pragma omp parallel for for (int k = 2; k < nk - 2; k++) |