While I have no issue using mpirun -np N to control the number of processes when using a Krylov solver, I’m finding that if I try to use this with either superlu_dist or mumps, it just runs with all the cores on my system, no matter what I set for N. How do I control this?