Stability and temporal error estimate of scalar auxiliary variable schemes for the magnetohydrodynamics equations with variable density