revert a part of the code to a previous status in elliptic interface …

…laws
pmgbergen · Feb 21, 2019 · 6ea62d6 · 6ea62d6
1 parent e2bbc34
commit 6ea62d6
Showing 1 changed file with 6 additions and 3 deletions.
diff --git a/src/porepy/numerics/interface_laws/elliptic_interface_laws.py b/src/porepy/numerics/interface_laws/elliptic_interface_laws.py
@@ -715,9 +715,12 @@ def assemble_matrix_rhs(
         self.discr_master.enforce_neumann_int_bound(
             g_master, data_edge, matrix, False, master_ind
         )
-        self.discr_slave.enforce_neumann_int_bound(
-            g_slave, data_edge, matrix, True, slave_ind
-        )
+
+        # Consider this terms only if the grids are of the same dimension
+        if g_master.dim == g_slave.dim:
+            self.discr_slave.enforce_neumann_int_bound(
+                g_slave, data_edge, matrix, True, slave_ind
+            )
 
         return matrix, rhs