[petsc-dev] assuming sizeof(PetscInt) <= sizeof(PetscScalar)

Matthew Knepley knepley at gmail.com
Sat Feb 8 17:52:26 CST 2014


On Sat, Feb 8, 2014 at 5:40 PM, Barry Smith <bsmith at mcs.anl.gov> wrote:

>
>    It appears people have been assuming sizeof(PetscInt) <=
> sizeof(PetscScalar) right and left in the code (to take advantage of
> VecScatter functionality). This is most definitely not true with
> --with-64-bit-indices --with-precision=single
>
> MatGetBrowsOfAoCols_MPIAIJ
> PCGAMGProlongator_Classical_Standard
> maybe others
>
> VecScatters are for moving Vec elements around, they are not a replacement
> for MPI calls!
>
> Please don't do this in the future.


Please please please use PetscSF if you have regular layouts, and
PetscSection for irregular layouts.

  Matt


>
>    Barry
>
>


-- 
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20140208/a7214d19/attachment.html>


More information about the petsc-dev mailing list