Re: [petsc-dev] Copy scatter and change block size

2017-11-25 Thread Pierre Jolivet
> On 25 Nov 2017, at 7:49 PM, Stefano Zampini wrote: > > Pierre, > > If you have access to the indices used to create the "unblocked" vecscatter, > you can then call ISCreateBlock and obtain the vecscatter you want. OK, that’s a better start. > To obtain the

Re: [petsc-dev] Copy scatter and change block size

2017-11-25 Thread Stefano Zampini
Pierre, If you have access to the indices used to create the "unblocked" vecscatter, you can then call ISCreateBlock and obtain the vecscatter you want. To obtain the garray, MatMPIAIJGetSeqAIJ > On Nov 25, 2017, at 9:26 PM, Pierre Jolivet > wrote: > > OK, so I

Re: [petsc-dev] Copy scatter and change block size

2017-11-25 Thread Pierre Jolivet
OK, so I guess I need to tinker around with VecScatterCreate. As a follow-up question, is there an easy way to get the “garray” from http://www.mcs.anl.gov/petsc/petsc-current/src/mat/impls/aij/mpi/mmaij.c.html#line105

Re: [petsc-dev] Copy scatter and change block size

2017-11-25 Thread Smith, Barry F.
Pierre, There is not. The block size is baked directly into the scatter data structure so not trivial to change in a copy. Barry > On Nov 25, 2017, at 7:41 AM, Pierre Jolivet > wrote: > > Hello, > Given a VecScatter inctx that operates on vectors x and

[petsc-dev] Copy scatter and change block size

2017-11-25 Thread Pierre Jolivet
Hello, Given a VecScatter inctx that operates on vectors x and y, I’d like to create a VecScatter inctx_bs with the same pattern, but that operators on vectors x_bs and y_bs with different block sizes than the original x and y (and that are thus of sizes n * bs). Is there an easy way to do this

Re: [petsc-dev] composite vs. shell

2017-11-25 Thread Pierre Jolivet
Thank you Barry. I made the appropriate changes in my code but it’s good to know about the new behaviour. Thanks, Pierre > On 25 Nov 2017, at 2:34 PM, Smith, Barry F. wrote: > > > Pierre, > >

Re: [petsc-dev] composite vs. shell

2017-11-25 Thread Smith, Barry F.
Pierre, https://bitbucket.org/petsc/petsc/pull-requests/809/only-propagate-operators-into-inner-pcs-in/diff Sorry for the long delay in responding. Barry > On Nov 10, 2017, at 1:36 PM, Jed Brown wrote: > > "Smith, Barry F." writes: > >>> On Nov

Re: [petsc-dev] some change in DM plex output breaks next?

2017-11-25 Thread Lisandro Dalcin
Sorry about that, I updated PR 806 and pushed the fix to next. On 25 November 2017 at 03:01, Smith, Barry F. wrote: > > ftp://ftp.mcs.anl.gov/pub/petsc/nightlylogs/archive/2017/11/24/examples_full_next.log > > not ok diff-dm_impls_plex_tests-ex1_gmsh_6 > # 8d7 > #