From de98fa0710b68754536badf726661621ef5f2ab3 Mon Sep 17 00:00:00 2001 From: Martin Kronbichler Date: Tue, 24 Oct 2017 12:36:25 +0200 Subject: [PATCH] Avoid VectorView in matrix-free operators. --- include/deal.II/matrix_free/operators.h | 9 ++------- 1 file changed, 2 insertions(+), 7 deletions(-) diff --git a/include/deal.II/matrix_free/operators.h b/include/deal.II/matrix_free/operators.h index c72ec5a5c4..cc65ed759c 100644 --- a/include/deal.II/matrix_free/operators.h +++ b/include/deal.II/matrix_free/operators.h @@ -24,7 +24,6 @@ #include #include -#include #include #include #include @@ -1197,14 +1196,10 @@ namespace MatrixFreeOperators // copy the vector content to a temporary vector so that it does not get // lost - VectorView view_src_in(subblock(src,i).local_size(), - subblock(src,i).begin()); - const Vector copy_vec = view_src_in; + LinearAlgebra::distributed::Vector copy_vec(subblock(src,i)); subblock(const_cast(src),i). reinit(data->get_dof_info(mf_component).vector_partitioner); - VectorView view_src_out(subblock(src,i).local_size(), - subblock(src,i).begin()); - static_cast&>(view_src_out) = copy_vec; + subblock(const_cast(src),i).copy_locally_owned_data_from(copy_vec); } } -- 2.39.5