Optimizations that produced a linear transition state were retained so that the barrier to linearity could be computed. For all structures that were minima, reaction energies were computed, as well as shifts in C-H and C-N bond distances and corresponding harmonic stretching frequencies compared to HCN. All functionals were assigned a percentile to assess performance for reaction energies, bond distances, shifts in bond distances, harmonic vibrational frequencies, and shifts in harmonic vibrational frequencies. For some of the Psi4 data, confirmatory results were also obtained with the MOLPRO or QCHEM software packages. To assess the trends in the DFT data, the functionals were grouped according to the following categories, as adopted from Morgante and Peverati: Group 0 (LDA, HF, LC, SE) -local spin density approximations, Hartree Fock, low-cost methods, and semiempirical methods; Group 1 (GGA) -generalized gradient approximations and nonseparable gradient approximations; Group 2 (mGGA) -meta-GGA and meta-NGA functionals; Group 3 (GH-GGA) -global hybrid GGA and NGA functionals; Group 4 (GH-mGGA) -global hybrid meta-GGA and meta-NGA functionals; Group 5 (RSH) -long range corrected hybrid functionals; Group 6 (DH) -double-hybrid functionals. These groupings are in general accord with commonly assigned rungs of the DFT ladder but do not presuppose a strict hierarchy.