Dear Denis
I checked the the system residual assembly process. For a triangulation with 28 cells, the dof index 17 is belong to cells 2, 5, 8 and 27. The corresponding cell_rhs values for this index are, cell cell_rhs value 2 2.88462e+00 5 2.88462e+00 8 -2.88462e+00 27 -1.44231e+00 This dof is not constrained. When the following assembly method is used: *for* (*unsigned* *int* i = 0; i < dofs_per_cell; ++i) system_rhs(data.local_dof_indices[i]) += data.cell_rhs(i); the system_rhs(17) = 1.4423076, however, using the following method, constraints_update.distribute_local_to_global(data.cell_rhs, data.local_dof_indices, system_rhs); the system_rhs(17) = -1.33e-15. This is the source of the problem because I have used the second method in the previous code and the first one (similar to step-44) in the new code. I cant understand why the results are different. Regards P. On Wednesday, February 1, 2017 at 8:24:51 PM UTC+3:30, Daniel Arndt wrote: > > Pasha, > > 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 >> > You are providing way too few information to expect useful answers. What > are you doing in your code? > For which part in your code do you want to change to WorkStream? > What exactly did you change? How is this related to whether you are doing > adaptive refinement or not? > > Best, > Daniel > -- 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.
