RectangularDomain: fix compute function

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;
 }

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);