Re: [petsc-users] Report Bug TaoALMM class

2022-11-11 Thread Barry Smith
I noticed this in the TAOALMM manual page. It is also highly recommended that the subsolver chosen by the user utilize a trust-region strategy for globalization (default: TAOBQNKTR) especially if the outer problem features bound constraints. I am far from an expert on these topics. > On

Re: [petsc-users] Report Bug TaoALMM class

2022-11-11 Thread Barry Smith
I am still working to understand. I have a PETSc branch barry/2022-11-11/fixes-for-tao/release where I have made a few fix/improvements to help me run and debug with your code. I made a tiny change to your code, passing Hessian twice, and ran with ./test_tao_neohooke

Re: [petsc-users] Use of LSC preconditioning

2022-11-11 Thread Alexander Lindsay
That makes sense. Thanks for the quick reply! On Fri, Nov 11, 2022 at 2:59 PM Matthew Knepley wrote: > On Fri, Nov 11, 2022 at 5:57 PM Alexander Lindsay < > alexlindsay...@gmail.com> wrote: > >> Under what conditions can I use LSC preconditioning for field split >> problems with Schur? Let's

Re: [petsc-users] Use of LSC preconditioning

2022-11-11 Thread Matthew Knepley
On Fri, Nov 11, 2022 at 5:57 PM Alexander Lindsay wrote: > Under what conditions can I use LSC preconditioning for field split > problems with Schur? Let's imagine that all I've done is called > SNESetJacobian with global A and P and provided the index sets for 0 and 1. > Based off of the

[petsc-users] Use of LSC preconditioning

2022-11-11 Thread Alexander Lindsay
Under what conditions can I use LSC preconditioning for field split problems with Schur? Let's imagine that all I've done is called SNESetJacobian with global A and P and provided the index sets for 0 and 1. Based off of the documentation on the man page

Re: [petsc-users] On PCFIELDSPLIT and its implementation

2022-11-11 Thread Edoardo alinovi
Hi guys, anyone with any suggestion to make this thing work? 梁 Il Ven 11 Nov 2022, 10:24 Edoardo alinovi ha scritto: > Hi Barry, > > FYI, in test.F90 I noted that "ui" starts for 1 and not from 0. I fixed it > but the situation does not change much. I attached the new file in this > email. >

Re: [petsc-users] Report Bug TaoALMM class

2022-11-11 Thread Barry Smith
> On Nov 4, 2022, at 7:43 AM, Stephan Köhler > wrote: > > Barry, > > this is a nonartificial code. This is a problem in the ALMM subsolver. I > want to solve a problem with a TaoALMM solver what then happens is: > > TaoSolve(tao)/* TaoALMM solver */ >| >| >