We report on our progress in computing the excitation spectrum in Lattice QCD. We focus on the isospin 0, 1 and 2 channels using the stochastic LapH algorithm for the quark propagators. For the isospin-0 channel, a new glueball operator constructed from the stochastic LapH operator is included in the variational basis along with the isoscalar meson and pi-pi operators. A representative signal for each channel is presented and the feasibility of extending the calculations to larger lattices is discussed.