Re: [petsc-users] FieldSplit, multigrid and blocksize

Hi Nicolas, for us the solution has been to "manually" create a MatNest with, ie, block A00 containing only u-u coupling and block A11 containing p-p coupling. Thus, we are able to assign block size of 3 for A00 and block size of 1 for A11. The other thing we did, is to be able to number

Re: [petsc-users] FieldSplit, multigrid and blocksize

Thank you all. These are the answers I was looking for! Best regards, Nicolas 2016-12-07 15:08 GMT+01:00 Mark Adams : > Note, for best performance with ML and GAMG you want to give it the near > kernel for the 00 block. These are the 6 "rigid body modes" or zero energy > modes.

Re: [petsc-users] FieldSplit, multigrid and blocksize

Note, for best performance with ML and GAMG you want to give it the near kernel for the 00 block. These are the 6 "rigid body modes" or zero energy modes. PETSc provides some tools to do that (eg, MatNullSpaceCreateRigidBody). On Wed, Dec 7, 2016 at 4:45 PM, Lawrence Mitchell <

Re: [petsc-users] FieldSplit, multigrid and blocksize

On 07/12/16 13:43, Karin wrote: > Thanks Barry. > I must emphasize that my unknowns are not numbered in a regular way : > I am using a P2-P1 finite element and the middle nodes do not carry a > pressure DOF. So the global numbering is somewhat like : >

Re: [petsc-users] FieldSplit, multigrid and blocksize

Thanks Barry. I must emphasize that my unknowns are not numbered in a regular way : I am using a P2-P1 finite element and the middle nodes do not carry a pressure DOF. So the global numbering is somewhat like :

Re: [petsc-users] FieldSplit, multigrid and blocksize

> On Dec 7, 2016, at 7:06 AM, Karin wrote: > > Dear PETSc gurus, > > I am using FieldSplit to solve a poro-mechanics problem. Thus, I am dealing > with 3 displacement DOF and 1 pressure DOF. > In order to precondition the 00 block (aka the displacement block), I am >

[petsc-users] FieldSplit, multigrid and blocksize

Dear PETSc gurus, I am using FieldSplit to solve a poro-mechanics problem. Thus, I am dealing with 3 displacement DOF and 1 pressure DOF. In order to precondition the 00 block (aka the displacement block), I am using a multigrid method (ml or gamg). Nevertheless, I have the feeling that the