From c0ff199014ea97969744df5a1fedc5c2bc70ae96 Mon Sep 17 00:00:00 2001
From: Anton Khirnov <anton@khirnov.net>
Date: Tue, 7 Aug 2018 16:22:37 +0200
Subject: ell_relax: fill in boundary ghosts for fixval boundaries as well

---
 ell_relax.c | 6 +++++-
 1 file changed, 5 insertions(+), 1 deletion(-)

(limited to 'ell_relax.c')

diff --git a/ell_relax.c b/ell_relax.c
index e0eadd4..679f340 100644
--- a/ell_relax.c
+++ b/ell_relax.c
@@ -118,8 +118,12 @@ static void boundaries_apply_fixval(double *dst, ptrdiff_t dst_stride,
                                     size_t boundary_size, ptrdiff_t boundary_stride,
                                     int is_upper)
 {
+    if (!is_upper)
+        boundary_stride *= -1;
+
     for (size_t i = 0; i < boundary_size; i++) {
-        *dst = *src;
+        for (int j = 0; j <= FD_STENCIL_MAX; j++)
+            dst[j * boundary_stride] = *src;
         dst += dst_stride;
         src++;
     }
-- 
cgit v1.2.3