Hi George, I think I don't quite understand what you mean. Can you write out mathematically what you're after and I'll try to understand if it's something we can extract from Dolfin-Adjoint? In addition, I need to apologise in advance because I'm off sailing all of next week so I shan't be very responsive. I will be back online in a week's time. Regards, David On Mon, 10 Apr 2017 at 18:29 George Ovchinnikov <lives9@gmail.com> wrote:
Dear colleagues, sorry for hijacking the thread, but could you explain to me, how to compute residue of adjoint functional, if possible, with firedrake-adjoint?
I've seen in the documentation that
J = Functional(inner(u, u)*dx*dt[FINISH_TIME])for (variable, solution) in compute_adjoint(J):
provides the solution of adjoint problem, but for residual I need not only solution, but also a set of adjoint equations.
Sincerely, George
On 10.04.2017 17:32, Buesing, Henrik wrote:
Yes! I opened a github issue here <https://github.com/firedrakeproject/firedrake/issues/1042>. I guess the markdown messes with the python code. Sorry, for that.
Henrik
--
Dipl.-Math. Henrik Büsing
Institute for Applied Geophysics and Geothermal Energy
E.ON Energy Research Center
RWTH Aachen University
Mathieustr. | Tel +49 (0)241 80 49907 <+49%20241%208049907>
52074 Aachen, Germany | Fax +49 (0)241 80 49889 <+49%20241%208049889>
http://www.eonerc.rwth-aachen.de/GGE
hbuesing@eonerc.rwth-aachen.de
*Von:* firedrake-bounces@imperial.ac.uk [ mailto:firedrake-bounces@imperial.ac.uk <firedrake-bounces@imperial.ac.uk>] *Im Auftrag von *David Ham *Gesendet:* Montag, 10. April 2017 11:57 *An:* firedrake <firedrake@imperial.ac.uk> <firedrake@imperial.ac.uk> *Betreff:* Re: [firedrake] firedrake_adjoint
Hi Henrik,
That sounds like a bug. Can you create a minimal failing example and file a github issue please?
Regards,
David
On Mon, 10 Apr 2017 at 10:48 Buesing, Henrik < HBuesing@eonerc.rwth-aachen.de> wrote:
Dear all,
If I include
*from firedrake*_*adjoint import **
in my firedrake program, my forward problem does not converge any more. If I do nothing else than loading firedrake_adjoint, should the forward problem not run as before? Note, that I set all the solver options from the command line and not in the program itself. Thus, in the program I just have
*solver = NonlinearVariationalSolver(problem,options_prefix="")*
Thank you! Henrik
--
Dipl.-Math. Henrik Büsing
Institute for Applied Geophysics and Geothermal Energy
E.ON Energy Research Center
RWTH Aachen University
------------------------------------------------------
Mathieustr. 10 | Tel +49 (0)241 80 49907 <+49%20241%208049907>
52074 Aachen, Germany | Fax +49 (0)241 80 49889 <+49%20241%208049889>
------------------------------------------------------
http://www.eonerc.rwth-aachen.de/GGE
hbuesing@eonerc.rwth-aachen.de
------------------------------------------------------
--
Dr David Ham
Department of Mathematics
Imperial College London
_______________________________________________ firedrake mailing listfiredrake@imperial.ac.ukhttps://mailman.ic.ac.uk/mailman/listinfo/firedrake
-- Dr David Ham Department of Mathematics Imperial College London