> say the Neumann bc are imposed at face with boundary indicator 0 > Neumann_boundary['0']=pt; > > And then pass this FunctionMap to the estimate function.
This is the way to go. There's a slight subtlety here: you assign pt to indicator '0', i.e. the character that represents the zero in the screen, but this is not the numerical value zero (which you could write as (char)0). That's not wrong, but you need to make sure you use the same code in all places, i.e. passing it to the estimate function, when assembling the matrix and rhs, etc. > As a result I > did not get any run time exceptions and the code run normally although > as far as I can understand there seems to be some problems with the > convergence of the solver. That must be an unrelated problem. W. ------------------------------------------------------------------------- Wolfgang Bangerth email: [email protected] www: http://www.math.tamu.edu/~bangerth/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
