<!DOCTYPE html>
<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
</head>
<body>
<pre class="moz-quote-pre" wrap="">Dear PETSc/Tao team,
this is still an open issue and <span _d-id="4988"
class="--l --r sentence_highlight">I haven't heard anything else so far that I'm wrong.
</span><span _d-id="4988" class="--l --r sentence_highlight"></span></pre>
<span style="white-space: pre-wrap">Kind regards,</span><br>
Stephan Köhler<br>
<br>
<div class="moz-cite-prefix">Am 18.07.23 um 02:21 schrieb Matthew
Knepley:<br>
</div>
<blockquote type="cite"
cite="mid:CAMYG4GmafqenDV=i9tJZ=KEY=nnbdinaRCZBTztRMvQ-2acr-A@mail.gmail.com">
<pre class="moz-quote-pre" wrap="">Toby and Hansol,
Has anyone looked at this?
Thanks,
Matt
On Mon, Jun 12, 2023 at 8:24 AM Stephan Köhler <
<a class="moz-txt-link-abbreviated" href="mailto:stephan.koehler@math.tu-freiberg.de">stephan.koehler@math.tu-freiberg.de</a>> wrote:
</pre>
<blockquote type="cite">
<pre class="moz-quote-pre" wrap="">Dear PETSc/Tao team,
I think there might be a bug in the Tao ALMM class: In the function
TaoALMMComputeAugLagAndGradient_Private(), see, eg.
<a class="moz-txt-link-freetext" href="https://petsc.org/release/src/tao/constrained/impls/almm/almm.c.html#TAOALMM">https://petsc.org/release/src/tao/constrained/impls/almm/almm.c.html#TAOALMM</a>
line 648 the gradient seems to be wrong.
The given function and gradient computation is
Lc = F + Ye^TCe + Yi^T(Ci - S) + 0.5*mu*[Ce^TCe + (Ci - S)^T(Ci - S)],
dLc/dX = dF/dX + Ye^TAe + Yi^TAi + 0.5*mu*[Ce^TAe + (Ci - S)^TAi],
but I think the gradient should be (without 0.5)
dLc/dX = dF/dX + Ye^TAe + Yi^TAi + mu*[Ce^TAe + (Ci - S)^TAi].
Kind regards,
Stephan Köhler
--
Stephan Köhler
TU Bergakademie Freiberg
Institut für numerische Mathematik und Optimierung
Akademiestraße 6
09599 Freiberg
Gebäudeteil Mittelbau, Zimmer 2.07
Telefon: +49 (0)3731 39-3173 (Büro)
</pre>
</blockquote>
<pre class="moz-quote-pre" wrap="">
</pre>
</blockquote>
<br>
<pre class="moz-signature" cols="72">--
Stephan Köhler
TU Bergakademie Freiberg
Institut für numerische Mathematik und Optimierung
Akademiestraße 6
09599 Freiberg
Gebäudeteil Mittelbau, Zimmer 2.07
Telefon: +49 (0)3731 39-3188 (Büro)</pre>
</body>
</html>