Hamed, [...] > However, to be able to use ConstraintMatrix::distribute_local_to_global, I > can assemble system_matrix in local level instead of global level. Having > that said, I can not do so for the system_rhs because solution and > old_solution vectors are global vectors. > You can obtain local DoF values from a global FE vector using FEValuesBase::get_function_values [1], create the right-hand side locally as you wish and use ConstraintMatrix::distribute_local_to_global to create the global right-hand side. This is outlined in step-21 for example.
Best, Daniel [1] https://www.dealii.org/8.4.1/doxygen/deal.II/classFEValuesBase.html#a3f84ee49f0fbb32bfe01b2f4a473bb47 [2] https://www.dealii.org/8.4.1/doxygen/deal.II/step_21.html#TwoPhaseFlowProblemassemble_system -- 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]. For more options, visit https://groups.google.com/d/optout.
