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

_______________________________________________
DOLFIN-dev mailing list
[email protected]
http://www.fenics.org/mailman/listinfo/dolfin-dev

Reply via email to