New protocols for efficient tomography of bosonic quantum states achieve optimal Gaussian approximations in both high and low fidelity regimes, with the first tolerant Gaussianity testing method and proof that non-Gaussian measurements are necessary.
