Dear Daniel, As you can see my assembly in the first post, my system_matrix and system_rhs are not made directly from cell_system_matrix or cell_system_rhs, but they are the result of some manipulations on mass_matrix, laplase_matrix and nl_matrix (the same approach as step_25) that have already been assembled (in the red part of the code).
To simplify the question, lets say we want to apply any sort of constraints in the assembly part of step-25. I was wondering how to do use ConstraintMatrix::distribute_local_to_global, while there is no cell_system_matrix and cell_system_rhs. Thanks, Hamed -- 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. For more options, visit https://groups.google.com/d/optout.