Re: [deal.II] Problem with "constraints.distribute_local_to_global(local_rhs, local_dof_indices, system_rhs)"

2020-03-05 Thread Magdalini Ntetsika
I see. Thank you Daniel! Best, Magda On Thursday, March 5, 2020 at 1:39:26 PM UTC-8, Daniel Arndt wrote: > > Magdalini, > > The problem only appears with inhomogeneous boundary conditions, but > periodic boundary conditions don't set any inhomogeneity. > > Best, > Daniel > > Am Do., 5. März 2020

Re: [deal.II] Problem with "constraints.distribute_local_to_global(local_rhs, local_dof_indices, system_rhs)"

2020-03-05 Thread Daniel Arndt
Magdalini, The problem only appears with inhomogeneous boundary conditions, but periodic boundary conditions don't set any inhomogeneity. Best, Daniel Am Do., 5. März 2020 um 16:13 Uhr schrieb Magdalini Ntetsika < ntets...@berkeley.edu>: > Hi Wolfgang, > > thank you for the prompt reply. I see,

Re: [deal.II] Problem with "constraints.distribute_local_to_global(local_rhs, local_dof_indices, system_rhs)"

2020-03-05 Thread Magdalini Ntetsika
Hi Wolfgang, thank you for the prompt reply. I see, I'll look up both. But I have a quick question: what if I have periodic boundary conditions all around my domain? It seems to be working in that case, even if the boundary values are nonzero, I get to have the correct system_rhs vector. I gues

Re: [deal.II] Problem with "constraints.distribute_local_to_global(local_rhs, local_dof_indices, system_rhs)"

2020-03-05 Thread Wolfgang Bangerth
On 3/5/20 9:24 AM, Magdalini Ntetsika wrote: 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, sy

[deal.II] Problem with "constraints.distribute_local_to_global(local_rhs, local_dof_indices, system_rhs)"

2020-03-05 Thread Magdalini Ntetsika
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{ cons