← Back to team overview

ufl team mailing list archive

Re: Dirichlet boundary conditions

 

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.

--
Martin


Follow ups

References