Current State of MHD Simulation of Massive Star Winds In magnetic massive stars, the competition between highly non-linear radiative forces and magnetic field leads to a highly complex interaction. Such an interplay cannot usually be modeled analytically, instead numerical modeling become a necessary tool. In this talk, I will discuss how numerical magnetohydrodynamic (MHD) simulations are employed to understand the nature of such magnetized massive star winds. I will also briefly outline the challenges involved in such simulations, especially in 3D.