<div dir="ltr">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><div><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?</div></div><br><div class="gmail_quote"><div dir="ltr">On Sun, Aug 27, 2017 at 12:54 AM Barry Smith <<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
  Toon,<br>
<br>
   We don't have any such way of doing this. Note that sparse factorization packages store the "factors" in specialized data structures that are unlikely to be acessable to perform such operations as you describe below.<br>
<br>
   I think what you want to do is an over optimization not worth the coding effort.<br>
<br>
  Barry<br>
<br>
> On Aug 25, 2017, at 5:01 AM, Toon Weyens <<a href="mailto:weyenst@gmail.com" target="_blank">weyenst@gmail.com</a>> wrote:<br>
><br>
> Dear all,<br>
><br>
> For a Bounday Element Method problem I require the solution of a system of linear equations with multiple right-hand sides. Though this is a dense system, I still want to do it via Petsc. Would the best way to do this be through something such as<br>
><br>
> -ksp_type preonly -pc_type lu -pc_factor_mat_solver_package mumps  ?<br>
><br>
> Furthermore, I refine my grid in a regular way, which leads to a new system of equations of size 2N x 2N where N is the original number of unknowns before refining it. The upper-left matrix A_11 in this new system is identical to the unrefined matrix (multiplied by 0.5). The other three blocks A_12, A_21 and A_22 are new:<br>
><br>
>        A_11 A_12<br>
>        A_21 A_22<br>
><br>
> The question is now whether knowledge about the unrefined matrix A_11 can be used to speed up calculation of the refined system?<br>
><br>
> I was thinking, for example, about using the LU decomposition of A_11 to calculate the LU decomposition of the entire matrix A using the well-known formula's<br>
><br>
> L_21 U_11 = A_21<br>
> L_11 U_12 = A_12<br>
> L_22 U_22 = A_22 - L_21 U_12<br>
><br>
> where L_21 and U_12 are full matrices, L_11, L_22 are lower triangular matrices and U_11, U_22 are upper triangular.<br>
><br>
> Is there any way to do this in Petsc?<br>
><br>
> Or is there a better thing I can do?<br>
><br>
> Thanks very much in advance!<br>
><br>
> Regards.<br>
><br>
> Dr. Toon Weyens<br>
> ITER Organization<br>
<br>
</blockquote></div>