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;
mg_mat
>
> Yeah, I was testing it in release mode. After running in debug mode I
> realized i didn't have hypre installed. Now it works.
>
Thanks
Prakarsh
--
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
-
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:
>
Franco,
If you forget about the static_cast in the error message just above
>> you can actually see what happens:
>> The violated condition was:
>> global_matrix.m() == global_vector.size()
>>
>>
>
> I am lost in this, I have made some stupid error here but I don't see how.
> I am posting th
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
> No unfortunately we don't have such things. However, most of the new
> effort on linear algebra is in the LinearAlgebra namespace.
>
>
Thanks, Bruno, I've started looking into that namespace. There are several
options to make it work in a coherent way, I'll ask for comments on these
soluti