Re: [petsc-users] issues with VecSetValues in petsc 3.19

2023-05-10 Thread Jed Brown
Edoardo alinovi writes: > Hello Barry, > > Welcome to the party! Thank you guys for your precious suggestions, they > are really helpful! > > It's been a while since I am messing around and I have tested many > combinations. Schur + selfp is the best preconditioner, it converges within > 5 iters

Re: [petsc-users] issues with VecSetValues in petsc 3.19

2023-05-05 Thread Edoardo alinovi
Hello Barry, Welcome to the party! Thank you guys for your precious suggestions, they are really helpful! It's been a while since I am messing around and I have tested many combinations. Schur + selfp is the best preconditioner, it converges within 5 iters using gmres for inner solvers but it is

Re: [petsc-users] issues with VecSetValues in petsc 3.19

2023-05-05 Thread Barry Smith
To expand on what Matt said slightly. When you have a preconditioner based on (possibly nested) sub solves one generally "tunes" the solves to minimize time to solution. We recommend doing this by first using very accurate subsolves (when possible using direct solves inside); this tells us

Re: [petsc-users] issues with VecSetValues in petsc 3.19

2023-05-05 Thread Matthew Knepley
On Fri, May 5, 2023 at 5:13 AM Edoardo alinovi wrote: > Hi Matt, > > I have some more questions on the fieldsplit saga :) > > I am running a 1M cell ahmed body case using the following options: > > "solver": "fgmres", > "preconditioner": "fieldsplit", > "absTol": 1e-6, > "relTo

Re: [petsc-users] issues with VecSetValues in petsc 3.19

2023-05-05 Thread Edoardo alinovi
Hi Matt, I have some more questions on the fieldsplit saga :) I am running a 1M cell ahmed body case using the following options: "solver": "fgmres", "preconditioner": "fieldsplit", "absTol": 1e-6, "relTol": 0.0, "options":{ "pc_fieldsplit_type": "multiplicative",

Re: [petsc-users] issues with VecSetValues in petsc 3.19

2023-04-21 Thread Matthew Knepley
On Thu, Apr 20, 2023 at 11:07 AM Edoardo alinovi wrote: > Hi Matt, > > thanks for sharing the literature. > > Would you suggest any monolitic approach for the mpiaij/mpibaij matrix > instead of fieldsplit? > I have not seen them be better, but Vanka-type smoothers can work for this system. > I

Re: [petsc-users] issues with VecSetValues in petsc 3.19

2023-04-20 Thread Edoardo alinovi
Hi Matt, thanks for sharing the literature. Would you suggest any monolitic approach for the mpiaij/mpibaij matrix instead of fieldsplit? I did some blind search using gamg/hypre and they look terribile. I guess i am missing a trick, probaly they are not the way to go? Thanks! Il Lun 17 Apr 20

Re: [petsc-users] issues with VecSetValues in petsc 3.19

2023-04-17 Thread Matthew Knepley
On Mon, Apr 17, 2023 at 6:37 AM Edoardo alinovi wrote: > Sure thing, the solver I am working on is this one: > https://gitlab.com/alie89/flubio-code-fvm. > > It is a 3D, collocated, unstructured, finite volume solver for > incompressibility NS. I can run steady, unsteady and I can use SIMPLE, PIS

Re: [petsc-users] issues with VecSetValues in petsc 3.19

2023-04-17 Thread Edoardo alinovi
Sure thing, the solver I am working on is this one: https://gitlab.com/alie89/flubio-code-fvm. It is a 3D, collocated, unstructured, finite volume solver for incompressibility NS. I can run steady, unsteady and I can use SIMPLE, PISO and Factional step method (both explicit and fully implicit mome

Re: [petsc-users] issues with VecSetValues in petsc 3.19

2023-04-17 Thread Matthew Knepley
On Mon, Apr 17, 2023 at 6:16 AM Edoardo alinovi wrote: > Do you mean the solver I am messing around? XD > Yes, and what physics it is targeting. THanks, Matt -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results t

Re: [petsc-users] issues with VecSetValues in petsc 3.19

2023-04-17 Thread Edoardo alinovi
Do you mean the solver I am messing around? XD

Re: [petsc-users] issues with VecSetValues in petsc 3.19

2023-04-17 Thread Matthew Knepley
On Mon, Apr 17, 2023 at 6:09 AM Edoardo alinovi wrote: > Thanks Matt, your always there when you need <3 > Glad it's working! Sometime you have to tell me what it is solving. Thanks, Matt -- What most experimenters take for granted before they begin their experiments is infinitely mor

Re: [petsc-users] issues with VecSetValues in petsc 3.19

2023-04-17 Thread Edoardo alinovi
Thanks Matt, your always there when you need <3

Re: [petsc-users] issues with VecSetValues in petsc 3.19

2023-04-17 Thread Edoardo alinovi
Aaah yes you are right. Do not ask me why, but I was not getting this with 3.18.5, odd.

Re: [petsc-users] issues with VecSetValues in petsc 3.19

2023-04-17 Thread Matthew Knepley
On Mon, Apr 17, 2023 at 6:00 AM Edoardo alinovi wrote: > Hey Matt, > > Thanks for the help. Here is the error: > > [0]PETSC ERROR: - Error Message > -- > [0]PETSC ERROR: Object is in wrong state > [0]PETSC ERROR: Not

Re: [petsc-users] issues with VecSetValues in petsc 3.19

2023-04-17 Thread Edoardo alinovi
Hey Matt, Thanks for the help. Here is the error: [0]PETSC ERROR: - Error Message -- [0]PETSC ERROR: Object is in wrong state [0]PETSC ERROR: Not for unassembled vector, did you call VecAssemblyBegin()/VecAssemblyEnd(

Re: [petsc-users] issues with VecSetValues in petsc 3.19

2023-04-17 Thread Matthew Knepley
On Mon, Apr 17, 2023 at 5:36 AM Edoardo alinovi wrote: > Hello Barry, Matt, Jed, > > I have just installed the latest and greatest version of petsc and I am > hitting a problem I did not have in previous releases. > > Here is the error: > > > > > *[1]PETSC ERROR: - Error Messa

[petsc-users] issues with VecSetValues in petsc 3.19

2023-04-17 Thread Edoardo alinovi
Hello Barry, Matt, Jed, I have just installed the latest and greatest version of petsc and I am hitting a problem I did not have in previous releases. Here is the error: *[1]PETSC ERROR: - Error Message --[1]PETSC