Multi-phase computations in geometrical optics