.
Otherwise the program may hang when PetscCheck is executed.
Please find a minimal example attached.
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
Dear PETSc/Tao team,
this is still an open issue andI haven't heard anything else so far that I'm
wrong.
Kind regards,
Stephan Köhler
Am 18.07.23 um 02:21 schrieb Matthew Knepley:
Toby and Hansol,
Has anyone looked at this?
Thanks,
Matt
On Mon, Jun 12, 2023 at 8:24 AM Stephan
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
overwrites also the value in Y.
Best 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)
OpenPGP_0xC9BF2C20DFE9F713.asc
things is to come up with an appropriate correction for the current code.
Barry
On Nov 3, 2022, at 12:46 PM, Stephan Köhler
wrote:
Barry,
so far, I have not experimented with trust-region methods, but I can
imagine that this "design feature" causes no problem for trust-r
ure". Within the
standard use of ALMM I do not see how the objective function would be used within a line
search. The TaoSolve_ALMM() code is self-correcting in that if a trust region check fails
it automatically rolls back the solution.
Barry
On Oct 28, 2022, at 4:27
to the PETSc
installation in the makefile.
Best 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
hat happens is that the point ((xk + dxk) + alpha * dxk) is
evaluated, see, e.g.,
https://petsc.org/release/src/tao/linesearch/impls/armijo/armijo.c.html
line 191.
Best regards
Stephan Köhler
--
Stephan Köhler
TU Bergakademie Freiberg
Institut für numerische Mathematik und Optimierung
Akademiestraß