From 86ad823b9ade211bfa9361b61571933aff1c9d24 Mon Sep 17 00:00:00 2001
From: Anton Khirnov <anton@khirnov.net>
Date: Tue, 23 Apr 2019 15:36:34 +0200
Subject: egs: merge residual calc and correct when possible

Also, merge the reflect boundary condition into residual calc+add.
Improves performance due to better locality.
---
 relax_test.c | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

(limited to 'relax_test.c')

diff --git a/relax_test.c b/relax_test.c
index 5c07107..50ad04c 100644
--- a/relax_test.c
+++ b/relax_test.c
@@ -144,7 +144,7 @@ int main(int argc, char **argv)
     res_old = findmax(ctx->residual->data, ctx->residual->stride[0] * ctx->domain_size[1]);
 
     for (int i = 0; i < maxiter; i++) {
-        ret = mg2di_egs_solve(ctx);
+        ret = mg2di_egs_solve(ctx, 0);
         if (ret < 0) {
             fprintf(stderr, "Error during relaxation\n");
             ret = 1;
-- 
cgit v1.2.3