<div dir="ltr">Hi Hong,<div><br></div><div>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.</div><div><br></div><div>Best,</div><div>Evan</div><div><br></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Oct 11, 2017 at 9:26 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:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><br><div class="gmail_extra"><br><div class="gmail_quote"><span class="">On Wed, Oct 11, 2017 at 11:14 AM, Evan Um <span dir="ltr"><<a href="mailto:evanum@gmail.com" target="_blank">evanum@gmail.com</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">Hi Hong,<div><br></div><div>Thanks for your kind email. I write another email to make sure I understand it correctly. Does the zipped file from <a href="https://www.mcs.anl.gov/petsc/developers/index.html#browsing" target="_blank">https://www.mcs.anl.gov/p<wbr>etsc/developers/index.html#bro<wbr>wsing</a> have the updated feature?</div><div><br></div><div>Best,</div><div>Evan</div><div><br></div><ul style="color:rgb(0,0,0);font-family:"Times New Roman";font-size:medium;background-color:rgb(213,234,255)"><li><a href="https://bitbucket.org/petsc/petsc/get/master.tar.gz" target="_blank">https://bitbucket.org/petsc/pe<wbr>tsc/get/master.tar.gz</a></li></ul></div></blockquote></span><div>I do not see it :-(</div><div>Get petsc-dev by</div><ul style="color:rgb(0,0,0);font-family:Times;font-size:medium;background-color:rgb(213,234,255)"><li><code>git clone <a href="https://bitbucket.org/petsc/petsc" target="_blank">https://bitbucket.org/petsc/<wbr>petsc</a></code></li></ul><div>Then do following:</div><div>git checkout master</div><div>git pull</div><div><br></div><div>Check file </div><div>src/mat/impls/aij/mpi/mumps/<wbr>mumps.c to see if ICNTL(35) is there or not.</div><span class="HOEnZb"><font color="#888888"><div><br></div><div>Hong</div></font></span><div><div class="h5"><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><font color="#000000" face="Times New Roman" size="3"><br></font></div><div> </div></div><div class="m_-3202138047207851770gmail-HOEnZb"><div class="m_-3202138047207851770gmail-h5"><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Oct 11, 2017 at 7:56 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">Evan:</div><div class="gmail_quote">The branch is already merged to master. Just use petsc master branch.</div><span class="m_-3202138047207851770gmail-m_-4903602127405685030HOEnZb"><font color="#888888"><div class="gmail_quote">Hong</div></font></span><div><div class="m_-3202138047207851770gmail-m_-4903602127405685030h5"><div class="gmail_quote"><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">Hi Matt,<div><br></div><div>Still unclear to me. I go to <a href="https://bitbucket.org/petsc/petsc/addon/pipelines/home#!/" target="_blank">https://bitbucket.org/petsc<wbr>/petsc/addon/pipelines/home#!/</a> and then <a href="https://bitbucket.org/petsc/petsc/branches/" target="_blank">https://bitbucket.org/pet<wbr>sc/petsc/branches/</a>. I don't see any tarred file or directory named "<span style="color:rgb(80,0,80);font-size:12.8px">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.</span></div><span class="m_-3202138047207851770gmail-m_-4903602127405685030m_6638224014225267443HOEnZb"><font color="#888888"><div><span style="color:rgb(80,0,80);font-size:12.8px"><br></span></div><div><span style="color:rgb(80,0,80);font-size:12.8px">Evan</span></div><div><span style="color:rgb(80,0,80);font-size:12.8px"><br></span></div><div><span style="color:rgb(80,0,80);font-size:12.8px"><br></span></div><div><span style="color:rgb(80,0,80);font-size:12.8px"><br></span></div></font></span></div><div class="m_-3202138047207851770gmail-m_-4903602127405685030m_6638224014225267443HOEnZb"><div class="m_-3202138047207851770gmail-m_-4903602127405685030m_6638224014225267443h5"><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Oct 11, 2017 at 1:54 AM, Matthew Knepley <span dir="ltr"><<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</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"><span>On Tue, Oct 10, 2017 at 10:30 PM, Evan Um <span dir="ltr"><<a href="mailto:evanum@gmail.com" target="_blank">evanum@gmail.com</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">Dear Hong,<div><br></div><div>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.</div></div></blockquote><div><br></div></span><div>Hong means a branch. You get the dev repository and checkout that branch</div><div><br></div><div> git checkout hzang/update-mumps-5.1.1-cntl</div><div><br></div><div> Thanks,</div><div><br></div><div> Matt</div><div><div class="m_-3202138047207851770gmail-m_-4903602127405685030m_6638224014225267443m_785853025285377277h5"><div> </div><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> . 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.</div><div><br></div><div>Best,</div><div>Evan</div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Oct 3, 2017 at 8:34 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">Evan,<div>ICNTL(35) and CNTL(7) are added to petsc-mumps interface in branch </div><div>hzhang/update-mumps-5.1.1-cntl</div><div><br></div><div>You may give it a try. Once it passes our regression tests, I'll merge it to petsc master branch.</div><span class="m_-3202138047207851770gmail-m_-4903602127405685030m_6638224014225267443m_785853025285377277m_6394383616296017593gmail-m_-844966978622877325HOEnZb"><font color="#888888"><div><br></div></font></span><div><span class="m_-3202138047207851770gmail-m_-4903602127405685030m_6638224014225267443m_785853025285377277m_6394383616296017593gmail-m_-844966978622877325HOEnZb"><font color="#888888">Hong</font></span><div><div class="m_-3202138047207851770gmail-m_-4903602127405685030m_6638224014225267443m_785853025285377277m_6394383616296017593gmail-m_-844966978622877325h5"><br><div class="gmail_extra"><br><div class="gmail_quote">On Sun, Sep 24, 2017 at 8:08 PM, 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">I'll check it.<span class="m_-3202138047207851770gmail-m_-4903602127405685030m_6638224014225267443m_785853025285377277m_6394383616296017593gmail-m_-844966978622877325m_-559390365903874404gmail-HOEnZb"><font color="#888888"><div>Hong</div></font></span></div><div class="m_-3202138047207851770gmail-m_-4903602127405685030m_6638224014225267443m_785853025285377277m_6394383616296017593gmail-m_-844966978622877325m_-559390365903874404gmail-HOEnZb"><div class="m_-3202138047207851770gmail-m_-4903602127405685030m_6638224014225267443m_785853025285377277m_6394383616296017593gmail-m_-844966978622877325m_-559390365903874404gmail-h5"><div class="gmail_extra"><br><div class="gmail_quote">On Sun, Sep 24, 2017 at 3:42 PM, Evan Um <span dir="ltr"><<a href="mailto:evanum@gmail.com" target="_blank">evanum@gmail.com</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">Hi Barry,<div><br></div><div>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.</div><div><br></div><div>Best,</div><div>Evan</div><div> </div></div><div class="m_-3202138047207851770gmail-m_-4903602127405685030m_6638224014225267443m_785853025285377277m_6394383616296017593gmail-m_-844966978622877325m_-559390365903874404gmail-m_8252441951155370891HOEnZb"><div class="m_-3202138047207851770gmail-m_-4903602127405685030m_6638224014225267443m_785853025285377277m_6394383616296017593gmail-m_-844966978622877325m_-559390365903874404gmail-m_8252441951155370891h5"><div class="gmail_extra"><br><div class="gmail_quote">On Sat, Sep 23, 2017 at 6:45 PM, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@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"><span><br>
> On Sep 23, 2017, at 8:38 PM, Evan Um <<a href="mailto:evanum@gmail.com" target="_blank">evanum@gmail.com</a>> wrote:<br>
><br>
> Dear PETSC Users,<br>
><br>
> 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?<br>
<br>
</span> You don't pass factored matrices you just pass the original matrix and use -pc_type lu -pc_factor_mat_solver_package mumps<br>
<span><br>
> Can PETSC call the latest MUMPS that supports block low rank approximation?<br>
<br>
</span> No, send us info on it and we'll see if we can add an interface<br>
<div class="m_-3202138047207851770gmail-m_-4903602127405685030m_6638224014225267443m_785853025285377277m_6394383616296017593gmail-m_-844966978622877325m_-559390365903874404gmail-m_8252441951155370891m_-5478992268712916641HOEnZb"><div class="m_-3202138047207851770gmail-m_-4903602127405685030m_6638224014225267443m_785853025285377277m_6394383616296017593gmail-m_-844966978622877325m_-559390365903874404gmail-m_8252441951155370891m_-5478992268712916641h5"><br>
<br>
><br>
> In advance, thank you very much for your comments.<br>
><br>
> Best,<br>
> Evan<br>
><br>
><br>
><br>
><br>
><br>
<br>
</div></div></blockquote></div><br></div>
</div></div></blockquote></div><br></div>
</div></div></blockquote></div><br></div></div></div></div></div>
</blockquote></div><br></div>
</blockquote></div></div></div><span class="m_-3202138047207851770gmail-m_-4903602127405685030m_6638224014225267443m_785853025285377277HOEnZb"><font color="#888888"><br><br clear="all"><div><br></div>-- <br><div class="m_-3202138047207851770gmail-m_-4903602127405685030m_6638224014225267443m_785853025285377277m_6394383616296017593gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div><div><br></div><div><a href="http://www.caam.rice.edu/~mk51/" target="_blank">https://www.cse.buffalo.edu/~k<wbr>nepley/</a><br></div></div></div></div></div>
</font></span></div></div>
</blockquote></div><br></div>
</div></div></blockquote></div><br></div></div></div></div>
</blockquote></div><br></div>
</div></div></blockquote></div></div></div><br></div></div>
</blockquote></div><br></div>