Yes, that part now works. Thank you, Wolfgang. Anton.
On Wednesday, March 17, 2021 at 3:13:17 PM UTC-7 Wolfgang Bangerth wrote: > On 3/17/21 10:41 PM, Anton Ermakov wrote: > > Thanks for the reply. I think I got it. Basically, I had to manually > place the > > local matrices into the global matrix, without taking into account the > hanging > > node constraints, which are later taken care of by > > > > /constraints.condense(system_matrix, system_rhs);/ > > > > In the case of the mass matrix, I had to use this: > > > > /for (unsigned int i = 0; i < dofs_per_cell; ++i)/ > > > > /for (unsigned int j = 0; j < dofs_per_cell; ++j)/ > > > > > /mass_matrix.add(local_dof_indices[i], local_dof_indices[j], cell_matrix(i, > j));/ > > > > Instead of > > > > /constraints.distribute_local_to_global(cell_matrix, local_dof_indices, > > mass_matrix);/ > > Ah yes, that makes sense. So do I understand right that whatever you have > now > actually works? > > Best > W. > > > -- > ------------------------------------------------------------------------ > Wolfgang Bangerth email: bang...@colostate.edu > www: http://www.math.colostate.edu/~bangerth/ > > -- 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/fdc7b151-885e-48d9-aace-8e25e8669597n%40googlegroups.com.
