Computational modeling of neoclassical and resistive magnetohydrodynamic tearing modes in tokamaks