Continuation methods utilizing boundary value solvers are an effective tool for computing unstable periodic orbits of dynamical systems. The standard implementation for these procedures is AUTO, which often requires very fine meshes on problems with multiple time scales. We will explore alternate strategies based on geometric singular perturbation theory for computing periodic orbits in these systems.