A Fast Phase Space Method for Computing Creeping Rays 1