An adaptive scalable fully implicit algorithm based on stabilized finite element for reduced visco-resistive MHD.