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 replaced, and therefore making it difficult to implement it in my problem. Thus, what would be the most efficient way of giving the operator-function apply_add access to those solution vectors? 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.
Are there better ways? Thanks! -- 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/855a446f-4bde-4a56-8ae2-1871488f0b65%40googlegroups.com.