[petsc-users] Error while calling a fieldsplit preconditioner containing a mg/gamg preconditioner

2014-11-21 Thread Luc Berger-Vergiat
Hi all, I am using a DMShell to create to use a fieldsplit preconditioner. I would like to try some of petsc's multigrid options: mg and gamg. When I call my preconditioner: -ksp_type gmres -pc_type fieldsplit -pc_fieldsplit_type schur -pc_fieldsplit_schur_factorization_type full

Re: [petsc-users] Error while calling a fieldsplit preconditioner containing a mg/gamg preconditioner

2014-11-21 Thread Matthew Knepley
On Fri, Nov 21, 2014 at 9:55 AM, Luc Berger-Vergiat lb2...@columbia.edu wrote: Hi all, I am using a DMShell to create to use a fieldsplit preconditioner. I would like to try some of petsc's multigrid options: mg and gamg. Lets start with gamg, since that does not need anything else from the

Re: [petsc-users] Error while calling a fieldsplit preconditioner containing a mg/gamg preconditioner

2014-11-21 Thread Luc Berger-Vergiat
Sure, so with gamg (which was my first choice too), I get the following: [0]PETSC ERROR: - Error Message -- [0]PETSC ERROR: [0]PETSC ERROR: Must call DMShellSetGlobalVector() or DMShellSetCreateGlobalVector()

Re: [petsc-users] Error while calling a fieldsplit preconditioner containing a mg/gamg preconditioner

2014-11-21 Thread Lawrence Mitchell
On 21 Nov 2014, at 16:55, Luc Berger-Vergiat lb2...@columbia.edu wrote: Sure, so with gamg (which was my first choice too), I get the following: [0]PETSC ERROR: - Error Message -- [0]PETSC ERROR:

Re: [petsc-users] Error while calling a fieldsplit preconditioner containing a mg/gamg preconditioner

2014-11-21 Thread Luc Berger-Vergiat
Yes good call, I pulled master and got the patch. I don't have errors anymore. Do this mean that when a dm is active, the dm is assumed to create hold the information about the vectors (residual, solution and internal work vectors) for the problem, whereas when the dm is not active or is not

Re: [petsc-users] Error while calling a fieldsplit preconditioner containing a mg/gamg preconditioner

2014-11-21 Thread Matthew Knepley
On Fri, Nov 21, 2014 at 3:12 PM, Luc Berger-Vergiat lb2...@columbia.edu wrote: Yes good call, I pulled master and got the patch. I don't have errors anymore. +1 Lawrence Do this mean that when a dm is active, the dm is assumed to create hold the information about the vectors (residual,