On 12/03/14 09:17, Patrick Farrell wrote:P.S. For working code where you need this, see
> R is useful for things other than boundary conditions, too. If you have a constrained
> optimisation problem with a single constraint (an integral constraint, say, like [0]),
https://bitbucket.org/dolfin-adjoint/da-applications/src/3907be10e58ab0f61e5d4ea21f1703fb145b15ec/stokes_topology/ex4.py?at=default
which solves the KKT equations associated with example 4 of [0] by means
of a combined continuation on log barrier terms and casting the problem
as a mixed complementarity problem to the PETSc VI solver.
I'd be dead impressed if firedrake could do this (even more so if it could
do this on GPUs).
Patrick
_______________________________________________
firedrake mailing list
firedrake@imperial.ac.uk
https://mailman.ic.ac.uk/mailman/listinfo/firedrake