Mumps should take a sparse vector for the RHS. https://github.com/JuliaSparse/MUMPS.jl/blob/master/src/MUMPS_4.10.0/SCILAB/examples/sparseRHS_example.sce