Dear all, I am having some trouble understanding how to have access to the different components (2 for complex values) in a matrix-free implementation and MPI. I need to have access to each component to perform the different operations to them in the *local_apply()* function.
I saw this documentation (https://www.dealii.org/current/doxygen/deal.II/classFEEvaluation.html, @Handling multi-component systems) which get the components using BlockVectors. However, I define in my class *FESystem<dim> fe; *and after, *fe(FE_Q<dim>(degree), 2) * to have 2 components and use LinearAlgebra::distributed::Vector<number> <https://www.dealii.org/current/doxygen/deal.II/classLinearAlgebra_1_1distributed_1_1Vector.html> for the solution, dst, scr. In this case, is it possible to have access to the different components using blocks() as it is done in the example ( https://github.com/dealii/dealii/blob/master/tests/matrix_free/matrix_vector_stokes_noflux.cc, line 67) using for example phi.read_dof_values(src.block(0)) ? I would really appreciate it if I can get help on this matter. Regards, H. -- 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/320fbf7d-9532-4c2f-9952-46e3f501316en%40googlegroups.com.
