Re: [petsc-users] Confusing Schur preconditioner behaviour

2019-03-26 Thread Cotter, Colin J via petsc-users
Hi Dave, Thanks for the tip - you were right, and this works better for higher resolutions now. all the best --Colin From: Dave May Sent: 19 March 2019 11:25:11 To: Cotter, Colin J Cc: petsc-users@mcs.anl.gov Subject: Re: [petsc-users] Confusing Schur

Re: [petsc-users] Confusing Schur preconditioner behaviour

2019-03-19 Thread Cotter, Colin J via petsc-users
Hi Dave, >If you are doing that, then you need to tell fieldsplit to use the Amat to >define the splits otherwise it will define the Schur compliment as >S = B22 - B21 inv(B11) B12 >preconditiones with B22, where as what you want is >S = A22 - A21 inv(A11) A12 >preconditioned with B22. >If your

[petsc-users] Confusing Schur preconditioner behaviour

2019-03-18 Thread Cotter, Colin J via petsc-users
Dear petsc-users, I'm solving a 2x2 block system, for which I can construct the Schur complement analytically (through compatible FEM stuff), which I can pass as the preconditioning matrix. When using gmres on the outer iteration, and preonly+lu on the inner iterations with a Schur

Re: [petsc-users] SLEPc: Setting options for mass matrix KSP in generalised eigenvalue problem

2019-01-16 Thread Cotter, Colin J via petsc-users
ype.SHIFT) instead of F.setType(F.Type.PRECOND). Jose > El 16 ene 2019, a las 10:32, Cotter, Colin J via petsc-users > escribió: > > Dear PETSc mailing list, > A SLEPc question, sorry if this is obvious and I missed something. > > How do I set the KSP for the ma

[petsc-users] SLEPc: Setting options for mass matrix KSP in generalised eigenvalue problem

2019-01-16 Thread Cotter, Colin J via petsc-users
Dear PETSc mailing list, A SLEPc question, sorry if this is obvious and I missed something. How do I set the KSP for the mass matrix M in the generalised eigenvalue problem (I'm looking for the eigenvalue with largest magnitude) Ax = lambda*Mx ? I can only find how to set the KSP for a