We employ semilocal density functionals (LSDA, PBE GGA, and meta-GGAs), LSDA+U, a nonlocal range-separated hybrid functional (HSE06), and the random phase approximation (RPA), to assess their performances for the groundstate magnetism and electronic structure of a strongly-correlated metal, rutile VO 2 . Using recent quantum Monte Carlo results as the benchmark, all tested semilocal and hybrid functionals as well as RPA (with PBE inputs) predict the correct magnetic ground states for R-VO 2 . The observed paramagnetism could arise from temperature-disordered local spin moments, or from the thermal destruction of these moments. All semilocal functionals also give the correct ground-state metallicity for R-VO 2 . However, in the ferromagnetic (FM) and anti-ferromagnetic (AFM) phases, LSDA+U and HSE06 incorrectly predict R-VO 2 to be a Mott-Hubbard insulator.