Hamed, Let me ask a more specific question. I was wondering how to apply any sort > of constraints ( that may be periodic BC or so) in an assembly function > like mine in which some manipulations is done on mass_matrix and > laplase_matrix. > Following some posts in forum I found I can use void > ConstraintMatrix::condense (system_matrix, system_rhs) at the end of > assembly after all manipulations, but it seems this function works with > dealii sparsematrices not Trilinos ones. > The recommended way is to use ConstraintMatrix::distribute_local_to_global and ConstraintMatrix::distribute. Furthermore, ConstraintMatrix::condense for matrices is not possible in parallel and only implemented for the serial SparseMatrix and BlockSparseMatrix. For discussing these issues, have a look at "Constraints on degrees of freedom" [1].
Best, Daniel [1] https://www.dealii.org/8.4.1/doxygen/deal.II/group__constraints.html -- 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.
