Re: [deal.II] Parallization Step-40

2024-10-14 Thread James Shi
Hi, These lines in laplace_cell() (especially the reinit function) may cause MPI to stuck, LA::MPI::Vector evant_solution; evant_solution.reinit(locally_owned_dofs, locally_relevant_dofs, mpi_communicator); I cannot figure out why you need them. Is this a typo? Junxiang Wang 于2024年10月10日周四 22:

Re: [deal.II] Parallization Step-40

2024-10-10 Thread Wolfgang Bangerth
On 10/10/24 08:47, Junxiang Wang wrote: The program didn't proceed after I made a subroutine called              laplace_cell(cell,                                  cell_matrix,                                  cell_rhs); within the assemble_system() function. On a given MPI process, yo

[deal.II] Parallization Step-40

2024-10-10 Thread Junxiang Wang
Dear Prof. Timo Heister and Wolfgang Bangerth I am trying to solve my multiphysics problem by deal.ii and it's paralization. I started from STEP-40 and encountered a problem regarding paralization and MPI. The program didn't proceed after I made a subroutine called laplace_cel