Computational modeling of magnetoconvection: Effects of discretization method, grid refinement and grid stretching