Re: [petsc-dev] sor smoothers

2013-09-19 Thread Jed Brown
Mark F. Adams mfad...@lbl.gov writes: Jed, where you going to give me a PetscObjectGetId? It is in 'jed/object-id', which is now in 'next'. pgp116CD7jI0C.pgp Description: PGP signature

Re: [petsc-dev] coarse grid in gamg (for the OPTE paper)

2013-09-19 Thread Jungho Lee
To fill others in on the situation here; this was regarding a variational inequality example, that uses reduced space method for which routines are defined in src/snes/impls/vi/rs/virs.c. When MatGetSubMatrix was called to extract rows and columns of the original matrix corresponding to the

Re: [petsc-dev] coarse grid in gamg (for the OPTE paper)

2013-09-19 Thread Dmitry Karpeyev
On Thu, Sep 19, 2013 at 1:28 PM, Jungho Lee ju...@mcs.anl.gov wrote: To fill others in on the situation here; this was regarding a variational inequality example, that uses reduced space method for which routines are defined in src/snes/impls/vi/rs/virs.c. When MatGetSubMatrix was called to

Re: [petsc-dev] coarse grid in gamg (for the OPTE paper)

2013-09-19 Thread Barry Smith
On Sep 19, 2013, at 3:17 PM, Dmitry Karpeyev dkarp...@gmail.com wrote: On Thu, Sep 19, 2013 at 1:28 PM, Jungho Lee ju...@mcs.anl.gov wrote: To fill others in on the situation here; this was regarding a variational inequality example, that uses reduced space method for which routines

Re: [petsc-dev] coarse grid in gamg (for the OPTE paper)

2013-09-19 Thread Jed Brown
Barry Smith bsm...@mcs.anl.gov writes: That is what it is doing and apparently it doesn't result in a good preconditioner; I don't know why off hand. One thing is it no longer knows about the block structure. How is the near null space being specified? pgpnJypvyI__v.pgp Description:

Re: [petsc-dev] coarse grid in gamg (for the OPTE paper)

2013-09-19 Thread Mark F. Adams
On Sep 19, 2013, at 3:53 PM, Barry Smith bsm...@mcs.anl.gov wrote: It is essentially 3 Laplacians If you mean literally 3 scalar laplacians packed in one matrix for some reason then the one constant vector is fine. Is this a block diagonal matrix with 3 big blocks, if you order it

Re: [petsc-dev] coarse grid in gamg (for the OPTE paper)

2013-09-19 Thread Barry Smith
It is essentially 3 Laplacians so I think the default null space of 3 constant vectors is fine. The problem is without the block information presumably GAMG is using only a single constant vector over all variables? So maybe we need to construct a 3 vector null space which just marks in

Re: [petsc-dev] coarse grid in gamg (for the OPTE paper)

2013-09-19 Thread Jed Brown
Barry Smith bsm...@mcs.anl.gov writes: It is essentially 3 Laplacians so I think the default null space of 3 constant vectors is fine. The problem is without the block information presumably GAMG is using only a single constant vector over all variables? So maybe we need to

Re: [petsc-dev] coarse grid in gamg (for the OPTE paper)

2013-09-19 Thread Dmitry Karpeyev
On Thu, Sep 19, 2013 at 2:53 PM, Barry Smith bsm...@mcs.anl.gov wrote: It is essentially 3 Laplacians so I think the default null space of 3 constant vectors is fine. The problem is without the block information presumably GAMG is using only a single constant vector over all variables?

[petsc-dev] Push restrictions on branches

2013-09-19 Thread Jed Brown
Bitbucket added support to restrict push access based on branch names (glob matching). For example, that would allow us to have a smaller group of people with access to merge to 'maint' or 'master'. Is this a feature we should start using in petsc.git? One tangible difference from the current

Re: [petsc-dev] coarse grid in gamg (for the OPTE paper)

2013-09-19 Thread Barry Smith
On Sep 19, 2013, at 4:28 PM, Mark F. Adams mfad...@lbl.gov wrote: On Sep 19, 2013, at 3:53 PM, Barry Smith bsm...@mcs.anl.gov wrote: It is essentially 3 Laplacians If you mean literally 3 scalar laplacians packed in one matrix for some reason then the one constant vector is fine.

Re: [petsc-dev] coarse grid in gamg (for the OPTE paper)

2013-09-19 Thread Jed Brown
Barry Smith bsm...@mcs.anl.gov writes: If you mean literally 3 scalar laplacians packed in one matrix for some reason then the one constant vector is fine. Really, but (for periodic or N. bc) the null space has three vectors, you don't need all of them in your near null space? If the

Re: [petsc-dev] Push restrictions on branches

2013-09-19 Thread Patrick Sanan
This sounds like a great idea from the perspective of someone like me, who would only be pushing to maint/master/next in error at this point. It would also allow slightly easier access to features which are new/experimental enough to not be in next yet. On Sep 19, 2013, at 4:35 PM, Jed Brown

Re: [petsc-dev] coarse grid in gamg (for the OPTE paper)

2013-09-19 Thread Mark F. Adams
On Sep 19, 2013, at 7:05 PM, Barry Smith bsm...@mcs.anl.gov wrote: On Sep 19, 2013, at 4:28 PM, Mark F. Adams mfad...@lbl.gov wrote: On Sep 19, 2013, at 3:53 PM, Barry Smith bsm...@mcs.anl.gov wrote: It is essentially 3 Laplacians If you mean literally 3 scalar laplacians packed