From 9faa628ffb96d1f5e9e2a97e0e150f7534ce7003 Mon Sep 17 00:00:00 2001
From: Anton Khirnov <anton@khirnov.net>
Date: Sun, 2 Jun 2019 17:42:45 +0200
Subject: Add 1st order dissipation.

---
 src/apply_dissipation.c | 10 ++++++++++
 1 file changed, 10 insertions(+)

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++)
-- 
cgit v1.2.3