[petsc-dev] array ownership of MatCreateMPIAIJWithSplitArrays and MatCreateSeqAIJWithArrays
stefano.zampini at gmail.com
Sat Jul 23 14:36:31 CDT 2016
I’m interfacing a FEM library with PETSc.
They already use Hypre and so they have a valid distributed CSR format (except for the diagonal entry of the diagonal block :-))
For efficiency reasons, I was thinking of using MatCreateMPIAIJWithSplitArrays and MatCreateSeqAIJWithArrays;
However, If using these functions as they are now I will need to keep track of the allocated arrays, even with referencing objects, and this could be annoying and cumbersome.
A trivial fix would be to add an extra PetscCopyMode parameter to these routines, so that PETSc can take ownership of the arrays.
Thoughts? Any objection?
More information about the petsc-dev