I review some computational methods for calculating vibrational spectra. They all use iterative eigensolvers to compute eigenvalues of a Hamiltonian matrix by evaluating matrix-vector products (MVPs). A direct-product basis can be used for molecules with five or fewer atoms. This is done by exploiting the structure of the basis and the structure of a direct product quadrature grid. I outline three methods that can be used for molecules with more than five atoms. The first uses contracted basis functions and an intermediate (F) matrix. The second uses Smolyak quadrature and a pruned basis. The third uses a tensor rank reduction scheme.