be computed self-consistently, assuming the same normalization at the recombination epoch. Based on WMAP7, this implies $\sigma_8=0.63$ and $0.855$ for DGP and $f(R)$ models, respectively. Adopting these $\sigma_8$ values for these models we get $\gamma=0.315\pm 0.091$ and $0.55\pm 0.098$, respectively, for SGP and $f(r)$ (Cinzia Di Porto, priv. comm). In the DGP model, the expected value for $\gamma$ at $z=0$ is 0.664 [@Wu09], which is ruled by our constrain at more than the 3 $\sigma$ level. The $f(R)$ model cannot be ruled out at high confidence level by the constraint derived here.

Our results are in agreement with the analysis of [@st10]. Peculiar velocities from supernovae, although very sparse, yield bulk flows that are consistent with WMAP7 [@Dai11; @colin11], as we do. The results are in agreement with the WMAP7. The analysis of [@Bil11] of the 2MASS dipole from galaxy fluxes is also in agreement with the WMAP7 LCDM. The direction of the bulk flow is robust and agrees with the direction of the motion of the local group (LG) after correcting for the Virgocentric infall [@st10]. But we disagree strongly with the bulk flows [@feldwh10] who find a significant large bulk flow at $r=100\hmpc$, using the untrimmed SFI++ survey, 