> As I refine the mesh (1D), I'd need to scale the right-hand-side > appropriately as in that case, the solution would grow as 1/h.
If that is the case then you must have forgotten to multiply it with fe_values.JxW(q) Check your code and compare it with how step-3, for example, does that. > Also, how would I go about applying a non-homogeneous Neumann boundary > condition quickly at the right end? I know that Step-7 shows an example > but I was just wondering if there was a simpler or another way. That's the recommended way to do it. Best W. ------------------------------------------------------------------------- Wolfgang Bangerth email: [email protected] www: http://www.math.tamu.edu/~bangerth/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
