We calculate the parameters describing elastic I ¼ 1, P-wave ππ scattering using lattice QCD with 2 þ 1 flavors of clover fermions. Our calculation is performed with a pion mass of m π ≈ 320 MeV and a lattice size of L ≈ 3.6 fm. We construct the two-point correlation matrices with both quark-antiquark and twohadron interpolating fields using a combination of smeared forward, sequential and stochastic propagators. The spectra in all relevant irreducible representations for total momenta j ⃗ Pj ≤ ffiffi ffi 3 p 2π L are extracted with two alternative methods: a variational analysis as well as multiexponential matrix fits. We perform an analysis using Lüscher's formalism for the energies below the inelastic thresholds, and investigate several phase shift models, including possible nonresonant contributions. We find that our data are well described by the minimal Breit-Wigner form, with no statistically significant nonresonant component. In determining the ρ resonance mass and coupling we compare two different approaches: fitting the individually extracted phase shifts versus fitting the t-matrix model directly to the energy spectrum. We find that both methods give consistent results, and at a pion mass of am π ¼ 0.18295ð36Þ stat obtain g ρππ ¼ 5.69ð13Þ stat ð16Þ sys , am ρ ¼ 0.4609ð16Þ stat ð14Þ sys , and am ρ =am N ¼ 0.7476ð38Þ stat ð23Þ sys , where the first uncertainty is statistical and the second is the systematic uncertainty due to the choice of fit ranges.