Re: [petsc-users] grid sequencing & DMComposite

2015-09-11 Thread Matthew Knepley
On Fri, Sep 11, 2015 at 8:15 AM, Gideon Simpson wrote: > Ok, so here’s some behavior I don’t understand. I am working with a > DMComposite structure, and I do the following, > > SNESComputeFunction(snes,U_refine,r); > SNESGetDM(snes, _refine); > >

Re: [petsc-users] grid sequencing & DMComposite

2015-09-11 Thread Matthew Knepley
On Fri, Sep 11, 2015 at 8:37 AM, Gideon Simpson wrote: > There doesn’t seem to be a global analog of DMCompositeGetLocalVectors. > Do I need to do manual indexing? > Its simple, lets just write it:

Re: [petsc-users] grid sequencing & DMComposite

2015-09-11 Thread Gideon Simpson
There doesn’t seem to be a global analog of DMCompositeGetLocalVectors. Do I need to do manual indexing? -gideon > On Sep 11, 2015, at 9:36 AM, Matthew Knepley wrote: > > On Fri, Sep 11, 2015 at 8:15 AM, Gideon Simpson

Re: [petsc-users] grid sequencing & DMComposite

2015-09-11 Thread Lawrence Mitchell
> On 11 Sep 2015, at 14:37, Gideon Simpson wrote: > > There doesn’t seem to be a global analog of DMCompositeGetLocalVectors. Do I > need to do manual indexing? I think you use DMCompositeGetAccess. Lawrence signature.asc Description: Message signed with OpenPGP

Re: [petsc-users] grid sequencing & DMComposite

2015-09-11 Thread Matthew Knepley
On Fri, Sep 11, 2015 at 8:48 AM, Lawrence Mitchell < lawrence.mitch...@imperial.ac.uk> wrote: > > > On 11 Sep 2015, at 14:37, Gideon Simpson > wrote: > > > > There doesn’t seem to be a global analog of DMCompositeGetLocalVectors. > Do I need to do manual indexing? > > I

Re: [petsc-users] grid sequencing & DMComposite

2015-09-11 Thread Gideon Simpson
Yea, that works as advertised. -gideon > On Sep 11, 2015, at 9:48 AM, Lawrence Mitchell > wrote: > > >> On 11 Sep 2015, at 14:37, Gideon Simpson wrote: >> >> There doesn’t seem to be a global analog of DMCompositeGetLocalVectors.

Re: [petsc-users] grid sequencing & DMComposite

2015-09-11 Thread Gideon Simpson
Ok, so here’s some behavior I don’t understand. I am working with a DMComposite structure, and I do the following, SNESComputeFunction(snes,U_refine,r); SNESGetDM(snes, _refine); DMCompositeGetLocalVectors(dm_refine, , ); DMCompositeScatter(dm_refine, r, rp, rQ);

[petsc-users] grid sequencing & DMComposite

2015-09-10 Thread Gideon Simpson
I’m using a DMCompsosite along with grid sequencing and here’s what I’m trying to accomplish. After running the SNES solve, I’d like to evaluate the residual on the refined grid, on each piece of the DMComposite. How do I get a DM for the refined grid which I can then use with

Re: [petsc-users] grid sequencing & DMComposite

2015-09-10 Thread Barry Smith
SNESGetDM() this will return not the original DM you set but the refined one. > On Sep 10, 2015, at 3:47 PM, Gideon Simpson wrote: > > I’m using a DMCompsosite along with grid sequencing and here’s what I’m > trying to accomplish. After running the SNES solve,