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, well it actually sets the rhs to zero. And this is the current behavior I observe in step 11 of the tutorial. Even though the documentation says it adds 1 for the neumman boundary it is acutally adding zero. I wonder why the solution still "converges"???
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. -Sebastian _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
