Hi Denis,
Wouldn't it be easier to just create a copy of the vector by
ghosted_vector.reinit(locally_owned_set, locally_relevant_set, comm);
ghosted_vector = solution;
ghosted_vector.update_ghost_values();
The reason why there is a VectorView in the MatrixFreeOperators::Base
class is that I want to modify the original vector which I don't assume
you need in your code. (It would actually be dangerous if you used a
vector with the locally relevant dofs as ghosts because that gets the
local indices wrong as compared to the storage of MatrixFree-compatible
vectors.)
Best,
Martin
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups "deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
For more options, visit https://groups.google.com/d/optout.
