Dear colleagues, I am trying to solve my problem with different right_hand_sides.
For each right hand side, we can solve by the codes: LinearAlgebraPETSc::MPI::Vector completely_distributed_solution(locally_owned_dofs, mpi_communicator); SolverControl solver_control; PETScWrappers::SparseDirectMUMPS solver(solver_control, mpi_communicator); solver.solve(system_matrix, completely_distributed_solution, system_rhs); constraints.distribute(completely_distributed_solution); locally_relevant_solution = completely_distributed_solution; I think I can write the codes to assemble the different right hand sides to solve the system with the same system matrix. This give me the different solutions. However, I do not know if I can sum up all the completely_distributed_solution from the different right-hand-sides in order to get a vector of solution. Could you please tell me how to do this? I would like to thank you very much for your support. Best regards, Kien -- 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/d2272747-cbc2-475a-85bc-7fb02f9cce57%40googlegroups.com.