[petsc-users] Using factored complex matrices from MUMPS as a preconditioner in PETSC

[Evan Um]

Hi Hong,

I was able to check out the dev version. Now, I see ICNTL(35) inside
mumps.c  I am going to test the functionality. If I see any issue, I will
report it to you. Thank you very much for your kind help.

Best,
Evan



On Wed, Oct 11, 2017 at 9:26 AM, Hong <hzhang at mcs.anl.gov> wrote:

>
>
> On Wed, Oct 11, 2017 at 11:14 AM, Evan Um <evanum at gmail.com> wrote:
>
>> Hi Hong,
>>
>> Thanks for your kind email. I write another email to make sure I
>> understand it correctly. Does the zipped file from
>> https://www.mcs.anl.gov/petsc/developers/index.html#browsing have the
>> updated feature?
>>
>> Best,
>> Evan
>>
>>
>>    - https://bitbucket.org/petsc/petsc/get/master.tar.gz
>>
>> I do not see it :-(
> Get petsc-dev by
>
>    - git clone https://bitbucket.org/petsc/petsc
>
> Then do following:
> git checkout master
> git pull
>
> Check file
> src/mat/impls/aij/mpi/mumps/mumps.c to see if  ICNTL(35) is there or not.
>
> Hong
>
>>
>>
>>
>> On Wed, Oct 11, 2017 at 7:56 AM, Hong <hzhang at mcs.anl.gov> wrote:
>>
>>> Evan:
>>> The branch is already merged to master. Just use petsc master branch.
>>> Hong
>>>
>>> Hi Matt,
>>>>
>>>> Still unclear to me. I go to https://bitbucket.org/petsc
>>>> /petsc/addon/pipelines/home#!/ and then https://bitbucket.org/pet
>>>> sc/petsc/branches/. I don't see any tarred file or directory named "update-mumps-5.1.1-cntl".
>>>> Could you explain how to download the modified version from the site a
>>>> little bit in detail? Thank you very much for your help.
>>>>
>>>> Evan
>>>>
>>>>
>>>>
>>>>
>>>> On Wed, Oct 11, 2017 at 1:54 AM, Matthew Knepley <knepley at gmail.com>
>>>> wrote:
>>>>
>>>>> On Tue, Oct 10, 2017 at 10:30 PM, Evan Um <evanum at gmail.com> wrote:
>>>>>
>>>>>> Dear Hong,
>>>>>>
>>>>>> I just tried to check PETSC develeoper website but couldn't find the
>>>>>> updated dev version with name of hzang/update-mumps-5.1.1-cntl.
>>>>>>
>>>>>
>>>>> Hong means a branch. You get the dev repository and checkout that
>>>>> branch
>>>>>
>>>>>   git checkout  hzang/update-mumps-5.1.1-cntl
>>>>>
>>>>>   Thanks,
>>>>>
>>>>>    Matt
>>>>>
>>>>>
>>>>>> . Could you please let me know the location of the updated dev
>>>>>> version? Where do I need to visit to check out the dev version with the new
>>>>>> control switches? Thank you very much for your help.
>>>>>>
>>>>>> Best,
>>>>>> Evan
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>> On Tue, Oct 3, 2017 at 8:34 AM, Hong <hzhang at mcs.anl.gov> wrote:
>>>>>>
>>>>>>> Evan,
>>>>>>> ICNTL(35) and CNTL(7) are added to petsc-mumps interface in branch
>>>>>>> hzhang/update-mumps-5.1.1-cntl
>>>>>>>
>>>>>>> You may give it a try. Once it passes our regression tests, I'll
>>>>>>> merge it to petsc master branch.
>>>>>>>
>>>>>>> Hong
>>>>>>>
>>>>>>>
>>>>>>> On Sun, Sep 24, 2017 at 8:08 PM, Hong <hzhang at mcs.anl.gov> wrote:
>>>>>>>
>>>>>>>> I'll check it.
>>>>>>>> Hong
>>>>>>>>
>>>>>>>> On Sun, Sep 24, 2017 at 3:42 PM, Evan Um <evanum at gmail.com> wrote:
>>>>>>>>
>>>>>>>>> Hi Barry,
>>>>>>>>>
>>>>>>>>> Thanks for your comments. To activate block low rank (BLR)
>>>>>>>>> approximation in MUMPS version 5.1.1, a user needs to turn on the
>>>>>>>>> functionality (i.e. ICNTL(35)=1) and specify the tolerance value (e.g.
>>>>>>>>> CNTL(7)=1e-4). In PETSC, I think that we can set up ICNTL and CNTL
>>>>>>>>> parameters for MUMPS. I was wondering if we can still use BLR approximation
>>>>>>>>> for a preconditioner for Krylov solvers.
>>>>>>>>>
>>>>>>>>> Best,
>>>>>>>>> Evan
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> On Sat, Sep 23, 2017 at 6:45 PM, Barry Smith <bsmith at mcs.anl.gov>
>>>>>>>>> wrote:
>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>> > On Sep 23, 2017, at 8:38 PM, Evan Um <evanum at gmail.com> wrote:
>>>>>>>>>> >
>>>>>>>>>> > Dear PETSC Users,
>>>>>>>>>> >
>>>>>>>>>> > My system matrix comes from finite element modeling and is
>>>>>>>>>> complex and unstructured. Its typical size is a few millions-by a few
>>>>>>>>>> millions. I wondering if I can use MUMPS parallel direct solver as a
>>>>>>>>>> preconditioner in PETSC. For example, I want to pass factored matrices to
>>>>>>>>>> Krylov iterative solvers such as QMR. Is there any PETSC+MUMPS example code
>>>>>>>>>> for the purpose?
>>>>>>>>>>
>>>>>>>>>>   You don't pass factored matrices you just pass the original
>>>>>>>>>> matrix and use -pc_type lu -pc_factor_mat_solver_package mumps
>>>>>>>>>>
>>>>>>>>>> > Can PETSC call the latest MUMPS that supports block low rank
>>>>>>>>>> approximation?
>>>>>>>>>>
>>>>>>>>>>   No, send us info on it and we'll see if we can add an interface
>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>> >
>>>>>>>>>> > In advance, thank you very much for your comments.
>>>>>>>>>> >
>>>>>>>>>> > Best,
>>>>>>>>>> > Evan
>>>>>>>>>> >
>>>>>>>>>> >
>>>>>>>>>> >
>>>>>>>>>> >
>>>>>>>>>> >
>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>
>>>>>>>>
>>>>>>>
>>>>>>
>>>>>
>>>>>
>>>>> --
>>>>> 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.caam.rice.edu/~mk51/>
>>>>>
>>>>
>>>>
>>>
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20171011/21fe3280/attachment-0001.html>