The most commonly encountered form of echo planar spectroscopy involves oscillating gradients in one spatial dimension during readout. Data are consequently not sampled on a Cartesian grid. A fast gridding algorithm applicable to this particular situation is presented. The method is optimal, i.e., it performs as well as the full discrete Fourier transform for band limited signals while allowing for use of the fast Fourier transform. The method is demonstrated for reconstruction of data that are partially undersampled in the time domain. The advantages of undersampling are lower hardware requirements or fewer interleaves per acquisition. The method is of particular interest when large bandwidths are needed (e.g., for high field scanning) and for scanners with limited gradient performance. The unavoidable artifacts resulting from undersampling are demonstrated to be acceptable for spectroscopy with long echo times.