Hello everyone. I am working on a problem very similar to the Laplace equation solved in step-3. I have made some changes to step-3 as I need to. Now, I want to solve this problem in parallel. I looked at step-40, which solves the Laplace equation in parallel. However, it uses a constraints object to apply BCs. Right now, I don't have any hanging nodes to handle, so I used a std::map object to apply BCs. I first copy all the local contributions to the global matrix and then apply the BCs. While going through the documentation of step-40, I read that *"Copying local contributions into the global matrix must include distributing constraints and boundary values. In other words, we cannot (as we did in step-6 <https://www.dealii.org/current/doxygen/deal.II/step_6.html>) first copy every local contribution into the global matrix and only in a later step take care of hanging node constraints and boundary values."*
If my understanding is correct, then I think the following line of code does this job in step-40: *constraints.distribute_local_to_global(cell_matrix,cell_rhs,local_dof_indices,* * system_matrix,system_rhs);* My doubt is whether it is possible to do this(*Copying local contributions into the global matrix and distributing boundary values.*) with the std::map object or do I need to use a constraint object even though I don't have any hanging nodes to handle. Thanks and regards Wasim -- 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/76972977-2325-4cc3-9822-560e520ba99bn%40googlegroups.com.
