message: adding the function call 'hypre_ParCSRMatrixSetNumNonzeros' so that the returned matrix has the right number of non-zeros (NNZ)
diff: @@ -1376,6 +1376,7 @@ HypreParMatrix * ParMult(HypreParMatrix *A, HypreParMatrix *B)
  {
  hypre_ParCSRMatrix * ab;
  ab = hypre_ParMatmul(*A,*B);
  + hypre_ParCSRMatrixSetNumNonzeros(ab);
  hypre_MatvecCommPkgCreate(ab);
  