> 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
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
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
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
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
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,
>
>
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
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
> #