<div dir="ltr"><div class="gmail_quote"><div dir="ltr"><div>Dear Petsc team,</div><div><br></div><div>I am (still) trying to solve Biot's poroelasticity problem :</div><div> <img width="237" height="54" class="m_461426253857996911gmail-CToWUd" style="margin-right: 0px;" alt="Images intégrées 1" src="cid:ii_158f8fd5ddeb7876"></div><div><br></div><div>I am using a mixed P2-P1 finite element discretization. The matrix of the discretized system in binary format is attached to this email.</div><div><br></div><div>I am using the fieldsplit framework to solve the linear system. Since I am facing some troubles, I have decided to go back to simple things. Here are the options I am using : </div><div><br></div><div>-ksp_rtol 1.0e-5<br>-ksp_type fgmres<br>-pc_type fieldsplit<br>-pc_fieldsplit_schur_<wbr>factorization_type full<br>-pc_fieldsplit_type schur<br>-pc_fieldsplit_schur_<wbr>precondition selfp<br>-fieldsplit_0_pc_type lu<br>-fieldsplit_0_pc_factor_mat_<wbr>solver_package mumps<br>-fieldsplit_0_ksp_type preonly<br>-fieldsplit_0_ksp_converged_<wbr>reason<br>-fieldsplit_1_pc_type lu<br>-fieldsplit_1_pc_factor_mat_<wbr>solver_package mumps<br>-fieldsplit_1_ksp_type preonly<br>-fieldsplit_1_ksp_converged_<wbr>reason<br><br></div><div>On a single proc, everything runs fine : the solver converges in 3 iterations, according to the theory (see Run-1-proc.txt [contains -log_view]).</div><div><br></div><div>On 2 procs, the solver converges in 28 iterations (see Run-2-proc.txt).</div><div><br></div><div>On 3 procs, the solver converges in 91 iterations (see Run-3-proc.txt).</div><div><br></div><div>I do not understand this behavior : since MUMPS is a parallel direct solver, shouldn't the solver converge in max 3 iterations whatever the number of procs?</div><div><br></div><div><br></div><div>Thanks for your precious help,</div><div>Nicolas</div></div>
</div><br></div>