<div dir="ltr">Hi Sherry,<div><br></div><div>Did you guys have a standard PETSc example that you run to test the interface? Should I be asking someone else?</div><div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Aug 28, 2017 at 10:24 AM, Toon Weyens <span dir="ltr"><<a href="mailto:weyenst@gmail.com" target="_blank">weyenst@gmail.com</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">Hey, as indicated here: <a href="https://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/Mat/MatSTRUMPACKSetHSSMinSize.html" target="_blank">https://www.mcs.anl.gov/<wbr>petsc/petsc-dev/docs/<wbr>manualpages/Mat/<wbr>MatSTRUMPACKSetHSSMinSize.html</a><div><br></div><div>there is an example in <a href="https://www.mcs.anl.gov/petsc/petsc-dev/src/ksp/ksp/examples/tutorials/ex52.c.html" target="_blank">src/ksp/ksp/examples/<wbr>tutorials/ex52.c.html</a></div><div><br></div><div>However, in my 3.7.6 distribution ex52.c does not have lines 29-31:</div><div><pre width="80"><a> 29: </a><font color="#A020F0">#if defined(PETSC_HAVE_STRUMPACK)</font>
<a> 30: </a>  <a href="https://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/Sys/PetscBool.html#PetscBool" target="_blank">PetscBool</a>      flg_strumpack=<a href="https://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/Sys/PETSC_FALSE.html#PETSC_FALSE" target="_blank">PETSC_FALSE</a>;
<a> 31: </a><font color="#A020F0">#endif</font></pre></div></div><br><div class="gmail_quote"><div dir="ltr">On Mon, Aug 28, 2017 at 4:01 PM Matthew Knepley <<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Mon, Aug 28, 2017 at 9:17 AM, Toon Weyens <span dir="ltr"><<a href="mailto:weyenst@gmail.com" target="_blank">weyenst@gmail.com</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">Very interesting! I might try this afterwards, with H2lib.<br><br>Also, I must say that the strumpack interface looks promising. What did you mean by "only for the sparse part", by the way? Also, any idea when it is going to be released in an official petsc release, preferably with some documentation?</div></blockquote><div><br></div></div></div></div><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><div>I thought STRUMPACK was in 3.7.6</div><div><br></div><div>  Matt</div></div></div></div><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>Thanks!</div></div><br><div class="gmail_quote"><div dir="ltr">On Mon, Aug 28, 2017 at 12:15 PM Jose E. Roman <<a href="mailto:jroman@dsic.upv.es" target="_blank">jroman@dsic.upv.es</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
> El 28 ago 2017, a las 10:15, Toon Weyens <<a href="mailto:weyenst@gmail.com" target="_blank">weyenst@gmail.com</a>> escribió:<br>
><br>
> Thank you Barry, that explains why I couldn't find information about it. I am now going to implement this straight-forward implementation as a first step.<br>
><br>
> In the long term it would in any case be useful to have a solver that uses H-matrices, such as H2lib.<br>
><br>
> Is there any chance Petsc is thinking about moving to this kind of matrices as well in the future?<br>
><br>
><br>
<br>
You can use H2lib in a PETSc program by wrapping it in a shell matrix. We did this some time ago <a href="https://doi.org/10.1016/j.cam.2012.07.021" rel="noreferrer" target="_blank">https://doi.org/10.1016/j.cam.<wbr>2012.07.021</a><br>
<br>
Also, PETSc has an interface to STRUMPACK, but only to the sparse part and I guess what you need is the dense part.<br>
<br>
Jose<br>
<br>
</blockquote></div>
</blockquote></div></div></div><div dir="ltr"><div class="gmail_extra"><br><br clear="all"><span class="HOEnZb"><font color="#888888"><div><br></div>-- <br><div class="m_-551128285681966165m_6178142451908864997gmail_signature" data-smartmail="gmail_signature"><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">http://www.caam.rice.edu/~<wbr>mk51/</a><br></div></div></div>
</font></span></div></div></blockquote></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><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">http://www.caam.rice.edu/~mk51/</a><br></div></div></div>
</div>