def METHOD_NAME(self):
    """Test with NumPyEigensolver"""
    solver = NumPyEigensolver(k=4)
    solver.filter_criterion = self.vibrational_problem.get_default_filter_criterion()
    esc = ExcitedStatesEigensolver(self.qubit_mapper, solver)
    results = esc.solve(self.vibrational_problem)
    self._assert_energies(results.computed_vibrational_energies, self.reference_energies)