Merge pull request #4900 from davydden/mf_operator_fix

fix a minor bug in MatrixFreeOperators::Base::vmult_interface_down()
dealii · Aug 18, 2017 · 5a5d769 · 5a5d769
2 parents 3057d72 + 3683cf3
commit 5a5d769
Showing 1 changed file with 1 addition and 1 deletion.
diff --git 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.;
       }
   }