> Are there any suggestions on how this issue should be handled? I'm not > sure it is possible to automate this internally to ConstraintMatrix > (i.e., transparently to the user). > > One quick fix to interpolate_boundary_values, which doesn't require > major changes in the library, would be to check wether the given dof > is already constrained, in which case the dirichlet inhomogeineity > constraint is simply ignored.
I've been thinking about this procedure, too. Boundary values should usually be the weakest part (periodic b.c. or hanging nodes should dominate), and that would work fine for most user codes. The only thing we have to make very clear then is that hanging nodes should be inserted before boundary constraints. And of course document that dofs that already are constrained are not touched in the interpolate_boundary_values function. I can make this change if we all agree on that this is the best way to do things. Best, Martin _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
