Phase screens are widely used to simulate the effects of atmospheric turbulence on the phase fluctuations of a light beam. We compare three sampled-based phase-screen-generation methods (the fast-Fourier-transform, subharmonic, and covariance-matrix methods), in terms of accuracy and simulation time. We show that the covariance method generates the phase screens most accurately, with simulation time comparable to the other sampled-based methods.