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<float> local_solution;
mg_matrices[level].initialize_dof_vector(local_solution);
// local_solution = solution;
mg_matrices[level].set_old_solution(local_solution);
in the loop while creating those matrices. My next problem, though, is how
to fill the local_solution-vector based on the solution-vector. When
uncommenting the line above, I get exactly the same error as before. Are
there alternatives?
Am Freitag, 26. Juli 2019 19:36:49 UTC+2 schrieb Daniel Arndt:
>
> 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:
>> The parallel layout of the given vector is not compatible with the
>> parallel partitioning in MatrixFree. Use MatrixFree::initialize_dof_vector
>> to get a compatible vector.
>>
>> Did you initialize its layout via MatrixFree::initialize_dof_vector()
> <http:///> in fact?
>
> Best,
> Daniel
>
--
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].
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/ae0ba1a3-d2bc-4d95-a5a6-1a89331a0c4f%40googlegroups.com.
