From 50a5c30dd57a72d33a84517e78c36f5dd4123d67 Mon Sep 17 00:00:00 2001
From: Matthias Frey <matthias.frey@psi.ch>
Date: Wed, 15 Jul 2020 10:28:26 +0200
Subject: [PATCH] RectangularDomain: fix compute function

---
 src/Solvers/RectangularDomain.cpp | 4 ++++
 src/Solvers/RectangularDomain.h   | 4 +---
 2 files changed, 5 insertions(+), 3 deletions(-)

diff --git a/src/Solvers/RectangularDomain.cpp b/src/Solvers/RectangularDomain.cpp
index 741d3988c..a2be643a3 100644
--- a/src/Solvers/RectangularDomain.cpp
+++ b/src/Solvers/RectangularDomain.cpp
@@ -47,6 +47,10 @@ void RectangularDomain::compute(Vector_t hr){
     nxy_m = nr[0] * nr[1];
 }
 
+void RectangularDomain::compute(Vector_t hr, NDIndex<3> /*localId*/) {
+    compute(hr);
+}
+
 int RectangularDomain::getNumXY(int /*z*/) {
     return nxy_m;
 }
diff --git a/src/Solvers/RectangularDomain.h b/src/Solvers/RectangularDomain.h
index 9d4997d38..2fd853b92 100644
--- a/src/Solvers/RectangularDomain.h
+++ b/src/Solvers/RectangularDomain.h
@@ -44,9 +44,7 @@ public:
     /// calculates intersection with the beam pipe
     void compute(Vector_t hr);
 
-    void compute(Vector_t /*hr*/, NDIndex<3> /*localId*/) {
-        // we have no intersections
-    }
+    void compute(Vector_t hr, NDIndex<3> /*localId*/);
 
     /// returns number of nodes in xy plane (here independent of z coordinate)
     int getNumXY(int z);
-- 
GitLab