A small number of active galactic nuclei are known to exhibit prominent double peak emission profiles that are well-fitted by a relativistic accretion disk model. We develop a Monte Carlo code to compute the linear polarization of a double peaked broad emission line arising from Thomson scattering. A Keplerian accretion disk is adopted for the double peak emission line region and the geometry is assumed to be Schwarzschild. Far from the accretion disk where flat Minkowski geometry is appropriate, we place an azimuthally symmetric scattering region in the shape of a spherical shell sliced with ${Delta}{mu}=0.1$. Adopting a Monte Carlo method we generate line photons in the accretion disk in arbitrary directions in the local rest frame and follow the geodesic paths of the photons until they hit the scattering region. The profile of the polarized flux is mainly determined by the relative location of the scattering region with respect to the emission source. When the scattering region is in the polar direction, the degree of linear polarization also shows a double peak structure. Under favorable conditions we show that up to 0.6% linear polarization may be obtained. We conclude that spectropolarimetry can be a powerful probe to reveal much information regarding the accretion disk geometry of these active galactic nuclei.