Dear all
I have a code that works fine. Now, I am working to change the code structure to use “workstream”. The new code works fine without adaptive refinement. However, when the mesh is refined and the triangulation has hanging node, a constant residual error remains during Newton iteration. I cant find source of the problem. There are two “constraints matrix”. On of them for “Newton update” and one for “hanging nodes”. The constraints matrix are merged and is used to determine error as: BlockVector<*double*> error_ud(dofs_per_block); *for* (*unsigned* *int* i = 0; i < dof_handler_ref.n_dofs(); ++i) *if* (!constraints_update.*is_constrained*(i)) error_ud(i) = system_rhs(i); error_update.norm = error_ud.l2_norm(); Regards Pasha -- 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.
