
<html>

  <head>


    <meta http-equiv="content-type" content="text/html; charset=ISO-8859-1">

  </head>

  <body bgcolor="#FFFFFF" text="#000000">

    Dear petsc devs,<br>

    <br>

    I have came across strange behavior trying to make use of PETSc

    routines (from mhyp.c) to convert <i>rectangular</i> matrices into

    hypre format in parallel (i.e. MATMPIAIJ). <br>

    Whenever I use MatHYPRE_IJMatrixCreate -> MatHYPRE_IJMatrixCopy

    -> hypre_ParCSRMatrixTranspose I get following error:<br>

    hypre_assert failed: AT_i[A_j[j]] < num_nonzerosAT<br>

    hypre_assert failed: AT_i[A_j[j]] < num_nonzerosAT<br>

    <br>

    it comes from csr_matop.c : hypre_CSRMatrixTranspose.<br>

    <br>

    The error message never appears if construct hypre matrix directly

    and transpose it. <br>

    <br>

    Attached is a minimal example that reproduces the problem (try to

    (un)comment line 58).<br>

    <br>

    I am wondering where is an origin of this error: in my code, in

    petsc or in hypre?<br>

    <br>

    Thanks,<br>

    Alexander<br>

  </body>

</html>