Edoardo alinovi <edoardo.alinovi at gmail.com> writes: > Ps: I am solving Navier-Stokes equation (pressure-correction equation in > particular). The matrix is sparse and in mpiaij fotmat. If you're just solving a pressure projection then you can and should use multigrid. For example, -pc_type gamg.