fix a minor bug in MatrixFreeOperators::Base::vmult_interface_down()

diff --git a/include/deal.II/matrix_free/operators.h b/include/deal.II/matrix_free/operators.h
index 660aefad1af564ab11251a543cf2ed54d32f0bcc..be1293e32adea2ca0d1c1e9fd0799b6a853e14e3 100644 (file)
--- a/include/deal.II/matrix_free/operators.h
+++ b/include/deal.II/matrix_free/operators.h
@@ -1312,7 +1312,7 @@ namespace MatrixFreeOperators
             .local_element(edge_constrained_indices[j][i])
               = edge_constrained_values[j][i].first;
           }
-        for ( ; c<dst.local_size(); ++c)
+        for ( ; c<subblock(dst,j).local_size(); ++c)
           subblock(dst,j).local_element(c) = 0.;
       }
   }