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 [email protected].
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/174a9f9e-491a-47f2-9e36-d737bca51d13%40googlegroups.com.
