Two level scheme solvers for nuclear spectroscopy