Fast Solvers And Uncertainty Quantification For Models Of Magnetohydrodynamics