The performance of 39 different LDA, GGA, meta-GGA, and hybrid density functionals has been evaluated, for calculating forward and reverse barrier heights of 10 gas-phase reactions involving hydrogen. The reactions are all relevant to astrochemistry. Special focus is put on the applicability of DFT for calculating the rates of corresponding surface hydrogenation reactions that are relevant to the chemistry of ice-coated interstellar grains. General trends in the performance of the density functionals for reactions involving H atoms, H-2, and OH are discussed. The OH+CO reaction is shown to be a very problematic case for DFT. The best overall performance is found for the hybrid density functionals, such as MPW1K, B97-1, B97-2, and B1B95. For several reactions, the HCTH GGA functionals and the VS98 and OLAP3 meta-GGA functionals also give results that are almost as good as those of the hybrid functionals.