In 2005 we obtained very precise interferometric measurements of the pole-on rapid rotator Vega (A0 V) with the longest baselines of the Center for High Angular Angular Resolution (CHARA) Array and the Fiber Linked Unit for Optical Recombination (FLUOR). For the analysis of these data, we developed a code for mapping sophisticated PHOENIX model atmospheres on to the surface of rotationally distorted stars described by a Roche-von Zeipel formalism. Given a set of input parameters for a star or binary pair, this code predicts the interferometric visibility, spectral energy distribution and high-resolution line spectrum expected for the system. For the gravity-darkened Vega, our model provides a very good match to the K-band interferometric data, a good match to the spectral energy distribution – except below 160 nm – and a rather poor match to weak lines in the high dispersion spectrum where the model appears overly gravity darkened. In 2006, we used the CHARA Array and FLUOR to obtain high precision measurements of the massive, non-eclipsing, double-line spectroscopic binary Spica, a 4-day period system where both components are gravity darkened rapid rotators. These data supplement recent data obtained with the Sydney University Stellar Interferometer (SUSI). Our study follows the classic 1971 study by Herbison-Evans et al. who resolved Spica as a binary with the Narrabri Stellar Intensity Interferometer (NSII). We will report on our progress modeling the new interferometric and archival spectroscopic data, with the goal towards better constraining the apsidal constant.