<div dir="ltr">Dear all,<br><div><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</div><div><br></div><div>-ksp_type preonly -pc_type lu -pc_factor_mat_solver_package mumps  ?</div><div><br></div><div>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:</div><div><br></div><div>       A_11 A_12<br></div><div>       A_21 A_22</div><div><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</div><div><br></div><div>L_21 U_11 = A_21</div><div>L_11 U_12 = A_12</div><div>L_22 U_22 = A_22 - L_21 U_12</div><div><br></div><div>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.</div><div><br></div><div>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.</div><div><br></div><div>Dr. Toon Weyens<br>ITER Organization</div></div>