Just subjective, but I’ve achieved decent parallel performance as measured by strong scaling with MUMPS going up to about 10M DoF for 2D problems. Specifically this was Stokes coupled with the heat equation.
I’ve not been able to demonstrate good strong scaling for 3D problems. Perhaps others have suggestions.
See also How to choose the optimal solver for a PDE problem? - #2 by nate