[firedrake] Problem with Jacobian

[Lawrence Mitchell]

> On 6 Nov 2015, at 12:37, Buesing, Henrik <HBuesing at eonerc.rwth-aachen.de> wrote:
> 
>>> Firedrake uses the UFL derivative function, but does not expose this
>>> extra argument in the interface.  It is straightforward to alter the
>>> definition in firedrake/ufl_expr.py to take this extra argument and
>>> pass it through.
> 
> [Buesing, Henrik]
> 
> @Lawrence: So how would I pass the "coefficient_derivatives" field through?


If you change the definition of derivative in firedrake/ufl_expr.py to take an optional coefficient_derivatives argument:

def derivative(F, u, du=None, coefficient_derivatives=None):
   ...


And just pass that value through to the ufl.derivative call:

   ufl.derivative(..., coefficient_derivatives=coefficient_derivatives)

If you do this in the firedrake source directly, you will be able to commit your change in the local git repository.  You can send us these changes by forking firedrake on github and proposing a pull-request that way.

Cheers,

Lawrence
-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 455 bytes
Desc: Message signed with OpenPGP using GPGMail
URL: <http://mailman.ic.ac.uk/pipermail/firedrake/attachments/20151106/8f772d6c/attachment.sig>