Calculation of highly excited vibrational levels: a prediagonalized Davidson scheme