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?
