[petsc-users] Bug report VecNorm

2023-12-09 Thread Stephan Köhler
. 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

Re: [petsc-users] Bug Report TaoALMM class

2023-12-09 Thread Stephan Köhler
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

[petsc-users] Bug Report TaoALMM class

2023-06-12 Thread Stephan Köhler
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

[petsc-users] Bug report LMVM matrix class

2023-05-09 Thread 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

Re: [petsc-users] Report Bug TaoALMM class

2022-11-04 Thread Stephan Köhler
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

Re: [petsc-users] Report Bug TaoALMM class

2022-11-03 Thread Stephan Köhler
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

[petsc-users] Bug report LMVM matrix class

2022-11-02 Thread Stephan Köhler
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

[petsc-users] Report Bug TaoALMM class

2022-10-28 Thread Stephan Köhler
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ß