[petsc-users] petsc and mumps for compute user-specified set of entries in inv(A)

[Hong]

Marius :
I added MatMumpsGetInverse(), see
https://bitbucket.org/petsc/petsc/commits/f5e16c35adb7810c6c977ea1c9cd95e4afb8512b
It works in sequential and parallel.

It is in the branch hzhang/mumps-invA, and will be merged to petsc-master
after it passed our regression tests.

Let me know if you see any problems.
Hong

On Sat, Jan 20, 2018 at 11:42 AM, Hong <hzhang at mcs.anl.gov> wrote:

> Marius :
> Current petsc-mumps interface supports ICNTL(30), e.g., runtime option '-mat_mumps_icntl_30
> 0' as default.
> However, no user ever has tested it. I tested it
> using petsc/src/mat/examples/tests/ex130.c with
>  '-mat_mumps_icntl_30 1' and got an error in MatSolve() -- additional
> work needs to be down here.
>
> I'll read mumps user manual and investigate it next week.
>
> Meanwhile, you may give it a try and figure out how to set other
> parameters, such as icntl_27 and allocate correct rhs and solution
> vectors.
>
> Hong
>
>> Hi,
>>
>> Is it possible to interface MUMPS to  compute user-specified set of
>> entries in inv(A) (ICNTL(30)) using petsc ?
>>
>> best,
>> marius
>>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20180207/9ee2cdff/attachment.html>