[petsc-users] Mat_AllocAIJ_CSR in petsc4py

[D D]

Why does MatAlloc_AIJ_CSR make two preallocation calls?

        CHKERR( MatSeqAIJSetPreallocationCSR(A, i, j, v) )
        CHKERR( MatMPIAIJSetPreallocationCSR(A, i, j, v) )

I am assuming memory if preallocated twice, once for sequential and another for MPI. So each matrix created with createAIJ will have a sequential and MPI structure.
Is this for convenience when switching from sequential to MPI?
Thanks,Dale
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20170426/e6d64aff/attachment.html>