[petsc-users] About MatMumpsSetIcntl function

[Matthew Knepley]

On Wed, Nov 30, 2022 at 8:58 AM 김성익 <ksi2443 at gmail.com> wrote:

> I'm working on FEM.
>   When I used mumps alone, I fount it efficient to use mumps with metis.
> So my purpose is using MUMPSsolver with METIS.
>
> I tried to set metis (by icntl_7 : 5) after global matrix assembly and
> just before kspsolve.
> However there is error because of 'pcfactorgetmatrix' and
> 'matmumpsseticntl'.
>
> How can I fix this?
>

Give the Icntrl as an option.

  Thanks,

     Matt


> Thanks,
> Hyung Kim
>
> 2022년 11월 30일 (수) 오후 10:44, Matthew Knepley <knepley at gmail.com>님이 작성:
>
>> On Wed, Nov 30, 2022 at 8:40 AM 김성익 <ksi2443 at gmail.com> wrote:
>>
>>> Following your comments,
>>>
>>> After matrix assembly end,
>>> PetscCall(KSPGetPC(ksp,&pc));
>>> PetscCall(KSPSetFromOptions(ksp));
>>> PetscCall(KSPSetUp(ksp));
>>> PetscCall(PCFactorGetMatrix(pc,&xGK));
>>>
>>> However there is another error as below.
>>> [0]PETSC ERROR: Object is in wrong state
>>> [0]PETSC ERROR: Not for factored matrix
>>>
>>
>> The error message is telling you that you cannot alter values in the
>> factored matrix. This is because
>> the direct solvers use their own internal storage formats which we cannot
>> alter, and you should probably
>> not alter either.
>>
>> What are you trying to do?
>>
>>   Thanks,
>>
>>      Matt
>>
>>
>>> [0]PETSC ERROR: See https://petsc.org/release/faq/ for trouble shooting.
>>> [0]PETSC ERROR: Petsc Release Version 3.18.1, unknown
>>> [0]PETSC ERROR: ./app on a arch-linux-c-debug named ubuntu by ksi2443
>>> Wed Nov 30 05:37:52 2022
>>> [0]PETSC ERROR: Configure options -download-mumps -download-scalapack
>>> -download-parmetis -download-metis
>>> [0]PETSC ERROR: #1 MatZeroEntries() at
>>> /home/ksi2443/petsc/src/mat/interface/matrix.c:6024
>>> [0]PETSC ERROR: #2 main() at /home/ksi2443/Downloads/coding/a1.c:339
>>> [0]PETSC ERROR: No PETSc Option Table entries
>>>
>>> How can I fix this?
>>>
>>>
>>> Thanks,
>>> Hyung Kim
>>>
>>>
>>> 2022년 11월 30일 (수) 오후 4:18, Jose E. Roman <jroman at dsic.upv.es>님이 작성:
>>>
>>>> You have to call PCFactorGetMatrix() first. See any of the examples
>>>> that use MatMumpsSetIcntl(), for instance
>>>> https://petsc.org/release/src/ksp/ksp/tutorials/ex52.c.html
>>>>
>>>> Jose
>>>>
>>>>
>>>> > El 30 nov 2022, a las 6:52, 김성익 <ksi2443 at gmail.com> escribió:
>>>> >
>>>> > Hello,
>>>> >
>>>> >
>>>> > I tried to adopt METIS option in MUMPS by using
>>>> > ' PetscCall(MatMumpsSetIcntl(Mat, 7, 5));'
>>>> >
>>>> > However, there is an error as follows
>>>> >
>>>> > [0]PETSC ERROR: Object is in wrong state
>>>> > [0]PETSC ERROR: Only for factored matrix
>>>> > [0]PETSC ERROR: See https://petsc.org/release/faq/ for trouble
>>>> shooting.
>>>> > [0]PETSC ERROR: Petsc Release Version 3.18.1, unknown
>>>> > [0]PETSC ERROR: ./app on a arch-linux-c-debug named ubuntu by ksi2443
>>>> Tue Nov 29 21:12:41 2022
>>>> > [0]PETSC ERROR: Configure options -download-mumps -download-scalapack
>>>> -download-parmetis -download-metis
>>>> > [0]PETSC ERROR: #1 MatMumpsSetIcntl() at
>>>> /home/ksi2443/petsc/src/mat/impls/aij/mpi/mumps/mumps.c:2478
>>>> > [0]PETSC ERROR: #2 main() at /home/ksi2443/Downloads/coding/a1.c:149
>>>> > [0]PETSC ERROR: No PETSc Option Table entries
>>>> >
>>>> > How can I fix this error?
>>>> >
>>>> > Thank you for your help.
>>>> >
>>>> >
>>>> > Hyung Kim
>>>>
>>>>
>>
>> --
>> What most experimenters take for granted before they begin their
>> experiments is infinitely more interesting than any results to which their
>> experiments lead.
>> -- Norbert Wiener
>>
>> https://www.cse.buffalo.edu/~knepley/
>> <http://www.cse.buffalo.edu/~knepley/>
>>
>

-- 
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener

https://www.cse.buffalo.edu/~knepley/ <http://www.cse.buffalo.edu/~knepley/>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20221130/045d163e/attachment.html>