The error is from compute_component_wise on level = 0, the result it returns
don't equal to dof_handler.n_dofs(level).
I added a MPI_Barrier between first and second loop but error is not eliminated.
Snipaste_2022-11-28_12-53-02.png
I do not see this error. Are you running in debug mode? And with how many MPI
processes?
As a side note: The detour by writing the data to a file and then reading it
again is very inefficient, in particular if you have many MPI processes. If
you want to use GridOut/GridIn, let the former write the data into an object
of type std::ostringstream, then use Utilities::MPI::broadcast() to send the
string you get out of std::ostringstream to the other processes, and give it
to a std::istringstream object. Then use GridIn on the std::istringstream
object. This way, no data transfer to or from disk is necessary -- it all
happens in memory.
Best
Wolfgang
--
------------------------------------------------------------------------
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.
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/35e96ca4-0b33-1870-cde9-2e42c94eddf4%40colostate.edu.