I found my mistake, I could not reuse the same partition for the main
matrix and the MG-matrices, thus I had to add
LinearAlgebra::distributed::Vector local_solution;
mg_matrices[level].initialize_dof_vector(local_solution);
//local_solution = solution;
Maxi,
[...] I tried adding a class member and initializing it during
> setup_system(), but when accessing it via read_dof_values(), I had a layout
> mismatch, resulting in
> The violated condition was:
> vec.partitioners_are_compatible(*dof_info.vector_partitioner)
> Additional information:
>
I followed step-37 for creating a MatrixFreeOperator, but now I also would
like to access a previous solution (which is needed during apply_add()),
and stored using a LinearAlgebra::distributed::Vector. In step-48 there is
an example which uses previous solutions, but here the solve-function is