The ShakeNBreak method was used to identify the lowest energy structures for point defects in CuI. Upon generation of a defect, a series of targeted bond distortions (the quantity and direction of which are determined using simple valence chemistry rules) is applied, followed by a small, random "rattling" of the structure, in order to break the local symmetry. These distorted structures, along with the "unperturbed" starting structure, are then relaxed using coarse calculation parameters to get a relatively cheap sampling of the complex defect potential energy surface. The lowest energy structure identified from this initial search is then relaxed using full accuracy parameters. This method allows the search for a defect ground state to begin from several, distinct points on the potential energy surface, decreasing the probability of getting stuck in a high energy local minimum.