Re: [petsc-users] BJACOBI with FIELDSPLIT

2019-03-18 Thread Matthew Knepley via petsc-users
Justin Chang > *Cc:* Smith, Barry F.; petsc-users@mcs.anl.gov > *Subject:* Re: [petsc-users] BJACOBI with FIELDSPLIT > > > Thanks, using fgmres it does work as expected. > > I thought gmres would do the same since I'm solving the subblocks > "exactly". > > &

Re: [petsc-users] BJACOBI with FIELDSPLIT

2019-03-18 Thread Matthew Knepley via petsc-users
On Mon, Mar 18, 2019 at 3:18 PM Rossi, Simone via petsc-users < petsc-users@mcs.anl.gov> wrote: > Dear all, > > I'm debugging my application in which I'm trying to use the FIELDSPLIT > preconditioner for solving a 2x2 block matrix. > > > Currently I'm testing the preconditioner on a decoupled

Re: [petsc-users] BJACOBI with FIELDSPLIT

2019-03-18 Thread Rossi, Simone via petsc-users
: [petsc-users] BJACOBI with FIELDSPLIT Thanks, using fgmres it does work as expected. I thought gmres would do the same since I'm solving the subblocks "exactly". Simone From: Justin Chang Sent: Monday, March 18, 2019 3:38:34 PM To: Rossi, Simone

Re: [petsc-users] BJACOBI with FIELDSPLIT

2019-03-18 Thread Rossi, Simone via petsc-users
cs.anl.gov Subject: Re: [petsc-users] BJACOBI with FIELDSPLIT Use -ksp_type fgmres if your inner ksp solvers are gmres. Maybe that will help? On Mon, Mar 18, 2019 at 1:33 PM Rossi, Simone via petsc-users mailto:petsc-users@mcs.anl.gov>> wrote: Thanks Barry. Let me know if you can s

Re: [petsc-users] BJACOBI with FIELDSPLIT

2019-03-18 Thread Justin Chang via petsc-users
; linear system matrix = precond matrix: > > Mat Object: (fieldsplit_1_) 1 MPI processes > > type: seqaij > > rows=35937, cols=35937 > > total: nonzeros=912673, allocated nonzeros=912673 > > total number of mallocs us

Re: [petsc-users] BJACOBI with FIELDSPLIT

2019-03-18 Thread Rossi, Simone via petsc-users
, Simone Cc: petsc-users@mcs.anl.gov Subject: Re: [petsc-users] BJACOBI with FIELDSPLIT Simone, This is indeed surprising, given the block structure of the matrix and the exact block solves we'd expect the solver to converge after the application of the preconditioner. Please send the output

Re: [petsc-users] BJACOBI with FIELDSPLIT

2019-03-18 Thread Smith, Barry F. via petsc-users
Simone, This is indeed surprising, given the block structure of the matrix and the exact block solves we'd expect the solver to converge after the application of the preconditioner. Please send the output of -ksp_view Barry Also if you are willing to share your test code we can try

[petsc-users] BJACOBI with FIELDSPLIT

2019-03-18 Thread Rossi, Simone via petsc-users
Dear all, I'm debugging my application in which I'm trying to use the FIELDSPLIT preconditioner for solving a 2x2 block matrix. Currently I'm testing the preconditioner on a decoupled system where I solve two identical and independent Poisson problems. Using the default fieldsplit type