← Back to team overview

ufl team mailing list archive

Re: Dirichlet boundary conditions

 

On Mon, Jun 23, 2008 at 11:54:01AM +0200, Martin Sandve Alnæs wrote:
> 2008/6/23 Martin Sandve Alnæs <martinal@xxxxxxxxx>:
> > 2008/6/23  <kent-and@xxxxxxxxx>:
> >>> If we should add notation for boundary conditions to UFL, then I think
> >>> it must also be supported through the UFC interface so that the
> >>> assembler may assemble the correct matrix/vector directly (without
> >>> zeroing out entries in a separate step).
> >>
> >> I don't understand what you mean with this. It has to be a separate step
> >> that is performed after the tabulate_tensor functions.
> >
> > The assembler cannot use UFL directly, and UFL is supposed to pass
> > through a form compiler. Somebody must convert the UFL to something
> > the C++ assembler can eat, and UFC is the interface we hide equation
> > details behind. Unless you have a better option?
> 
> However, I agree that it should probably be a separate step in the
> assembler after calling tabulate_tensor but before calling A.add.

Yes, and the important thing is that it's communicated to the
assembler using the UFC interface (which we may need to extend).

It should not be something that is communicated from UFL --> DOLFIN
without UFC.

-- 
Anders

Attachment: signature.asc
Description: Digital signature


References