On 3/14/19 11:42 AM, jane....@jandj-ltd.com wrote: > > as an addition to this post, For Dirichlet condition that are in the > weak form, would this have to be done? > > And for example, for an inhomogeneous normal component of the normal > stress condition in step-22 for Stokes (1st of the partial bc discussion > in the tutorial), do you also have to initialise local_rhs? I have a > similar condition for the top boundary with the Stokes equation but > can't quite figure out if the pressure should equal the normal component > of the normal stress or not (if the stress tensor/velocity is non-zero, > then it shouldn't equal the inhomogeneous condition input?)
I'm not entirely sure I know what you mean with "it" or "this" -- I suspect you are asking whether local_rhs should be set to zero every time we move to a new cell. The answer is yes. We think of these as each cell's contribution to the global matrix. We want to compute it separately on each cell, rather than carry over some state from the previous cell. Best W. -- ------------------------------------------------------------------------ Wolfgang Bangerth email: bange...@colostate.edu www: http://www.math.colostate.edu/~bangerth/ -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout.