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 <hzh...@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 >> > >
