Re: [petsc-users] multigrid preconditioning and adaptivity

2016-03-15 Thread Lukasz Kaczmarczyk
Hello Matt, On 7 Mar 2016, at 16:42, Matthew Knepley > wrote: Yes, that should be it. It would be nice to have some example that does this if you would be willing to contribute some version of your code. Following discussion from last week, I had

Re: [petsc-users] multigrid preconditioning and adaptivity

2016-03-07 Thread Lukasz Kaczmarczyk
On 7 Mar 2016, at 16:42, Matthew Knepley > wrote: On Mon, Mar 7, 2016 at 10:28 AM, Lukasz Kaczmarczyk > wrote: Many thanks all for help. I started to implement function for

Re: [petsc-users] multigrid preconditioning and adaptivity

2016-03-07 Thread Matthew Knepley
On Mon, Mar 7, 2016 at 10:28 AM, Lukasz Kaczmarczyk < lukasz.kaczmarc...@glasgow.ac.uk> wrote: > Many thanks all for help. I started to implement function for DM. > > I understand that minimal implementation is that for the DM i need to > have, is to have DMCoarsen and in each level for all DMs,

Re: [petsc-users] multigrid preconditioning and adaptivity

2016-03-07 Thread Lukasz Kaczmarczyk
Many thanks all for help. I started to implement function for DM. I understand that minimal implementation is that for the DM i need to have, is to have DMCoarsen and in each level for all DMs, set operators DMKSPSetComputeOperators and DMCreateInterpolation. Matrix matrix free P from

Re: [petsc-users] multigrid preconditioning and adaptivity

2016-03-07 Thread Mark Adams
You can just set the coarse grid matrix/operator instead of using Galerkin. If you have a shell (matrix free) P then you will need to create and set this yourself. Our Galerkin requires a matrix P. On Mon, Mar 7, 2016 at 9:32 AM, Lukasz Kaczmarczyk < lukasz.kaczmarc...@glasgow.ac.uk> wrote: >

Re: [petsc-users] multigrid preconditioning and adaptivity

2016-03-07 Thread Lukasz Kaczmarczyk
> On 7 Mar 2016, at 14:21, Lawrence Mitchell > wrote: > > On 07/03/16 14:16, Lukasz Kaczmarczyk wrote: >> >>> On 7 Mar 2016, at 13:50, Matthew Knepley >> > wrote: >>> >>> On Mon, Mar 7, 2016 at 6:58 AM, Lukasz

Re: [petsc-users] multigrid preconditioning and adaptivity

2016-03-07 Thread Lawrence Mitchell
On 07/03/16 14:16, Lukasz Kaczmarczyk wrote: > >> On 7 Mar 2016, at 13:50, Matthew Knepley > > wrote: >> >> On Mon, Mar 7, 2016 at 6:58 AM, Lukasz Kaczmarczyk >> > > wrote: >>

Re: [petsc-users] multigrid preconditioning and adaptivity

2016-03-07 Thread Lukasz Kaczmarczyk
On 7 Mar 2016, at 13:50, Matthew Knepley > wrote: On Mon, Mar 7, 2016 at 6:58 AM, Lukasz Kaczmarczyk > wrote: Hello, I run multi-grid solver, with adaptivity, works well,

Re: [petsc-users] multigrid preconditioning and adaptivity

2016-03-07 Thread Matthew Knepley
On Mon, Mar 7, 2016 at 6:58 AM, Lukasz Kaczmarczyk < lukasz.kaczmarc...@glasgow.ac.uk> wrote: > Hello, > > I run multi-grid solver, with adaptivity, works well, however It is some > space for improving efficiency. I using hierarchical approximation basis, > for which > construction of

Re: [petsc-users] multigrid preconditioning and adaptivity

2016-03-07 Thread Lukasz Kaczmarczyk
> On 7 Mar 2016, at 12:58, Lukasz Kaczmarczyk > wrote: > > Hello, > > I run multi-grid solver, with adaptivity, works well, however It is some > space for improving efficiency. I using hierarchical approximation basis, for > which > construction of

[petsc-users] multigrid preconditioning and adaptivity

2016-03-07 Thread Lukasz Kaczmarczyk
Hello, I run multi-grid solver, with adaptivity, works well, however It is some space for improving efficiency. I using hierarchical approximation basis, for which construction of interpolation operators is simple, it is simple injection. After each refinement level (increase of order of