Stefeno, Doe this look right to you? Thanks Barry

On Dec 17, 2018, at 3:18 AM, Pierre Jolivet via petsc-dev <[email protected]> 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>