Stefeno, Doe this look right to you?
Thanks Barry > On Dec 17, 2018, at 3:18 AM, Pierre Jolivet via petsc-dev > <petsc-dev@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>