[firedrake] Problem with Jacobian

[Buesing, Henrik]

> > Are you sure the problem isn't that he needs bounds for his solve?
> > AFAICS dS/dh=0 for h<hw or h>hn. Otherwise shouldn't rewriting this
> > with conditionals have worked?

@Stephan: So how can I rewrite the attached saturation.py with conditionals?

Henrik

-------------- next part --------------
A non-text attachment was scrubbed...
Name: saturation.py
Type: application/octet-stream
Size: 504 bytes
Desc: saturation.py
URL: <http://mailman.ic.ac.uk/pipermail/firedrake/attachments/20151106/c235fb0a/attachment.obj>