Bernhardt,
> How i obtain the correct boundary values for nodes on the boundary? > or > How can i find all nodes which belong to the actual boundary and set the > correct boundary conditions. > > The problem is that interpolate_boundary_values does not find all nodes on > the > boundary after adaptive refinement. It's not impossible that you're right, but it would also surprise me that 20 years after this function was written, it would still have bugs :-) Do you happen to know whether the missing DoFs happen to be those that are also hanging nodes? > >So to summarize, you are interpolating boundary conditions on boundaries > with > >boundary_id 0, 2, 80, and then you are outputting all DoFs on boundary 80. > > Yes, this problem just happens on boundary 80. > After i gave them as output, i recognized that not all DoFs which schould be > on the boundary are actually at the boundary. I think you'll have to try to narrow down what exactly happens. If you call VectorTools::interpolate (..., 80, ...); then my expectation is that the output argument contains constraints for every DoF located on boundary 80. You can find out which DoFs are located on boundary 80 using DoFTools::extract_boundary_dofs(), and then you could loop over all of these extracted DoFs and ask the AffineConstraints (ConstraintMatrix) class whether these DoFs are constrained using constraints.is_constrained(i); If you get a `false` answer for any of these DoFs, then there is indeed something wrong. You could then double-check whether these DoFs are maybe hanging nodes using DoFTools::extract_hanging_nodes(). Best W. -- ------------------------------------------------------------------------ Wolfgang Bangerth email: [email protected] 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 [email protected]. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/9335d65f-f451-8d44-04d2-c0d2adb05b99%40colostate.edu. For more options, visit https://groups.google.com/d/optout.
