<div dir="ltr"><div>Hi Stefano,</div><div><br></div><div>Thanks for your email. The error arises from the first assembly. </div><div><br></div><div>Do you have any example demonstrating how to use HYPRE mat in PETSc? I would like to play around, and fix any issue I encounter.</div><div><br></div><div>Thanks,</div><div><br></div><div>Fande, </div><br><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Oct 4, 2017 at 2:27 AM, Stefano Zampini <span dir="ltr"><<a href="mailto:stefano.zampini@gmail.com" target="_blank">stefano.zampini@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div dir="ltr">Does this error arise at the first assembly of the matrix or in a successive assembly?<div><br><div>In the second case, I guess the issue is with HYPRE, or the way we are passing the preallocation to HYPRE routines.</div><div>For some reasons that I'm not aware of, reassembling an HypreParCSR matrix may have issues.</div><div>My knowledge of the HYPRE code is not so deep to tell you why.</div></div></div><div class="gmail_extra"><div><div class="gmail-h5"><br><div class="gmail_quote">2017-10-04 1:03 GMT+03:00 Kong, Fande <span dir="ltr"><<a href="mailto:fande.kong@inl.gov" target="_blank">fande.kong@inl.gov</a>></span>:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div dir="ltr">Hi ALL,<div><br></div><div>There is a new Mat type: HYPRE in PETSc. I would like to try it out because it may save the memory if use a hypre preconditioner. When use a native PETSc type such as MPIAIJ, we need to keep an extra copy of Jacobian matrix usually for the hypre preconditioner. </div><div><br></div><div>I was wondering if or not Mat HYPRE is complete. There are some errors when I am tring to use it. The code works fine with MPIAIJ.</div><div><br></div><div><div><br></div><div>[0]PETSC ERROR: --------------------- Error Message ------------------------------<wbr>------------------------------<wbr>-- [0]PETSC ERROR: Error in external library [0]PETSC ERROR: Error in HYPRE_IJMatrixAddToValues(): error code 1 [0]PETSC ERROR: See <a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__www.mcs.anl.gov_petsc_documentation_faq.html&d=DwMFaQ&c=54IZrppPQZKX9mLzcGdPfFD1hxrcB__aEkJFOKJFd00&r=DUUt3SRGI0_JgtNaS3udV68GRkgV4ts7XKfj2opmiCY&m=E9mqRwrQ3kax5iMptBq__1jUmNZ3uujcOM1aXkMH9WI&s=sAvCaFdkzfe1gCJgHo1PgCvSYrVbl5iRbweHKpysIdk&e=" target="_blank">http://www.mcs.anl.gov/petsc/d<wbr>ocumentation/faq.html</a> for trouble shooting. [0]PETSC ERROR: Petsc Development GIT revision: v3.8-2-gd4ad237  GIT Date: 2017-10-03 12:37:25 -0600 [0]PETSC ERROR: ../../../moose_test-opt on a arch-darwin-c-debug-master named FN604208 by kongf Tue Oct  3 15:54:15 2017 [0]PETSC ERROR: Configure options --with-debugging=yes --with-shared-libraries=1 --with-mpi=1 --download-fblaslapack=1 --download-metis=1 --download-parmetis=1 --download-superlu_dist=1 --download-hypre=1 --download-mumps=1 --download-scalapack=1 PETSC_ARCH=arch-darwin-c-debug<wbr>-master [0]PETSC ERROR: #1 MatSetValues_HYPRE() line 1014 in /Users/kongf/projects/petsc/sr<wbr>c/mat/impls/hypre/mhypre.c [0]PETSC ERROR: #2 MatSetValues() line 1270 in /Users/kongf/projects/petsc/sr<wbr>c/mat/interface/matrix.c application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0 [unset]: write_line error; fd=-1 buf=:cmd=abort exitcode=1 : system msg for write_line failure : Bad file descriptor</div></div><span class="gmail-m_-6958367300722685860HOEnZb"><font color="#888888"><div><br></div><div><br></div><div>Fande,</div></font></span></div>
</blockquote></div><br><br clear="all"><div><br></div></div></div><span class="gmail-HOEnZb"><font color="#888888">-- <br><div class="gmail-m_-6958367300722685860gmail_signature">Stefano</div>
</font></span></div>
</blockquote></div><br></div></div>