<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?<br><div><br>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">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.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>