On Mon, Jun 8, 2009 at 4:13 PM, Garth N. Wells<gnw20@xxxxxxxxx> wrote:
Martin Sandve Alnæs wrote:
I'll take a look at improving the error message and checking,
but try this (if I understand what you're trying to do correctly):
kc = Function(mixed_element)
k1, c1 = split(kc)
dkc = TrialFunction(mixed_element)
dk, dc = split(dkc)
...
a1 = derivative(L1, kc, dkc)
Thanks. I managed to figure it out and it's now used in the DOLFIN
Cahn-Hilliard demo.
Garth
Great. It doesn't seem to be checked in though.