dolfin team mailing list archive

[Marie Rognes <meg@xxxxxxxxxxx>]

Kristian Oelgaard wrote:
> You need to compile the form first:
> (compiled_form, module, form_data) = jit(a)
> bc.apply(A, b, compiled_form)
>
> Have a look in src/pydolfin/assemble.py (LinearPDE)
>
>   

Right. Thanks!

> > In addition, if I try just using
> >
> >     zero = Function(mesh, 0.0)
> >     
>
> The function 'f' has to be a FFC function when creating your form. So you have
> to initialise it with an element.
>   

Yes, makes sense now. So what would be the prettiest way of writing the 
zero function for use with a form?


--
Marie E. Rognes
Ph.D Fellow, 
Centre of Mathematics for Applications, 
University of Oslo
http://folk.uio.no/meg