We show how to apply the absorbing Markov chain Monte Carlo algorithm of Novotny to simulate kinetically constrained models of glasses. We consider in detail one-spin facilitated models, such as the East model and its generalizations to arbitrary dimensions. We investigate how to maximise the efficiency of the algorithms, and show that simulation times can be improved on standard continuous time Monte Carlo by several orders of magnitude. We illustrate the method with equilibrium and aging results. These include a study of relaxation times in the East model for dimensions d=1 to d=13, which provides further evidence that the hierarchical relaxation in this model is present in all dimensions. We discuss how the method can be applied to other kinetically constrained models. ; Comment: 8 pages, 4 figures