From: Martin Kronbichler Date: Tue, 24 Oct 2017 10:36:25 +0000 (+0200) Subject: Avoid VectorView in matrix-free operators. X-Git-Tag: v9.0.0-rc1~889^2~3 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=de98fa0710b68754536badf726661621ef5f2ab3;p=dealii.git Avoid VectorView in matrix-free operators. --- 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); } }