<div dir="ltr"><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline">Marius :</span><div>I added MatMumpsGetInverse(), see</div><div><a href="https://bitbucket.org/petsc/petsc/commits/f5e16c35adb7810c6c977ea1c9cd95e4afb8512b">https://bitbucket.org/petsc/petsc/commits/f5e16c35adb7810c6c977ea1c9cd95e4afb8512b</a></div><div>It works in sequential and parallel.</div><div><br></div><div>It is in the branch hzhang/mumps-invA, and will be merged to petsc-master after it passed our regression tests.</div><div><br></div><div>Let me know if you see any problems.</div><div>Hong<br><div class="gmail_extra"><br><div class="gmail_quote">On Sat, Jan 20, 2018 at 11:42 AM, Hong <span dir="ltr"><<a href="mailto:hzhang@mcs.anl.gov" target="_blank">hzhang@mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">Marius :</div><div class="gmail_quote">Current petsc-mumps interface supports <span style="font-family:Verdana;font-size:12px">ICNTL(30), e.g., runtime option '</span><font face="Verdana"><span style="font-size:12px">-mat_mumps_icntl_30 0' as default.</span></font></div><div class="gmail_quote"><font face="Verdana"><span style="font-size:12px">However, no user ever has tested it. I tested it using petsc/src/mat/examples/<wbr>tests/ex130.c with</span></font></div><div class="gmail_quote"><span style="font-family:Verdana;font-size:12px"> '</span><font face="Verdana"><span style="font-size:12px">-mat_mumps_icntl_30 1' and got an error in MatSolve() -- additional work needs to be down here.</span></font></div><div class="gmail_quote"><font face="Verdana"><span style="font-size:12px"><br></span></font></div><div class="gmail_quote"><font face="Verdana"><span style="font-size:12px">I'll read mumps user manual and investigate it next week.</span></font><font face="Verdana"><span style="font-size:12px"><br></span></font></div><div class="gmail_quote"><font face="Verdana"><span style="font-size:12px"><br></span></font></div><div class="gmail_quote"><font face="Verdana"><span style="font-size:12px">Meanwhile, you may give it a try and figure out how to set other parameters, such as </span></font><font face="Verdana"><span style="font-size:12px">icntl_27 and allocate correct rhs and solution vectors.</span></font></div><span class="gmail-HOEnZb"><font color="#888888"><div class="gmail_quote"><font face="Verdana"><span style="font-size:12px"><br></span></font></div></font></span><div class="gmail_quote"><span class="gmail-HOEnZb"><font color="#888888"><font face="Verdana"><span style="font-size:12px">Hong<br></span></font></font></span><div><div class="gmail-h5"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div><div style="font-family:Verdana;font-size:12px"><div>Hi,</div>

<div> </div>

<div>Is it possible to interface MUMPS to  compute user-specified set of entries in inv(A) (ICNTL(30)) using petsc ?</div>

<div> </div>

<div>best,</div>

<div>marius</div></div></div>
</blockquote></div></div></div><br></div></div>
</blockquote></div><br></div></div></div>