Hi, I want to assemble my system matrix only once since it doesn't change throughout the time steps. For that my code is similar to step-57:
if (assemble_matrix){ constraints.distribute_local_to_global(local_matrix, local_rhs, local_dof_indices, system_matrix, system_rhs); } else{ constraints.distribute_local_to_global(local_rhs, local_dof_indices, system_rhs); } but I don't seem to get the same *system_rhs *when *assemble_system=false* as when *assemble_system=true*. To be more specific, it seems that there is something wrong with how constraints.distribute_local_to_global(local_rhs, local_dof_indices, system_rhs) creates the system_rhs from local_rhs and local_dofs_indices. I don't change anything about constraints etc throughout the time steps since I set those things up at the very beginning and then I just run several times without updating anything. Any idea about what the problem could be? Best, Magda -- 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/174a9f9e-491a-47f2-9e36-d737bca51d13%40googlegroups.com.