Hi All, Lawrence and I had a productive lunch in which we think we worked out how to do R efficiently. Conversation on this subject will now move to the issue (https://github.com/firedrakeproject/firedrake/issues/226) and all of your input is actively sought. David On 12 March 2014 09:21, Farrell, Patrick <patrick.farrell06@imperial.ac.uk>wrote:

On 12/03/14 09:17, Patrick Farrell wrote:

...

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]),

P.S. For working code where you need this, see

https://bitbucket.org/dolfin-adjoint/da-applications/src/3907be10e58ab0f61e5...

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

-- Dr David Ham Departments of Mathematics and Computing Imperial College London http://www.imperial.ac.uk/people/david.ham