Very interesting! I might try this afterwards, with H2lib. 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?
Thanks! On Mon, Aug 28, 2017 at 12:15 PM Jose E. Roman <[email protected]> wrote: > > > El 28 ago 2017, a las 10:15, Toon Weyens <[email protected]> escribió: > > > > 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. > > > > In the long term it would in any case be useful to have a solver that > uses H-matrices, such as H2lib. > > > > Is there any chance Petsc is thinking about moving to this kind of > matrices as well in the future? > > > > > > You can use H2lib in a PETSc program by wrapping it in a shell matrix. We > did this some time ago https://doi.org/10.1016/j.cam.2012.07.021 > > Also, PETSc has an interface to STRUMPACK, but only to the sparse part and > I guess what you need is the dense part. > > Jose > >
