← Back to team overview

dolfin team mailing list archive

Re: Functions for applying boundary conditions to the RHS vector only now accept the solution vector x. Where Dirichlet boundary conditions are applied, b = x - bc.The is useful for applying boundary conditions within Newton iterations. [...]

 

> On Mon, Dec 05, 2005 at 08:29:27AM +0100, Johan Hoffman wrote:
>
>> Isn't this the natural way to deal with Neumann bc? Just add a surface
>> integral to the variational form?
>
> Yes. Maybe there can be some automation for specifying boundary
> conditions on a higher level, but the correct way to do this is
> probably to let each module specify a list of possible boundary
> conditions, since the available boundary conditions will depend on the
> equation (like "outflow boundary conditions" for Navier-Stokes).
> The module would choose how to implement the boundary conditions, like
> strong or weak.
>
> /Anders
>

Yes, this should probably be handled on module-level.

/Johan




References