Well I renumbered the DoFs because this discussion <https://groups.google.com/g/dealii/c/MXQO5Nhcpws/m/kteyeUOcCAAJ#:~:text=II%20User%20Group-,Well%20it%20takes%20several%20steps%20to%20Renumber%20on%20the%20coarsest%20level.,-1)%20For%20a> . I use trilinos because MatrixFree method relys on LinearAlbegra::distributed Vector or Block Vector. When solve the coarsest level of MatrixFree multigrid problem directly, trilinos direct solver only handles non-blocked vectors, so I do transfer between BlockVectors and Vectors. The ordering of DoFs is destructed when fe degree rises to 3 and higher, causing the transfer error.
Best, Wayne 在2023年2月24日星期五 UTC+8 19:54:45<Wolfgang Bangerth> 写道： > On 2/24/23 04:37, 'yy.wayne' via deal.II User Group wrote: > > > > However another problem pops out. When test a vector-valued problem, > renumber > > DoFs > > component wise make reinit of vectors fail. If > DoFRenumbering::component_wise > > is called, > > > > 1. LinearAlgebra::distributed::Vector<double> > > v1(dof_handler.locally_owned_dofs(), mpi_communicator); > > v1.reinit(dof_handler.locally_owned_dofs(), mpi_communicator); > > Success > > 2. LinearAlgebra::distributed::Vector<double> v2; > > v2.reinit(dof_handler.locally_owned_dofs(), mpi_communicator); > > Fails. > > > > Is that relates to trilinos date type? > > If you do not arrange the different components of a vector-valued finite > element into a block vector (i.e., if you use a single vector for all > components of the solution), then there is no need to call > DoFRenumbering::component_wise(). > > Best > W. > > -- > ------------------------------------------------------------------------ > Wolfgang Bangerth email: [email protected] > www: http://www.math.colostate.edu/~bangerth/ > > > -- 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/6ea6f3fa-9a29-4929-801d-d2102c7e5b4bn%40googlegroups.com.
