[petsc-dev] Faulty logic in MatHYPRESetPreallocation?

[Smith, Barry F.]

  Pierre,

    Could you make a pull request on this?

    Thanks

     Barry


> On Dec 17, 2018, at 3:18 AM, Pierre Jolivet via petsc-dev <petsc-dev at mcs.anl.gov> wrote:
> 
> Hello,
> I am not sure the hypre routines are called in the right order in MatHYPRESetPreallocation.
> Indeed, attached, you’ll find two Massif logs, one with the current order (massif.out.bad), and the other (massif.out.better) with the attached patch applied.
> In the former, you can see at the snapshot #31 that there is a call to hypre_CTAlloc in hypre_AuxParCSRMatrixInitialize (aux_parcsr_matrix.c:163).
> This means that hypre is allocating memory as if it was not aware of the number of nonzero per row, even though I supply in my application the two arrays dnnz and onnz.
> 
> Am I missing something?
> 
> Thanks,
> Pierre
> 
> <massif.out.bad><massif.out.better><patch-hypre.txt>