Thanks a lot to both of you! First, I modified my code in a way that the "colored iterators" (still one color) only contains locally owned cells. Unfortunately, this didn't help to avoid the crash and I still got the same random error as before.
Next, I tested Martin's suggestion: > The solution to this problem is to either initialize the vector with the > appropriate ghost entries (locally relevant dofs), i.e., the reinit > method/constructor that takes two index sets, an MPI_Comm and a > vector_writable flag. Alternatively, you can use > parallel::distributed::Vector that you initialize with the same two index > sets, locally owned dofs as first argument and locally relevant dofs as the > second element (i.e., ghosts are the locally relevant minus the locally > owned dofs). Trilinos matrices and preconditioners should know how to > collaborate with these vectors. And it works fine. To be clear, instead of using dealii::TrilinosWrappers::MPI::Vector dst(locally_owned_dofs,MPI_COMM_WORLD ); the correct way would be dealii::TrilinosWrappers::MPI::Vector dst; dst.reinit(locally_owned_dofs,locally_relevant_dofs,MPI_COMM_WORLD,true); This attempt also works fine in a much more sophisticated code of mine. And maybe just to note it, it runs for both cases, using the whole range of iterators as well as having a range consisting only of locally owned iterators. In my eyes this is fulfilled by the fact that cell_action() checks, if the given cell is locally owned or not. Again, thanks a lot for all your kind responses and your offered help. Best, Julius P.S. I will upload the corrected code after the holidays as it is on my local machine in the university. -- 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.
