\vB(r)=\frac{3}{4\pi r^3}\int_{x<r}  \vv({\vx})  \dd^3 x  \; .$$ where $\vv(\vx)$ is the 3D peculiar velocity field as a function of the comoving coordinate $\vx$. Beneath this innocuous expression lie a multitude of nuisances. An unbiased estimate of $\vB$ requires knowledge of $\vv$ sampled uniformly overall the volume. However, observational probes of peculiar velocities measurements are available only for a few thousand galaxies with a patchy coverage of the local Universe. Further, peculiar velocity probes such as the TF relation allow us to constrain only the radial component of the peculiar velocities of galaxies.

Recently compiled data on peculiar velocities has triggered renewed interest in the analysis of large scale flows, including the bulk flow [@DN10; @l10; @Erdo06]. [@feldwh10] report an unusually large bulk flow of $416\pm 78\kms$ in a sphere of $ 100\hmpc$ which is at odds with the $\Lambda$CDM model with the best fit parameters of the Seven-Year Wilkinson Microwave Anisotropy Probe (WMAP7) [e.g. @jaro10; @wmap7]. Here we provide an alternative estimate of the bulk using a single data set of TF measurements of galaxies, trimmed at faint magnitudes to ensure the linearity of the TF relation. The estimate is based on a method which we term  for All Space 