[petsc-dev] PetscLayoutSetWeights?

[Jed Brown]

Barry Smith <bsmith at mcs.anl.gov> writes:
>    So we then freely change the thread offsets in the layout based on
>    the matrix entries? 

I would prefer to only allow it before PetscLayoutSetUp.

>    This means the Vec and Mat Layouts have to be same? This is true if
>    the Vecs are obtained with MatGetVecs() after the matrix is created
>    but not if the Vecs are created independently. Of for the MatMult
>    do you just use the thread offsets from the Mat layout and ignore
>    the Vec layout if it is different?

If they are different, I would ignore the Vec layout.  (Favor load
balance and NUMA distribution of the matrix over NUMA distribution of
the Vec.)
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: application/pgp-signature
Size: 835 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20130923/a4d60d31/attachment.sig>