Anna,
Yes I would like to create a LA::MPI::BlockVector solution_at_receiver. And
then copy it to a
BlockVector localized_solution_at_receiver(solution_at_receiver)
And output localized_solution_at_receiver only on the 0 processor.
I am having trouble initializing and filling solution_at_receiver.
If I am simply write
solution_at_receiver.reinit(4);
for (unsigned int i=0; i<4; ++i)
solution_at_receiver.block(i).reinit(n_Rc*dim);
I am getting error message, that there is no such function reinit with such
arguments.
Well, this doesn't work because all you say that this (parallel!) should
have so and so many elements, but you don't say (i) which MPI
communicator the vector should used, nor (ii) which process is going to
own which elements of each of these vectors.
I think you may want to take a look at how something like this is
implemented in a different program:
https://github.com/geodynamics/aspect/blob/master/source/postprocess/point_values.cc#L34
Best
W.
--
------------------------------------------------------------------------
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
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 dealii+unsubscr...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.
