[petsc-dev] Mat Type: HYPRE

[Kong, Fande]

On Wed, Oct 4, 2017 at 8:07 AM, Stefano Zampini <stefano.zampini at gmail.com>
wrote:

>
>
> 2017-10-04 16:54 GMT+03:00 Kong, Fande <fande.kong at inl.gov>:
>
>> Hi Stefano,
>>
>> Thanks for your email. The error arises from the first assembly.
>>
>
> That's strange. Can you provide a minimal working example to reproduce the
> problem?
>

I am working on MOOSE. I will try a native example to reproduce this issue.


>
>> 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.
>>
>
> Take a look at src/mat/examples/tests/ex115.c
>

Thanks,

Fande,


>
>
>>
>> Thanks,
>>
>> Fande,
>>
>>
>> On Wed, Oct 4, 2017 at 2:27 AM, Stefano Zampini <
>> stefano.zampini at gmail.com> wrote:
>>
>>> Does this error arise at the first assembly of the matrix or in a
>>> successive assembly?
>>>
>>> In the second case, I guess the issue is with HYPRE, or the way we are
>>> passing the preallocation to HYPRE routines.
>>> For some reasons that I'm not aware of, reassembling an HypreParCSR
>>> matrix may have issues.
>>> My knowledge of the HYPRE code is not so deep to tell you why.
>>>
>>> 2017-10-04 1:03 GMT+03:00 Kong, Fande <fande.kong at inl.gov>:
>>>
>>>> Hi ALL,
>>>>
>>>> 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.
>>>>
>>>> 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.
>>>>
>>>>
>>>> [0]PETSC ERROR: --------------------- Error Message
>>>> --------------------------------------------------------------
>>>> [0]PETSC ERROR: Error in external library [0]PETSC ERROR: Error in
>>>> HYPRE_IJMatrixAddToValues(): error code 1 [0]PETSC ERROR: See
>>>> http://www.mcs.anl.gov/petsc/documentation/faq.html
>>>> <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=>
>>>> 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-master [0]PETSC
>>>> ERROR: #1 MatSetValues_HYPRE() line 1014 in /Users/kongf/projects/petsc/src/mat/impls/hypre/mhypre.c
>>>> [0]PETSC ERROR: #2 MatSetValues() line 1270 in
>>>> /Users/kongf/projects/petsc/src/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
>>>>
>>>>
>>>> Fande,
>>>>
>>>
>>>
>>>
>>> --
>>> Stefano
>>>
>>
>>
>
>
> --
> Stefano
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20171004/8ff6e7be/attachment.html>