Hi Sebastian, > Unless something is wrong in my understanding I believe that if > VectorTools::create_boundary_right_hand_side is called with the > default argument for the boundary_indicators then it does nothing,
that is correct and probably not what is expected from the function. At least step-11 does not what it says it does right now. > I wonder why the solution still > "converges"??? Good question, maybe someone else can jump in here? The errors are a little smaller after fixing this error. > The issue is quite simple to fix, but it requires a design decision: > either remove the default argument or modify the code so that if an > empty set is passed then all boundary faces are assumed. I checked in that an empty set assumes all boundary faces. Can you test that it works as expected now, please? -- Timo Heister http://www.math.tamu.edu/~heister/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
